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RESUMEN 



Se presenta el estudio de dos edificios esbeltos de 21 pisos con base en de marcos ductiles 
de acero estructural, disenados para distorsiones de entrepiso cercanas al limite maximo A=1.2% 
estipulado por el Reglamento del Distrito Federal vigente en Mexico (RCDF‘93). Los edificios no 
cumplen con la relacion de esbeltez limite H/L< 2.5 estipulada en el apartado seis de las Normas 
Tecnicas Complementarias por Sismo vigentes (NTCS'95), por lo que el factor reductivo de 
fuerzas sismicas con fines de diseno (Q') utilizado para el diseno fue afectado por un factor 
reductivo de 0.8 (Q -0.8Q), segun se especifica en la seccion 4.1 de las NTCS‘95. Las estructuras 
se disenaron para la zona III definida por las NTCS‘95, es decir, la zona de lago o de terreno 
blando. Se realizaron analisis dinamicos no lineales de marcos ante la accion de acelerogramas 
registrados y artificiales asociados al espectro de diseno. Se comparan los niveles y distribucion 
de fluencia estructural y las demandas de ductilidad maximas de cada nivel con respecto al factor 
de comportamiento sismico (Q) utilizado en el diseno. Se comparan, ademas, las distorsiones 
maximas de diseno con las distorsiones dinamicas maximas obtenidas de los analisis no lineales. 



Del estudio se desprende que aunque las demandas de ductilidad de entrepiso maximas 
se encuentran dentro de niveles razonables y que la fluencia de elementos estructurales 
corresponde al mecanismo ultimo de falla de viga debibcolumna fuerte, una interpretation 
erronea de las NTCS‘95 con respecto a como deben revisarse las restricciones de desplazamientos 
laterales para fines de diseno puede conducir a distorsiones dinamicas maximas de entrepiso del 
orden de 2.4%, que sobrepasan notablemente el limite maximo de distorsion de entrepiso 
permitido por el reglamento. A estos niveles de distorsion se pueden esperar danos considerables 
en elementos no estructurales, asi como problemas de funcionamiento y operation del edificio 
asociado a la conducta humana durante sismos. 



SUMMARY 



The study of two 21 story, slender, special moment resisting steel buildings is presented. 
Both buildings were designed according to the seismic provisions of RCDF-93 code of Mexico City 
for a maximum story drift ratio close to the limiting drift A=1.2% established by the code. As both 
buildings do not satisfy the limiting slenderness ratio H/L< 2.5 established by RCDF‘93 for 
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regular structures, then, the Q 'factor allowed by the code to reduce seismic forces for design 
purposes was affected by a 0.8 reduction factor, as established by RCDF-93 (Q-0.8Q'). Both 
buildings were designed for the lakebed region of Mexico City (soft soil conditions). Nonlinear 
dynamic analyses of representative frame models of subject structures were conducted. Several 
recorded and simulated accelerograms associated to the design spectra of RCDF-93 for the 
lakebed region were used for the nonlinear dynamic analyses. Response modification factors (Q) 
used in the design of subject buildings were compared against story ductility computed from the 
nonlinear dynamic analyses. Story drift ratios associated to the original design according to 
RCDF-93 were compared with peak dynamic story drift angles computed from nonlinear dynamic 
analyses. Structural yielding was studied and associated to hysteretic, deformation and strength 
demands. 

The study shows that although peak story ductility demands are within reasonable 
bounds, and the yielding mapping suggests the potential formation of a weak-beam, strong- 
column failure mechanism, peak dynamic drift angles are much larger than anticipated in the 
design (up to 2.4%). For drift angles higher than 2%, one should expect considerable damage in 
non-structural elements, important problems in the normal operation of the building, as well as 
potential panic in people, particularly people living in the lake-bed zone of Mexico City, where 
strong earthquakes are felt for a time interval of two or three seconds. Most residents of Mexico 
City that felt the M s =8.1 September 19, 1985 earthquake are particularly sensitive and still 
afraid of moderate and strong earthquakes felt in the city. 



1. INTRODUCTION 



Las Normas Tecnicas Complementarias para Diseno por Sismo (NTCS‘95) del 
Reglamento del Distrito Federal vigente (RCDF-93) establece, en su apartado seis, las once 
condiciones que debe cumplir una estructura para ser considerada como regular y, de esta 
manera, utilizar directamente el factor reductivo de fuerzas sismicas con fines de diseno (Q 1 ) en 
funcion de los factores de comportamiento sismico (Q) del apartado cinco. En caso de que no se 
cumpla una, varias o todas las condiciones de regularidad del apartado seis, Q 1 debera 
multiplicarse por un factor reductivo de 0.8, segun se especifica en la seccion 4.1 de las NTCS‘95. 
Esta disposicion reglamentaria esta basada exclusivamente en la intuicion y no se cuenta con 
estudios que la validen. Ademas, las normas de referenda no indican como proceder en caso de 
que se violen simultaneamente varias condiciones de regularidad, en cuyo caso el castigo de 0.8 
del factor reductivo puede no ser suficiente. 

A1 respecto, los comentarios a las NTCS-87, que aplican a las NTCS-95 puesto que no han 
sufrido modificaciones en este criterio, dicen lo siguiente : “La reduccion de Q’ establecida para 
edificios que no cumplen con requisitos de regularidad refleja el resultado de la experiencia y de 
la intuicion. En efecto, se espera rnejor comportamiento ante sismo de los edificios que satisfacen 
dichos requisites; ademas, un examen de los efectos de temblores intensos recientes confirma en 
terminos generales la menor vulnerabilidad de tales inmuebles. Sin embargo, no existen analisis 
ni pruebas de laboratorio que fundamenten cuantitativamente la reduccion de Q’ en un 20% para 
los edificios menos regulares, asi que esta disposicion no necesariamente implica que dichas 
estructuras desarrollan menores factores de ductilidad” (Rosenblueth y Gomez 1991). 

Por lo tanto, es evidente que ni siquiera puede asegurase que el factor Q' este relacionado 
con la esbeltez, o con alguna de las otras condiciones de regularidad, sino que la reduccion de Q’ 
es simplemente una medida preventiva arbitraria. En tal caso, es importante reflexionar sobre 
las siguientes cuestiones (entre muchas otras) : ante la accion de sismos, <^que debilidades 
presenta una estructura por el hecho de ser esbelta? ^Podra reducirse la vulnerabilidad ante 
sismos de una estructura esbelta con solo modificar el valor propuesto del factor Q’? De hecho, 
<ipodra realmente disminuirse la vulnerabilidad de estructuras irregulares reduciendo el valor de 
Q’ cuando no se satisface ninguna condicion de regularidad? ^ Seri a mas prudente tomar medidas 
preventivas diferentes en el diseno de estructuras irregulares dependiendo de la condicion de 
irregularidad que no se satisfaga, como lo hace, por ejemplo, en el reglamento UBC (Uniform 
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Building Code)? <^Es adecuado el factor de 0.80 para estructuras que violan la condicion de 
esbeltez u otras de las condiciones de regularidad? La contestacion a estas interrogantes requiere 
de muchos estudios detallados, de los cuales el estudio que aqui se presenta es solamente un 
modesto esfuerzo que intenta dar respuesta a algunas de ellas. 

Es un hecho que los efectos de las irregularidades pueden afectar el comportamiento de la 
estructura de diferentes maneras, y la influencia de estos no ha sido cuantificada suficientemente 
en estudios previos a nivel mundial realizados en las decadas de los anos setenta, ochenta y 
principios de los noventa, sobre todo en el intervalo de comportamiento inelastico, aunque en el 
ultimo lustro se han realizando algunos estudios para condiciones de irregularidad tales como 
piso debil, estructuras esbeltas, estructuras alargadas, irregularidades importantes en planta y/o 
elevacion, asi como irregularidades en la distribucion de masas, rigideces y resistencia (por 
ejemplo, Valmundsson y Nau 1997, Ruiz et al. 1997, Tena et al. 1996, Tena 1997, Tena-Colunga y 
del Valle 1998, Tena 2000, Tena 2001 y Tena, 2002). El estudio comparativo de las coincidencias 
y diferencias principales entre las disposiciones de diseno sismico para estructuras irregulares de 
los reglamentos mexicanos y el reglamento norteamericano UBC se ha presentado recientemente 
en la literatura (Tena-Colunga y del Valle 1998, Tena-Colunga 1999). 

Con respecto al diseno de estructuras irregulares, cabe senalar que las Normas Tecnicas 
Complementarias para Diseno por Sismo para el Distrito Federal de 2004 (NTCS-2004 2004), que 
estan por entrar en vigor, contemplan algunos cambios en el diseno de estructuras irregulares. 
Entre los mas importantes esta la definicion de estructuras fuertemente irregulares, que 
esencialmente son aquellas con un piso suave o una excentricidad torsional muy pronunciada. 
Para estas estructuras, Q' debe multiplicarse por un factor de 0.7. Otra modificacion es que si 
una estructura no satisface una condicion de regularidad, Q' se multiplicara por 0.9 y por 0.8 
cuando no cumpla con dos o mas de las condiciones de regularidad del apartado 6 de las NTCS- 
2004, que son esencialmente los mismos que se definen en las NTCS-95 y NTCS-87. 

En el presente trabajo se resume brevemente los estudios realizados en estructuras 
esbeltas de 21 pisos a base de marcos ductiles de acero estructural, disenados para distorsiones 
de entrepiso cercanas al limite maximo A=1.2% estipulado por el RCDF-93 que no cumplen con la 
relacion de esbeltez limite H/L< 2.5 estipulada en el apartado seis de las NTCS-95, y que se 
presentan en detalle en trabajos mas extensos (Cheja 1997, Tena et al. 1996 y Tena 1997), donde 
ademas se presenta el estudio realizado para edificios de ocho pisos con relaciones de esbeltez 
similares. El diseno de cada estructura se hizo utilizando el factor de comportamiento sismico 
adecuado conforme a las NTCS-95, como se presenta brevemente en las siguientes secciones. Este 
estudio es tambien representative de la propuesta de NTCS-2004, ya que de acuerdo con ellas, el 
diseno sismico por condicion de irregularidad seria el mismo que se establece en las NTCS-95, 
ademas que los criterios especificos para el diseno de marcos ductiles de acero estructural son 
esencialmente los mismos, sobre todo los que rigen el diseno de vigas y columnas. 



2. ESTRUCTURAS EN ESTUDIO 



En este estudio se eligio como estructura esbelta prototipo al proyecto arquitectonico de 
una de las torres de unos condominios de lujo ubicados en la calle Ruben Dario en la colonia 
Polanco, en una de las zonas residenciales mas exclusivas de la ciudad de Mexico. Se considero 
que este tipo de proyecto es representative de estructuras esbeltas comunes que pueden 
construirse en los proximos anos en distintas zonas de la ciudad de Mexico, incluyendo la zona del 
lago. 



El conjunto condominal cuenta con un terreno con dimensiones de 64x57 m. y un area de 
3648 m 2 , y esta compuesto por tres torres (A, B y C) de 21 pisos y 66 m. de altura sobre la 
superficie del terreno, y de tres sotanos comunes donde se alojan las areas de estacionamiento y 
servicios (SI, S2 y S3). Por su parte, las torres A y B son iguales entre si, mientras que la torre C, 
aunque presenta la mismas altura que las dos anteriores, es un edificio mas regular tanto en 
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planta como en elevacion. Los detalles del proyecto se presentan en Cheja (1997); sin embargo, a 
continuacion se resumen las caracteristicas generales de la torre A, cuya geometria fue 
seleccionada para el presente estudio. 

La torre A esta conformada por un mezanine en el primer piso, y diez departamentos de 
dos niveles cada uno en los pisos subsecuentes. El mezanine tiene una altura de 4.00 m., la planta 
baja de cada departamento es de 3.40 m. y la planta alta de 2.80 m. Los elevadores y las escaleras 
de servicio solo tienen entrada en la planta baja de cada condominio, y una escalera interior 
ubicada entre los ejes C, D, 8 y 9 (figura l) permite el acceso a la planta alta. Todas las losas 
tienen un saliente entre los ejes 7 y 8 de forma sensiblemente triangular y, excepto en la azotea, 
otro entre los ejes 8 y 9 donde se aloja una de las jardineras. La jardinera ubicada entre los ejes C 
y D solo aparece en la planta baja de cada departamento, ya que en dicha area existe doble altura 
y, por lo tanto, en la losa de entrepiso es espacio vacio. El area entre los ejes D, E, 6 y 7 esta 
destinada para las escaleras de servicio. Los claros principales tienen longitudes entre 3.80 m. el 
mas pequeno y 6.00 m. el mas grande. Las dimensiones gruesas de la planta son de 15.9 x 20.0 m, 
siendo esta simetrica respecto a un eje central paralelo al lado mayor de 20.0 m. 




Se eligio la geometria de la torre “A” como prototipo por ser adecuada para el estudio de 
estructuras esbeltas, ya que no cumple con la condicion de regularidad por esbeltez (condicion 2 
NTCS'95). La condicion de esbeltez no se cumple en ambas direcciones en planta (figura l), ya 
que para la direccion X, la relacion de esbeltez es 66/15.9=4.15, mientras que para la direccion Y 
dicha relacion resulta igual a 66/20.0=3.30; ambas mayores al limite de 2.5. Ademas, la torre A 
tampoco cumple con las condiciones de regularidad 6, 7, y 11 de las NTCS‘95, segun se demuestra 
con detalle en Cheja (1997). 

2.1 Diseno sismico de los edificios IR6 e IR6B 
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Se realizo el diseno sismico de dos modelos, denominados IR6 e IR6B, con base en el 
proyecto arquitectonico de la torre “A”, con la salvedad que en el presente estudio se considera 
que la estructura tiene condicion de empotramiento perfecto en su base, lo que obedece a que se 
pretende estudiar exclusivamente la validez de las disposiciones debidas a la irregularidad de la 
estructura, tales como la esbeltez; despreciando otras acciones que pueden ser significativas como 
la interaccion suekrestructura. Cabe senalar que el que no se consideren la interaccion suekr 
estructura y otros efectos no implica de ninguna manera que no se consideren importantes, mas 
bien, el razonamiento asociado a esta decision es que el incluir otros efectos (como la interaccion) 
no permitiria discernir claramente en que medida afectan las condiciones de irregularidad 
(particularmente la esbeltez), por si solas, al diseno y a la respuesta dinamica no lineal de 
estructuras disenadas conforme a las NTCS'95, que es realmente uno de los objetivos centrales 
que se persiguen con los estudios que se estan conduciendo sobre este tema. 

Los modelos IR6 e IR6B se disenaron considerando que pertenecen al grupo “B” segun el 
RCDF-93, y como tales fueron considerados los factores de carga para el diseno. Los valores del 
espectro de diseno utilizados para el analisis sismico son los correspondientes a la zona III (figura 
2), es decir, la zona de mayores ordenadas espectrales de diseno, y en la cual resulta mas 
interesante evaluar a las NTCS'95. La estructuracion considerada son marcos ductiles a 
momento, por lo que el factor de comportamiento sismico empleado fue Q=4. Debido a que no 
cumplen con las condiciones de regularidad antes destacadas, el factor reductivo de fuerzas 
sismicas con fines de diseno (Q') se multiplied por el factor reductivo de 0.8 especificado en la 
seccion 4.1 de las NTCS'95, y dadas las caracteristicas dinamicas de los modelos (T<Tb), en este 
caso Q'=0.8Q=4x0.8=3.2. Para las restricciones de desplazamientos laterales que indica el 
articulo 209 del RCDF-93, se acepto que para las estructuras en cuestion los elementos no 
estructurales estan totalmente desligados de la estructura principal, por tanto, se considero como 
limite de la distorsion de entrepiso de diseno el valor de A=1.2%. 




Figura 2. Espectros de diseno para estructuras del grupo B, zona III del RCDF-93 



La diferencia basica entre el diseno de IR6 e IR6B estriba exactamente en como se revisan 
las restricciones de desplazamientos laterales para fines de diseno (Ad), ya que para algunos la 
redaccion actual de las NTCS'95 no es totalmente clara al respecto. De acuerdo con la seccion 4.1 
de las NTCS'95, “las deformaciones laterales se calcularan multiplicando por Q las causadas por 
las fuerzas sismicas reducidas cuando se emplee el metodo estatico de analisis que se detalla en 
la seccion 8 de las presentes normas o el analisis modal de la seccion 9”. Si llamamos Aq- a las 
deformaciones laterales causadas por las fuerzas sismicas reducidas por Q’, entonces, de acuerdo 
con las NTCS'95, A d = QAq-, que es el criterio con que se diseno IR6B. En cambio, IR6 se diseno 

considerando que A D =Q'A q - , ya que existe quien interpreta que el castigo por el efecto de 

irregularidad es reducir Q por 0.8, por lo que consideran logico revisar las deformaciones 
suponiendo que Q =Q, que no es lo que se estipula en las NTCS'95 y, de hecho, las NTCS'2004. 



Arturo Tena Colunga y Carlos Cheja Mochon 



6 



Por tanto, en este caso concreto, se considero que A D = 3.2A Q -para IR6 y A D =4A Q -para IR6B. De 

hecho, de acuerdo con las NTCS'95, las estructuras irregulares deben disenarse para resistir 
fuerzas sismicas reducidas mayores que las estructuras regulares y, si se observa con detalle, el 
criterio de revision para deformaciones laterales, aunque pareciera ser el mismo que para las 
estructuras regulares es, de hecho, mas estricto. Este criterio se mantiene en la propuesta de 
NTCS-2004. 

En la estructuracion con base en marcos de acero ductiles se utilizaron secciones cajon 
para las columnas, pues sus propiedades son mas convenientes para dar estabilidad ante carga 
lateral y evitar problemas de pandeo y, ademas se usan frecuentemente en la practica mexicana. 
Se seleccionaron secciones I (perfiles tipo W) para las vigas principals y secundarias, verificando 
que cumplieran los criterios para secciones compactas correspondientes. Como sistema de piso se 
propuso una losa de lamina de acero acanalada y concreto reforzado con malla electrosoldada de 5 
cm. de espesor, soportada por trabes secundarias. El diseno estructural de IR6 e IR6B se llevo a 
cabo con apoyo del programa ETABS. Se realizaron las siguientes hipotesis en el modelado de las 
estructuras. El modulo de elasticidad del acero se tomo como E = 2.039x10® Kg./cm^. Se 
utilizaron las propiedades geometricas nominales de los perfiles W en las vigas y las propiedades 
geometricas efectivas en las columnas de section cajon. La zona rigida que se considero tanto 
para las columnas como para las trabes es del 50% de la longitud total de apoyo entre los 
elementos, segun se recomienda para este tipo de edificios. Se realizo un diseno conforme a un 
analisis modal espectral elastico tridimensional conforme a las NTCS'95 y al RCDF'93. En el 
analisis estructural fueron tornados en cuenta los efectos secundarios (P'A) que se inducen a la 
estructura por los desplazamientos laterales debidos al sismo, los cuales deben de incluirse en el 
diseno de estructuras esbeltas. 

El diseno preliminar de los elementos estructurales principales de acero se hizo conforme 
al reglamento norteamericano “Load and Resistance Factor Design” (LRFD'93), utilizando el 
modulo STEELER del ETABS. Posteriormente, se verified que las columnas y trabes mas 
esforzadas de dicho diseno, cumplieran con las especificaciones de las Normas Tecnicas 
Complementarias Para Diseno y Construction de Estructuras Metalicas (NTCEM'95, 1995); como 
se ilustra en detalle en Cheja (1997). Los detalles de los disenos se presentan en Cheja (1997), 
Tena et al. (1996) y Tena (1997) y se resumen en las tablas 1 y 2, donde se presentan las secciones 
finales de diseno para columnas y trabes, respectivamente, y en las figuras 3 y 4, donde se 
presentan los niveles de esfuerzo de trabes y columnas asociados a los disenos de IR6 e IR6B, 
respectivamente . 



Tabla 1. Secciones finales de diseno de columnas 



Modelo IR6 


Modelo IR6B 


Niveles 


Columnas 


Dimensiones (cm.) 


Niveles 


Columnas 


Dimensiones (cm.) 


1-2 


C1-C5, C8, C9, 
C12, C13, C16 a 
C20 


70x60x2.54 


1-4 


C1-C20 


70x70x2.86 


1-2 


C6, C7, CIO, Cll, 
C14, C15 


50x50x2.54 


3-10 


C1-C20 


50x50x2.54 


5-10 


C1-C20 


60x60x2.38 


11-16 


C1-C20 


45x40x2.54 


11-16 


C1-C20 


50x50x1.91 


17-21 


C1-C20 


40x40x1.27 


17-21 


C1-C20 


40x40x1.59 



Tabla 2. Secciones finales de diseno para trabes 
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Niaioveles 


Ceaw2556 


Modelo IR6B 


Trabes 


Perfil 


Trabes 


Perfil 


1-9 


T-l, T-3 a T-15 


W24xl03 


T-l, T-3 a T-15 


W24xl31 


T-2 


W24xll7 


T-2 


W24xl46 


T-16 a T-31 


W24x84 


T16 a T-31 


W24xl03 


10-12 


T-l, T-3 a T-15 


W24xl03 


T-l, T-3 a T-15 


W24xl31 


T-2 


W24xll7 


T-2 


W24xl46 


T-16 a T-31 


W24x76 


T-16 a T-31 


W24x94 


13 


T-l, T-3 a T-15 


W24x84 


T-l, T-3 a T-15 


W24xll7 


T-2 


W24xl03 


T-2 


W24xl31 


T-16 a T-31 


W24x68 


T-16 a T-31 


W24x94 


14-16 


T-l, T-3 a T-15 


W24x76 


T-l, T-3 a T-15 


W24xll7 


T-2 


W24x84 


T-2 


W24xll7 


T-16 a T-31 


W24x68 


T-16 a T-31 


W24x94 


17 


T-l, T-3 a T-15 


W24x62 


T-l a T-15 


W24x94 


T-2, T-16 a T-31 


W24x68 


T-16 a T-31 


W24x94 


18-21 


T-l, T-3 a T-31 


W24x62 


T-l a T-15 


W24x94 


T-2 


W24x68 


T-16 a T-31 


W24x76 



En las figuras 3 y 4 se presentan las distorsiones de entrepiso de diseno calculadas para 
las dos direcciones principales de los modelos IR6 e IR6B respectivamente, y se comparan con los 
llmites establecidos por el RCDF-93. Tambien se ilustran los niveles de esfuerzo del diseno final 
de trabes y columnas. Puede observarse que en los niveles 6 a 14 se presentan distorsiones en la 
direccion X muy cercanas al limite de 1.2%; y que tanto las columnas de los primeros niveles, 
como las trabes paralelas al eje Y de los niveles intermedios, se encuentran sometidas a un alto 
nivel de esfuerzos. Se observa que la direccion critica ambos modelos es la direccion X. Es 
evidente tambien que muchos de los elementos se encuentran sobredisenados, en gran medida, 
como consecuencia de tomar secciones tipo, ya que como se aclaro anteriormente, se pretende que 
el diseno sea lo mas simple, realista y construible posible, con distorsiones cercanas al limite 
maximo permitido por el RCDF-93, pero no necesariamente optimo en cuanto a peso y/o nivel de 
esfuerzos. En este respecto, parece que el diseno de IR6B es ligeramente mas optimo en cuanto al 
nivel de esfuerzos de los elementos estructurales que el de IR6, conforme a las acciones 
estipuladas por el Reglamento. 




Figura 3. Niveles de esfuerzo y distorsiones de diseno, modelo IR6 
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Figura 4. Niveles de esfuerzo y distorsiones de diseno, modelo IR6B 



2.2 Caracteristicas dinamicas elasticas de los modelos IR6 e IR6B 



Los resultados del analisis modal del diseno final de los modelos IR6 e IR6B se resumen 
escuetamente en la tabla 3, donde se indica tambien el peso total de cada estructura. Ambos 
modelos son flexibles, pero un poco “rigidos” (valga la aparente contradiccion) para ser 
estructuras de marcos de acero de 21 niveles y 66 m de altura, debido a que se considero que los 
muros de bloque de concreto que se ubican en el area del cubo de elevadores participan en la 
rigidez lateral inicial de las estructuras. Las estructuras son mas flexibles en la direccion E-W, en 
la cual son mas esbeltas. Los modos de traslacion estan moderadamente acoplados, debido a las 
importantes excentricidades de masa que existen por el uso de areas de los entrepisos y a la 
ubicacion de los 



Tabla 3. Caracteristicas dinamicas de los modelos IR6 e IR6B 



MODO 


Direccion 


Modelo IR6 
(Wi6=5542.7 Ton) 


Modelo IR6B 
Wi6B=5704.5 Ton) 


T (s) 


Masa modal (%) 


T (s) 


Masa modal (%) 


N-S 


E-W 


0 


N-S 


E-W 


0 


1. Primer modo 

traslacion 


E-W 


1.923 


4.64 


63.02 


6.45 


1.646 


3.59 


63.10 


5.79 


2. Primer modo 

traslacion 


N-S 


1.882 


70.51 


4.51 


0.07 


1.615 


70.08 


3.50 


0.05 


3. Primer modo torsion 


Rotacion 


1.352 


0.11 


5.64 


69.45 


1.197 


0.07 


5.22 


69.04 


4. Segundo modo 
traslacion 


E-W 


0.660 


2.20 


10.49 


0.83 


0.566 


4.03 


9.47 


0.83 


5. Segundo modo 
traslacion 


N-S 


0.646 


11.33 


2.25 


0.14 


0.559 


9.92 


4.10 


0.30 


6. Segundo modo torsion 


Rotacion 


0.479 


0.03 


1.91 


10.88 


0.430 


0.02 


1.89 


11.32 



muros de bloque de concreto (figura l). Los modos de traslacion en direccion N-S presentan menor 
acoplamiento en rotacion que los modos de traslacion en direccion E-W. Los modos de torsion 
presentan mayor acoplamiento con traslacion en la direccion E-W. El peso total de IR6 es 
Wjg=5,542.7 ton y el de IR6B es Wjgg=5, 704.5 ton. El peso de los elementos estructurales 
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principales es aproximadamente 14% el peso total de la estructura para IR6 y 16% para IR6B. Se 
observa que el interpretar con rigor las disposiciones de las NTCS‘95 para revisar las 
restricciones de desplazamientos laterales para fines de diseno (Ad) en el diseno de estructuras 
irregulares guio al diseno de una estructura mas pesada y rigida (IR6B) con respecto al diseno 
resultante con la interpretacion que muchos otros estructuristas le dan a dichas normas (IR6). 
Sera importante evaluar si esta diferencia en la interpretacion puede guiar a diferencias de 
comportamiento notables entre ambos modelos ante la accion de sismos intensos. Este es un 
aspecto sobre el cual reflexiona el presente trabajo. 



3. ANALISIS NO LINEALES 



Los disenos de los modelos IR6 e IR6B, conforme a los lineamientos de las NTCS‘95 y 
NTCEM‘95 del RCDF-93, se evaluaron mediante la simulacion de su respuesta ante la accion de 
acelerogramas asociados al sismo de diseno realizando analisis dinamicos no lineales en el 
programa DRAIN-2DX (Prakash et al. 1992). 

Debido a que en la investigacion en curso uno de los objetivos es estudiar si los factores de 
comportamiento sismico para el diseno de estructuras irregulares en la zona III de las NTCS‘95 
son indicados, es necesario utilizar acelerogramas representativos cuyas caracteristicas esten 
relacionadas al sismo maximo esperado y contemplado por el presente reglamento y 
acelerogramas que pudieran estar contemplados en reglamentos futuros. Por lo tanto, se 
seleccionaron para el estudio de IR6 e IR6B cinco acelerogramas, dos registrados en el valle de 
Mexico durante el sismo del 19 de septiembre de 1985 (M g =8.l) y tres acelerogramas sinteticos 
generados para un sismo de M g =8.2 para tres estaciones de la red acelerografica del valle de 
Mexico. 



Los acelerogramas registrados durante el sismo del 19/09/85 que fueron seleccionados 
corresponden a los sitios Secretaria de Comunicaciones y Transportes en su componente E‘W 
(SCT-EW) y Tlahuac bombas en su componente N-S (TBOM-NS). Las caracteristicas de los 
registros difieren en cuanto a duracion, contenido de frecuencias, periodo dominante y respuestas 
espectrales maximas. Los registros de SCT tienen una duracion de un poco mas de 160s y 
aceleracion maxima del terreno de aproximadamente 0.17g en su componente E‘W, con respuesta 
espectral maxima de l.Og para un periodo dominante de 2.0s. El registro de Tlahuac bombas en 
direccion N-S tiene una duracion de 150s con espectro de respuesta con periodos dominantes 
cercanos a los 2.0s, con aceleraciones maximas del terreno de aproximadamente 0.14g y 
aceleraciones espectrales maximas de 0.75g. 

Los registros artificiales para estaciones de la red acelerografica de la ciudad de Mexico se 
generaron para un sismo postulado de M g =8.2, basado en las funciones de transferencia promedio 
calculadas para cada estacion con base en la informacion capturada por la red acelerografica a 
partir de 1988, conforme a los procedimientos expuestos en Perez Rocha (1998). Se seleccionaron 
tres estaciones cuyos niveles de respuesta espectral de pseudoaceleraciones resultaran similares 
o mayores al acelerograma registrado en SCT-EW durante el sismo del 19 de septiembre de 1985 
y que tuvieran periodos dominantes de respuesta mayores a 1.3s y menores a 2.2s dadas las 
caracteristicas de las estructuras a estudiar. Se selecciono, por tanto, la estacion 05 (Las Cibeles), 
ubicada en la colonia Condesa, con periodo dominante de 2.15s y aceleracion maxima espectral de 
0.86g; la estacion 56, ubicada en la calle de Cordoba, en la colonia Roma, con periodo dominante 
de 2.20s y aceleracion maxima espectral de 1.44g (figura 5), y la estacion 84 (Culhuacan), ubicada 
en San Francisco Culhuacan, al sur de la ciudad de Mexico, con periodo dominante de 1.38s y 
aceleracion maxima espectral de 1.03g. Por fines de espacio no se muestran todos los registros y 
sus caracteristicas espectrales, pero pueden consultarse en otros trabajos (por ejemplo, Cheja 
1997, Tena et al. 1996, Tena 1997). 
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Los analisis dinamicos no lineales se condujeron exclusivamente en los marcos mas 
esforzados en cada direccion (marco 10 en la direccion X, marco B en la direccion Y), es decir, se 
obtuvieron marcos dinamicamente equivalentes con el procedimiento explicado en detalle en 
Cheja (1997). Se procedio de esta manera, ya que en estudios paralelos que se condujeron en 
edificios esbeltos de acero de ocho niveles cuya geometrla se definio de manera que fuera analoga 
a la de IR6 e IR6B (Cheja 1997, Tena et <3/1996), se observaron respuestas dinamicas muy 
similares en analisis conducidos en marcos equivalentes con respectos a modelos que consideran 
la interaction entre los marcos que resisten a las fuerzas sismicas en una direccion dada, por lo 
que se decidio, para una primera evaluation, economizar en el esfuerzo y recursos de compute 
necesarios para la evaluation de la respuesta dinamica de IR6 e IR6B. 




Figura 5. Acelerograma artificial S56-EW 



En los analisis no lineales se considero un nivel de amortiguamiento viscoso equivalente 
<^=3%, ya que el nivel de amortiguamiento viscoso equivalente medido experimentalmente en 
algunas estructuras de acero varia entre C= 2% y C= 4%. Se considero que los elementos 
estructurales son capaces de desarrollar sus capacidades plasticas teoricas y que se presenta un 
comportamiento bilineal del acero estructural, donde inicialmente se considero una rigidez 
posterior a la fluencia del 2% de la rigidez elastica, que es la minima observada en elementos de 
acero estructural en pruebas experimentales y que ha sido reportada indirectamente en la 
literatura (por ejemplo, Mamaghani et al. 1995). En estudios subsecuentes se han considerado 
tambien pendientes posteriores a la fluencia de 5% y 10%, para evaluar su influencia en la 
respuesta dinamica no lineal, resultados que tambien se discuten en este trabajo. 

Se estudiaron los siguientes parametros dinamicos para los marcos 10 (direccion X) y B 
(direccion Y) de los modelo IR6 e IR6B: (a) curvas de histeresis de entrepiso, (b) mapeo de la 
fluencia de elementos, en tiempos especificos, en envolventes, y su monitoreo en el dominio del 
tiempo, (c) envolventes de distorsiones de entrepiso maximas (A), (d) envolventes de cortantes de 
entrepiso maximos (V/Wt), (e) demandas maximas de ductilidad de entrepiso (a), (f) distorsion de 
entrepiso asociada a la primera fluencia de cualquier elemento estructural que forma parte del 
entrepiso (Ate), (g) distorsion de entrepiso maxima asociada al medio ciclo de histeresis de carga 
maxima (A me ), (h) cociente del numero de ciclos donde se presento respuesta inelastica (Cine) entre 
el numero total de ciclos a que fue sujeta la estructura ante un acelerograma dado (Ctot), (i) 
cociente entre la rigidez de entrepiso secante minima inelastica (kimin) o “pico a pico” asociada a 
los ciclos histereticos de amplitud maxima y la rigidez elastica del entrepiso (kei) y, (j) cociente 
entre el promedio de las rigideces de entrepiso secantes (kiprom) asociadas a los medios ciclos 
histereticos inelasticos (kmc) y la rigidez elastica del entrepiso (k e i). Algunos de estos conceptos se 
ilustran en la figura 6. 

La information analizada es muy extensa, ya que se realizaron analisis para los marcos y 
modelos de referenda ante la action de todos los acelerogramas senalados anteriormente, 
considerando pendientes posteriores a la fluencia de 2%, 5% y 10%, y por razones de espacio en 
este trabajo solamente se presentan y discuten algunos resultados del estudio descrito. Las 
mayores respuestas para los marcos 10 y B de los modelos IR6 e IR6B se obtuvieron ante la 
action de los acelerogramas SCT85-EW y S56-EW; sin embargo, solo se presentan los resultados 
de los marcos 10 de los modelos IR6 (figura 8) e IR6B (figura 7) y de los marcos B de IR6 (figura 
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10) e IR6B (figura 9) ante la accion del acelerograma S56-EW, ya que, en general, las respuestas 
mas crlticas se asocian a este acelerograma. 




De las envolventes de respuesta maximas del marco 10 del modelo IR6B (el que cumple 
con todas las disposiciones de las NTCS-95 para edificios irregulares) se observa que su respuesta 
no lineal es reducida (figura 7), con demandas maximas de ductilidad de 1.3 asociadas a rigideces 
secantes minimas inelasticas (pico a pico) del 73% de su rigidez elastica (k2=2%, nivel 5), a un 
valor minimo observado del promedio de las rigideces secantes inelasticas de medio ciclo del 83% 
de su rigidez elastica (ka=2%, nivel l) y distorsiones de entrepiso maximas muy por debajo de las 
calculadas en su diseno. Se presento respuesta elastica en numerosos niveles. En aquellos 
entrepisos donde se presento respuesta inelastica, la distorsion de entrepiso asociada a la 
primera fluencia del mismo vario de Af e =0.31% (k2=2%, nivel l) a Af e =0.65% (k2=2%, nivel 8), 
obteniendose en algunos entrepisos respuestas elasticas para distorsiones tan altas como 
A=0.76% (k2=2%, nivel 12), lo que llama poderosamente la atencion. La distorsion de entrepiso 
maxima asociada al medio ciclo de histeresis de carga maxima fue Ame=0.97% en los niveles 4 a 6, 
la cual se encuentra por debajo del limite maximo permitido por el RCDF-93 (A=1.2%). 
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Figura 7. Envolventes de respuesta maxima, marco 10 del modelo IR6B, acelerograma S56-EW 
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Como se puede observar de la figura 4, el diseno de IR6B en esta direction, a pesar de ser 
mas flexible, esta asociado a niveles de esfuerzos de las vigas inferiores, en promedio, a aquellos 
calculados para la direction Y. La fluencia se presento exclusivamente en vigas, y observando al 
nivel de esfuerzos de diseno asociados al diseno, no se puede explicar el por que existe tanta 
variabilidad en la distorsion asociada a la primera fluencia entre los distintos entrepisos (Af e ), 
aspecto que se analizara con mayor detalle en estudios futuros. Ademas, para este edificio no se 
observan grandes diferencias en las respuestas dinamicas maximas del marco 10 con respecto a 
la hipotesis hecha con respecto a la pendiente posterior a la fluencia de los elementos 
estructurales (2%, 5% o 10%). 



El cortante basal que resiste el marco 10 ante la action del acelerograma S56-EW varia 
entre el 5.8% y 6.0% del peso total del edificio IR6B, por lo que crudamente se puede estimar que 
la resistencia a cortante basal de IR6B en direction X debe encontrarse entre el 29% (k2=2%) y 
30% (k2=10%) del peso total de la estructura, ya que en esta direction participan cinco marcos de 
caracteristicas muy similares en cuanto a estructuracion y rigidez. Por tanto, si se toma en 
cuenta que la estructura se diseno para Q'=3.2 y en la meseta del espectro de diseno para zona III 
(c=0.40, figura 2), entonces, el cortante sismico no factorizado en una direction de analisis seria 

c 

crudamente estimado como V = — VP =0. 125VP , por lo que la sobrerresistencia observada para 



este diseno en direction X, ante la action del acelerograma de referenda, debe variar crudamente 
entre 2.3 y 2.4, ignorando aun otras fuentes posibles de sobrerresistencia como la action 
tridimensional de la estructura y la contribution en la rigidez y resistencia del sistema de piso. 
En estudios futuros se pretende cuantificar de una manera mas precisa la sobrerresistencia que 
se presenta en esta y otras estructuras. 



De las envolventes de respuesta maximas del marco 10 del modelo IR6 (donde se estudia 
una variation en la interpretation de las disposiciones de las NTCS‘95 de como revisar las 
restricciones de desplazamientos laterales para fines de diseno, Ad, para edificios irregulares), 
figura 8, se observa que su respuesta no lineal, aunque reducida, es mas importante que la 
observada en el modelo IR6B. La demanda maxima de ductilidad de entrepiso es de 2.26 (k2=2%, 
nivel 5) asociada a rigideces secantes minimas inelasticas (pico a pico) del 43% de su rigidez 
elastica (k2=2%, nivel 5), a un valor minimo observado del promedio de las rigideces secantes 
inelasticas de medio ciclo del 68% de su rigidez elastica (k2=2%, nivel 12). Las distorsiones de 
entrepiso maximas sobrepasan a las calculadas en su diseno en los niveles 3 a 9, y al limite 
maximo permitido por el RDF-93 A=1.2%. La distorsion maxima corresponde al sexto entrepiso y 
es A=1.58%, calculada esta distorsion con respecto a la position inicial de equilibrio. Se presento 
respuesta elastica exclusivamente en los niveles 18 y 20. En aquellos niveles donde se presento 
respuesta inelastica, la distorsion de entrepiso en que se presento su primera fluencia vario de 
Af e =0.30% (k2=2%, nivel l) a Af e =0.97% (k2=2%, nivel 16), obteniendose en algunos entrepisos 
respuestas elasticas para distorsiones tan altas como A=1.07% (k2=2%, nivel 18), coincidiendo 
esta gran variation con lo observado en IR6B. La distorsion de entrepiso maxima asociada al 
medio ciclo de histeresis de carga maxima fue A me =2.2% en los niveles 5 y 6, valor que se 
encuentra muy por encima del limite maximo permitido por el RCDF-93 (A=1.2%). 



La fluencia se presento principalmente en vigas, aunque las columnas exteriores de los 
niveles 1, 3 y 4 experimentaron fluencias incipientes (se presentaron tres veces en 2500 pasos de 
integration analizados). Al igual que para IR6B, observando al nivel de esfuerzos asociados al 
diseno, no se puede explicar el por que existe tanta variabilidad en la distorsion asociada a la 
primera fluencia entre los distintos entrepisos (Afe). Al igual que para IR6B, no se observan 
grandes diferencias en las respuestas dinamicas maximas del marco 10 de IR6 con respecto a 
pendientes posteriores a la fluencia de los elementos estructurales de 2%, 5% o 10%. El cortante 
basal que resiste el marco 10 de IR6 ante la action del acelerograma S56-EW varia entre 5.1% 
(k2=2%) y 5.5% (k2=10%) del peso total del edificio IR6. Siguiendo un razonamiento analogo al 
presentado para IR6B, se puede estimar crudamente que la resistencia a cortante basal de IR6 en 
direction X debe encontrarse entre el 25% y 28% del peso total de la estructura, y como en este 
caso el cortante sismico no factorizado en una direction de analisis seria crudamente estimado 
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tambien como V = 0.125W , la sobrerresistencia observada para este diseno en direccion X, ante la 
action del acelerograma de referenda, debe variar crudamente entre 2.0 y 2.25. 
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Figura 8. Envolventes de respuesta maxima, marco 10 del modelo IR6, acelerograma S56-EW 



Se puede concluir que la respuesta dinamica del marco 10 del modelo IR6B, que se diseno 
completamente conforme a las disposiciones de las NTCS'95 del RCDF-93 para estructuras 
irregulares, es mas favorable en todos los rubros que la del modelo IR6, donde se estudio una 
variation en la interpretation de las disposiciones de las normas para la revision de las 
restricciones de desplazamientos laterales para fines de diseno, que como puede observarse a 
partir de la discusion y de las figuras 7 y 8, conduce a desplazamientos y distorsiones de 
entrepisos mayores a los permitidos por el Reglamento, que si bien no estan asociados a dano 
estructural o demandas inelasticas severas, a estos niveles de distorsion se pueden esperar danos 
considerables en elementos no estructurales, asi como problemas de funcionamiento y operation 
del edificio asociados a la conducta humana durante sismos que deben ser tornados en cuenta por 
el proyectista. 

En la direccion Y se presentan mayores demandas no lineales que en la direccion X ante 
la action del acelerograma S56-EW, principalmente en los niveles inferiores, a juzgar de 
comparar las respuestas de los marcos B de IR6B e IR6 (figuras 9 y 10) con las de los marcos 10 
discutidas anteriormente (figuras 7 y 8). De las envolventes de respuesta maximas del marco B 
del modelo IR6B se observa que su respuesta no lineal es reducida (figura 9), con demandas 
maximas de ductilidad de 2.1 asociadas a rigideces secantes minimas inelasticas (pico a pico) del 
54% de su rigidez elastica (k2=2%, nivel 5), a un valor minimo observado del promedio de las 
rigideces secantes inelasticas de medio ciclo del 91% de su rigidez elastica (k2=2%, niveles 3 y 4) y 
distorsiones de entrepiso maximas que sobrepasan las calculadas en su diseno en los niveles 1 a 8 
(k2=1.17%, nivel 5), pero que no sobrepasan el limite del limite maximo de distorsion (A=1.2%) 
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permitido por el Reglamento. Se presento respuesta elastica en los niveles superiores 15 a 21. En 
aquellos entrepisos donde se presento respuesta inelastica, la distorsion de entrepiso asociada a 
la primera fluencia del mismo vario de Af e =0.28% (k2=2%, nivel l) a Af e =0.61% (ka=2%, nivel 12), 
obteniendose en algunos entrepisos respuestas elasticas para distorsiones de hasta A=0.43% 
(k2=2%, niveles 15 y 18). La distorsion de entrepiso maxima asociada al medio ciclo de histeresis 
de carga maxima fue A me =0.97% en los niveles 4 a 6, la cual se encuentra por debajo del limite 
maximo permitido por el RCDF-93 (A=1.2%). La distorsion de entrepiso maxima asociada al 
medio ciclo de histeresis de carga maxima fue Ame=1.4% en los niveles 4 y 5, valor que se 
encuentra por encima del limite maximo permitido por el RCDF-93 (A=1.2%). Como se puede 
observar de la figura 4, el diseno de IR6B en esta direccion, es mas rigido y esta asociado a 
niveles de esfuerzos mayores para las vigas, en promedio, a aquellos calculados para la direccion 
X. 
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Figura 9. Envolventes de respuesta maxima, marco B del modelo IR6B, acelerograma S56-EW 



La fluencia se presento exclusivamente en vigas, y observando al nivel de esfuerzos 
asociados al diseno, aqui tampoco se puede explicar el por que existe tanta variabilidad en la 
distorsion asociada a la primera fluencia entre los distintos entrepisos (Af e ). Aunque no se 
observan grandes diferencias en las respuestas dinamicas maximas del marco B con respecto a la 
hipotesis hecha con respecto a la pendiente posterior a la fluencia de los elementos estructurales 
(2%, 5% o 10%), si son un poco mas apreciables que para el marco 10 (figura 7) y, en general, 
tienden a disminuir las demandas de ductilidad y a incrementar las rigideces secantes inelasticas 
y la distorsion de entrepiso asociada a la primera fluencia. 

El cortante basal que resiste el marco B ante la accion del acelerograma S56-EW varia 
entre el 7.5% (k2=2%) y 7.8% (k2=10%) del peso total del edificio IR6B, por lo que crudamente se 
puede estimar que la resistencia a cortante basal de IR6B en direccion Y debe encontrarse entre 
el 30% y 31.2% del peso total de la estructura, ya que en esta direccion participan cuatro marcos 
de caracteristicas muy similares en cuanto a estructuracion y rigidez. Por tanto, la 
sobrerresistencia observada para este diseno en direccion Y, ante la accion del acelerograma de 
referenda, debe variar crudamente entre 2.4 y 2.5. 



15 



Comportamiento sismico de edificios esbeltos con base en marcos ductiles de acero 





K. 

limn 



K. 

lprom 



A 



el 



Figura 10. Envolventes de respuesta maxima, marco B del modelo IR6, acelerograma S56-EW 



A1 igual que para los marcos 10, de las envolventes de respuesta maximas del marco B del 
modelo IR6 (figura 10) se observa que su respuesta no lineal es mas importante que la observada 
en el modelo IR6B (figura 9). La demanda maxima de ductilidad de entrepiso es de 2.5 (k2=2%, 
nivel 5) asociada a rigideces secantes minimas inelasticas (pico a pico) del 40% de su rigidez 
elastica (k2=2%, nivel 5), a un valor minimo observado del promedio de las rigideces secantes 
inelasticas de medio ciclo del 72% de su rigidez elastica (k2=5% y 10%, nivel 5). Las distorsiones 
de entrepiso maximas sobrepasan a las calculadas en su diseno en los niveles 1 a 14, y al limite 
maximo permitido por el RDF-93 A=1.2% en los niveles 3 a 7. La distorsion maxima corresponde 
al cuarto entrepiso y es A=1.63% (k2=10%) calculada esta distorsion con respecto a la position 
inicial de equilibrio. Se presento respuesta elastica exclusivamente en los niveles 19 a 21. En 
aquellos niveles donde se presento respuesta inelastica, la distorsion de entrepiso asociada a su 
primera fluencia vario de Af e =0.33% (k2=2%, nivel l) a Af e =0.53% (k2=2%, nivel 14), obteniendose 
en algunos entrepisos respuestas elasticas para distorsiones de A=0.43% (k2=2%, nivel 19), 
coincidiendo esta variacion con la observada en IR6B para el marco B. El intervalo de variacion 
observado en los marcos B es menor al observado en los marcos 10. La distorsion de entrepiso 
maxima asociada al medio ciclo de histeresis de carga maxima fue Ame=2.4% en los niveles 4 y 5 
(k2=10%), valor que se encuentra muy por encima del limite maximo permitido por el RCDF-93 
(A=1.2%). 



La fluencia se presento principalmente en vigas, excepto de una columna exterior del 
nivel 1, que experimento fluencia incipiente en su base (de una a cinco veces en 2500 pasos de 
integration analizados). Se pueden hacer las mismas observaciones que se han hecho para el 
marco B de IR6B con respecto a las respuestas dinamicas maximas, dependiendo de las 
pendientes posteriores a la fluencia de los elementos estructurales de 2%, 5% o 10%. 
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El cortante basal que resiste el marco B de IR6 ante la accion del acelerograma S56-EW 
varla entre 6.6% (k2=2%) y 7.4% (k2=10%) del peso total del edificio IR6, por lo que se puede 
estimar crudamente que su resistencia a cortante basal en direccion Y debe variar entre el 26% y 
30% del peso de IR6, y que la sobrerresistencia observada para este edificio debe variar 
crudamente entre 2.08 y 2.40. Por tanto, se puede concluir, al igual que marco el marco 10, que la 
respuesta dinamica del marco B del modelo IR6B es mas favorable en todos los rubros que la del 
marco B del modelo IR6, por las mismas razones expuestas anteriormente. 



4. COMENTARIOS FINALES 



Se presento el estudio de estructuras esbeltas de 21 niveles de acero estructural, 
disenadas para distorsiones de entrepiso de diseno cercanas al limite maximo A=1.2% permitido 
por el RCDF-93, para marcos ductiles. El diseno de cada estructura irregular se realizo utilizando 
que el factor reductivo de fuerzas sismicas con fines de diseno (Q 1 ) establecido por las NTCS-95. 
Se estudiaron dos modelos, el modelo IR6B que cumple con todas las disposiciones de diseno de 
las NTCS-95 para estructuras irregulares y el modelo IR6, donde se estudio una variacion en la 
interpretacion de las disposiciones de las NTCS-95 de como deben revisarse las restricciones de 
desplazamientos laterales para fines de diseno (Ad). Se efectuaran analisis dinamicos no lineales 
empleando acelerogramas representatives para la zona del lago de la ciudad de Mexico, cuyas 
caracteristicas esten relacionadas al sismo maximo esperado y contemplado por el reglamento 
vigente (acelerogramas registrados durante el sismo del 19/09/85) y acelerogramas artificiales 
que pudieran estar contemplados en reglamentos futuros pues, de hecho, el sismo de M s =8.2 
proveniente de las costas de Guerrero es uno de los escenarios sismicos que se estudiaron en la 
revision de los espectros de diseno para la propuesta de NTCS-2004, que aun no entran en vigor. 
Los resultados obtenidos para los modelos en estudio permiten concluir, de manera preliminar, lo 
siguiente. 

Conforme a lo observado en el estudio del modelo IR6B, las disposiciones de diseno para 
estructuras irregulares de las NTCS-95 permiten llegar a disenos razonablemente seguros para 
estructuras esbeltas con base en marcos de acero estructural aun cuando estos se disenen para 
altas distorsiones (cercanas al limite maximo A=1.2% estipulado por el RCDF-93) y factores de 
comportamiento sismico altos (Q=4), ya que aun para acelerogramas asociados a sismos 
postulados intensos, la respuesta no lineal es moderada, estable y concentrada en las vigas. Las 
demandas maximas de ductilidad de entrepiso (p) son pequenas e incluso, bastante inferiores al 
factor de comportamiento sismico Q asociado a su diseno. Cabe senalar que aunque a y Q no son 
exactamente el mismo concepto, si existe relacion entre ambos parametros, y en este respecto, 
este estudio y otros que conduce el primer autor estan encaminados tambien para establecer 
relaciones entre ambos parametros. Sin embargo, como se observo para el marco B, existe la 
posibilidad que las distorsiones de entrepiso maximas sobrepasen a las calculadas en su diseno, 
aunque en el caso particular de IR6B no se observo que sobrepasaran al limite maximo A=1.2% 
permitido por el RCDF-93. 

Como se indico en secciones anteriores, las disposiciones de diseno por irregularidad de la 
propuesta de NTCS-2004 que competen a los edificios en estudio se mantienen iguales, por lo que 
tambien deben esperarse disenos seguros con la futura norma para este tipo de edificios. 

Conforme a lo observado en el estudio del modelo IR6, la variacion estudiada en la 
interpretacion de las disposiciones de las NTCS-95 de como deben revisarse las restricciones de 
desplazamientos laterales para fines de diseno (Ad) tambien permiten llegar a disenos 
razonablemente seguros para estructuras esbeltas con base en marcos de acero estructural en 
cuanto a niveles de esfuerzos y demandas no lineales, aun cuando estos se disenen para altas 
distorsiones (cercanas al limite maximo A=1.2% estipulado por el RCDF-93) y factores de 
comportamiento sismico altos (Q=4), ya que aun para acelerogramas asociados a sismos 
postulados intensos, la respuesta no lineal es moderada, estable y concentrada en las vigas. Sin 
embargo, se observa que para este caso, las distorsiones de entrepiso maximas sobrepasan 
notablemente, tanto a las calculadas en su diseno, como al limite maximo de distorsion (A=1.2%) 
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permitido por el RCDF-93. A estos niveles de distorsion se pueden esperar danos considerables en 
elementos no estructurales, asi como problemas de funcionamiento y operacion del edificio 
asociados a la conducta humana durante sismos, por lo que debera procurarse ser mas claro y 
explicito en las futuras normas de diseno por sismo y en sus comentarios como deben revisarse 
las restricciones de desplazamientos laterales para fines de diseno en el diseno sismico de 
estructuras irregulares. Las demandas maximas de ductilidad de entrepiso (li) son moderadas, 
ligeramente mayores a las que se presentaron en IR6B e inferiores al factor de comportamiento 
sismico Q asociado a su diseno. 

Con respecto a la pendiente posterior a la fluencia supuesta en las curvas primarias de los 
elementos estructurales (2%, 5% o 10%), no se presentaron grandes diferencias en las respuestas 
dinamicas maximas de los modelos IR6 e IR6B, observandose diferencias mas apreciables para 
los marcos en direccion Y que para los marcos en direccion X. En general, el aumento en el valor 
de la pendiente posterior a la fluencia tiende a disminuir las demandas de ductilidad, a 
incrementar las rigideces secantes inelasticas, la distorsion de entrepiso asociada a la primera 
fluencia y aumentar la capacidad a cortante de los entrepisos (sobrerresistencia). Habra que 
evaluar en estudios futures la influencia de la pendiente posterior a la fluencia en sistemas 
estructurales que sean sujetos a una respuesta no lineal mas amplia, con el objeto de determinar 
si a estos niveles de excitacion se observan o no diferencias mayores en las respuestas dinamicas 
maximas con respecto a las valoradas en IR6 e IR6B. 

Ademas, se deben realizar estudios adicionales donde se evaluen otras fuentes de rigidez 
y sobrerresistencia adicional, como la participacion del sistema de piso en la rigidez de las vigas 
considerandolas como vigas de seccion compuesta, en la respuesta dinamica no lineal y ver de que 
manera afectan la respuesta de estos disenos cuando se sujetan a los acelerogramas empleados en 
este estudio u otros asociados a los espectros de diseno. 

En el presente estudio se proponen algunos parametros que permiten valorar mas 
claramente a la respuesta no lineal de la estructura, ya que en ocasiones los lazos histereticos, el 
mapeo de fluencias, determinacion de distorsiones maximas de entrepiso, calculo de demandas de 
ductilidad y de rigideces secantes “pico a pico” no son suficientes para comparer la intensidad de 
respuesta no lineal de una estructura cuando se sujeta a movimientos del terreno de distintas 
caracteristicas, o se quiere evaluar el efecto de distintas fuentes de sobrerresistencia. En este 
respecto, el presente estudio ilustra que varios de los parametros dinamicos propuestos resultan 
muy utiles para valorar la importancia y caracteristicas de la respuesta no lineal mas claramente 
que a partir de los lazos de histeresis, y que miden de una manera indirecta la fatiga (de bajo y 
alto ciclaje), entre los que destacan, (l) el cociente entre el numero de ciclos donde se presenta 
respuesta inelastica (Cine) y el numero total de ciclos a que fue sujeta la estructura ante un 
acelerograma dado (Ctot) y, (2) el cociente entre el promedio de las rigideces de entrepiso secantes 
(kiprom) asociadas a los medios ciclos histereticos inelasticos y la rigidez elastica del entrepiso (k e i). 
De hecho, este ultimo parametro resulta mas idoneo para valorar la respuesta inelastica de una 
estructura sujeta a un acelerograma dado que, por ejemplo, el cociente entre la rigidez de 
entrepiso secante minima inelastica (kimin) o “pico a pico” asociada a los ciclos histereticos de 
amplitud maxima y la rigidez elastica del entrepiso (kei). Este ultimo parametro se utiliza mucho 
en investigaciones experimentales, donde se tiene un mejor control de las amplitudes de 
deformacion y de los ciclos a que se sujetan los modelos experimentales, situacion que no se da 
cuando una estructura es sujeta a movimientos aleatorios del terreno. 

El primer autor considera que para que los reglamentos evolucionen a criterios de diseno 
por desempeno y por desplazamiento de una manera realista y transparente, se deben conducir 
estudios analiticos encaminados a determinar los niveles de distorsion de entrepiso que deben 
permitirse a edificios, dependiendo del material y sistema estructural utilizado. Resulta 
incongruente tener los mismos limites para todos los sistemas y materiales, independientemente 
de si se trata de estructuras a base de mures de mamposteria o marcos ductiles de acero, como se 
establece en el RDF-93. Existen algunos indices que pueden proponerse a partir de los estudios 
experimentales que se han hecho hasta la fecha en modelos de reducido tamano o escala; de 
hecho, en este tipo de estudios se basan los limites de distorsion que ahora se proponen en el 
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apendice A de las NTCS-2004, donde ya se diferencian las distorsiones que deben permitirse a 
distintos sistemas estructurales. Sin embargo, muchos de los limites propuestos se basan en 
estudios de sistemas muy simples y, por lo tanto, para aumentar su confiabilidad deben ser 
complementados por modelos analiticos que permitan estudiar un mayor numero de casos y 
variables, utilizando modelos mas complejos que sean representatives de estructuras reales, no 
de modelos simples de laboratorio. 

Por ello, en el presente estudio se determino la distorsion de entrepiso asociada a la 
primera fluencia de cualquier elemento estructural que forma parte del entrepiso (Af e ), ya que 
este indice puede ser muy util para proponer con bases mas firmes, en un futuro, limites de 
distorsion de entrepiso de diseno en los reglamentos que sean congruentes con el material y 
sistema estructural que pretenda utilizarse. Los estudios aqui conducidos muestran una gran 
variabilidad en el valor de la distorsion de entrepiso asociada a la primera fluencia, la que oscilo 
entre Af e =0.28% y Af e =0.97%, obteniendose en algunos entrepisos respuestas elasticas para 
distorsiones tan altas como A=1.07%. A1 parecer, lo anterior puede deberse a un nivel de esfuerzos 
de diseno disparejo en los elementos estructurales en los distintos entrepiso asociado al hecho de 
emplear secciones tipo en el diseno, aspecto que no se ha estudiado con profundidad y que se 
analizara con mayor detalle en estudios futuros. El primer autor ha realizado este mismo tipo de 
analisis en estructuras esbeltas con base en marcos de concrete reforzado, donde se ha observado, 
tanto una variabilidad menor en el valor de Af e , como un valor absolute menor de Af e con respecto 
al presentado en este estudio para IR6 e IR6B. 

Finalmente, quiere recalcarse que los estudios se realizaron utilizando el programa de 
analisis no lineal DRAIN-2DX, que permite modelar exclusivamente las fallas por flexion y la 
accion bidimensional. No se pueden obtener mas que ideas aproximadas de la redistribucion de 
esfuerzos tridimensional (como se hizo en los modelos de ocho niveles que no se presentan, pero 
no en IR6 e IR6B), ignorando efectos significativos como es la torsion dinamica no lineal y la 
accion multidireccional de las aceleraciones del terreno. De igual manera, no se pueden evaluar 
potenciales fallas por cortante, flexocortante y torsion simultaneamente con las fallas por flexion, 
flexocompresion y flexotension. Se considera importante realizar analisis adicionales en las 
estructuras en estudio tomando en cuenta la accion de estos efectos, cuando los recursos de 
compute y de software, permitan realizar este tipo de estudios en modelos tan complejos (aunque 
aparentemente simples) como son los modelos IR6 e IR6B presentados. Para ellos, se requieren 
programas de analisis tridimensionales que modelen con rigurosidad los efectos descritos en el 
intervalo de comportamiento no lineal y en cada paso de integracion en el tiempo. Programas 
como el SAP2000 ofrecen un modelado no lineal de baja aproximacion para este tipo de 
estructuras, sobre todo si toma uno en cuenta que la respuesta paso a paso no se calcula con 
rigurosidad en cada paso de integracion, sino con un procedimiento “equivalente”, segun se 
documenta en el manual de dicho software. 
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RESUMEN 



Se presenta un modelo de generacion de acelerogramas artificiales consistente con la 
sismicidad de la region oeste de la Argentina. En su formulacion se hace uso del concepto de 
espectro evolutivo, donde la evolucion en el canrpo frecuencial esta materializado por una 
discretizacion del espectro de Amplitudes de Fourier, definiendo un proceso uniformemente 
modulado en cada franja. 

La modulacion en el carnpo temporal se hace mediante una funcion tipo gamma cuyos 
parametros se determinan en funcion de los momentos temporales del registro, los que a su vez 
son estimados en base a una regresion lineal en funcion de la magnitud y la distancia epicentral; 
de esta manera se concreta la relacion con la sismicidad local. 

ABSTRACT 



A model for the artificial generation of accelerograms consistent with the local seismicity 
of the western region of Argentina is presented. It is based upon evolutionary spectrum concepts 
where the evolution in frequency content is taken into account by means of the discretization of 
the Fourier Amplitude Spectra and the definition of a uniform modulated process in each strip. 

The amplitude modulating function is considered to be of gamma type and its parameters 
are determined in function of the temporal moments of the record, which in turn are estimated by 
a multilinear regression analysis in term of magnitude and epicentral distance. In this way the 
relation to local seismicity is considered. 



1. INTRODUCTION 



La Republica Argentina cuenta con una amplia zona adyacente a la Cordillera de los 
Andes con una actividad sismica importante, evidenciada en las cronicas historicas y en los 
ultimos anos, en el catalogo instrumental. Se cuenta con algunos registros obtenidos 
principalmente en la zona cuyana, que es la de mayor actividad sismica, pero ellos son todavia 
insuficientes para las aplicaciones que el calculo estructural exige actualmente. Esta amplia area 
de la region oeste-central ha sido microzonificada, por lo que el conocimiento que en ella se tiene 
de las condiciones del suelo, tectonica y geomorfologia, es mas acabado que en otras regiones 
argentinas. Se estan realizando tambien grandes esfuerzos para encarar estudios de la 
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sismicidad en la Provincia de Salta, en el norte del pais, que es la segunda zona de importancia 
considerando la peligrosidad sismica, donde no se cuenta con ningun registro de sismos intensos, 
a pesar de que hay evidencias historicas de su ocurrencia. 

Por todo ello es que se considero factible encarar un trabajo de investigacion con el 
principal objetivo de definir un modelo de generation de sismos artificiales probables de ocurrir en 
principio, en la region cuyana argentina, para ser aplicados al calculo estructural. Este trabajo 
dio como resultado la realization de dos modelos : uno que reproduce un registro objetivo, con el 
que se realizaron numerosas experiencias numericas a fin de obtener historias de aceleraciones 
concordantes con las de la base de datos disponible, constituida en principio por registros del 
sismo que acontecio en Enero de 1985 en la provincia de Mendoza, cuya sismicidad frecuente se 
debe a movimientos de fallas superficiales cercanas. Basado en este modelo, se propone un 
modelo de generation de registros artificiales, que hace uso de la tecnica de regresion, 
considerando como datos basicos la magnitud y la distancia epicentral que a veces son los unicos 
datos con que puede contarse en algunos sitios, y es el que se explica a continuation. 



2. METODOLOGIA 



Las distintas propuestas y estrategias para encarar la solution de esta problematica, 
pueden englobarse en uno de los dos grupos siguientes : modelos empiricos y modelos sismologicos. 
Normalmente, los modelos de generation con aplicacion al calculo estructural, que tienen un 
sentido eminentemente ingenieril, siguen la linea de los metodos empiricos. Para ello es necesario 
contar con un buen numero de registros de buena calidad de varios terremotos originados en las 
distintas fuentes sismicas que afecten el sitio en estudio. Se ha sorteado esta dificultad en varios 
paises con pocos registros, usando modelos definidos para otras regiones, y modificando sus 
parametros a condiciones locales, normalmente haciendo regresiones en funcion de magnitud y 
distancia epicentral, que son las magnitudes macrosismicas mas accesibles en todas partes. Los 
modelos sismologicos tienen una aplicacion mas relacionada al ambito de la geofisica y geologia; 
son en general mas complejos y la determination de muchos parametros para definirlos hace 
necesario contar con information y conocimientos sismologicos, que la mayoria de las zonas 
sismicas del tercer mundo no cuentan aun. 

Para este trabajo de investigacion se conto con pocos registros, correspondientes a dos 
eventos sismicos acaecidos en las Provincias de Mendoza y San Juan. Se tiene tambien 
information complementaria valiosa, resultante de los estudios de microzonificacion llevados a 
cabo en las areas mas importantes de las mismas. Se penso conjugar en la definition del modelo 
ambas caracteristicas (empiricos y sismologicos); de esta manera se aprovecha la principal 
ventaja de los modelos sismologicos, pero sin perder de vista la aplicacion ingenieril del mismo. 



3. BASES CONCEPTUALES 



El modelo de regresion se basa en el modelo de simulation desarrollado por la autora 
(Orosco, 2003) que hace uso de un acelerograma objetivo a reproducir, por lo que sus bases 
conceptuales se identifican con el mismo. 

A fin de considerar la no estacionariedad de los registros sismicos tanto en el campo del 
tiempo como en el de la frecuencia, el modelo se basa en considerar el movimiento del suelo como 
constituido de varios componentes uniformemente modulados. Esta metodologia fue utilizada 
entre otros, por Der Kiureghian y Crempien (1989) y Carli y Faravelli (1990). Tambien en esta 
linea se inscribe la propuesta de Saragoni y Hart (1974) con la diferencia que en el modelo de 
Saragoni se discretiza el movimiento en el dominio del tiempo mientras que en el modelo de 
Crempien la funcion de modulation varia en intervalos determinados en el campo de la 
frecuencia, por lo que presenta continuidad temporal. 
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La aceleracion del suelo debido a un movimiento sismico a(t) , se representa mediante la 
sumatoria de k procesos uniformemente modulados x„„ it) , correspondientes a las k bandas 

ne k 

consideradas en el campo de las frecuencias. Luego 

a (0 = Tj*ne t (0 

k 



donde el proceso uniformemente modulado se define como 

= (p k (t)x e (t) 

por lo que la aceleracion del suelo 

ait) = Yj<PkK 

k 

X t , t it) define en cada banda de frecuencia, un proceso estacionario estadisticamente 

independiente, Gaussiano, de media nula y varianza unitaria, que se expresa mediante la 
conocida expresiom 



x e it) = ^ 2-JSico)A(o cos icDjt + (j) i ) 

1=1 

En la expresion anterior Sico) es la funcion de densidad espectral de potencia del proceso 
estacionario en cada franja. Con el objetivo de representar mejor los acelerogramas de la zona 
cuyana argentina, se experimento numericamente con varias funciones de densidad espectral de 
potencia (FDEP), partiendo de una simple funcion constante, utilizada en el modelo de Crempien 
para simular sismos chilenos de subduccion y los del oeste norteamericano (Der Kiureghian y 
Crempien, 1989; Crempien y Araneda, 1993). Cumpliendo con la exigencia de varianza unitaria, 
la FDEP constante se define como : 



S k (co) 



1 



2 K +I -oj k ) 



siendo nula fuera de la franja k considerada. 

La funcion espectral adoptada en este trabajo, incorpora en su definicion al espectro de 
Boore de la siguiente manera (Orosco, 2003; Queck et al., 1990): 

= ^ 7 A» 

2m 

donde I , denominada intensidad, se define como (Queck et al., 1990) 



/ = \\cpit)fdt 

o 

A a ico) es el espectro de Boore que se expresa como (Boore, 1990) 



A a ico) = C S^co) S 2 ico) S 3 ico) 
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donde C es un factor de escala, S x ( CO ) es el espectro de fuente, S 2 (co) tiene en cuenta las 
posibles amplificaciones por efectos locales de sitio y S 3 ( CO ) es el factor de atenuacion. 



En las expresiones (3) y (7) , <p k (t ) son funciones de modulacion en funcion del tiempo, 

cuyos parametros se determinan considerando los momentos espectrales del registro, los que se 
estiman para este modelo mediante una regresion multilineal, en funcion de la magnitud y la 
distancia epicentral. Se explica a continuacion, detalles del proceso de regresion realizado y la 
base de datos utilizada; en parrafos posteriores se retoma la discusion acerca de la 
determinacion de la funcion de modulacion. 



4. BASE DE DATOS 



Las tablas 1 y 2 resumen caracteristicas de los sismos utilizados en este estudio. Se los 
dividio en dos grupos, considerando la magnitud. 



Tabla 1 Registros del grupo U1 



Sismo 


M-H 


Registro 


A 


A 

max 


8 


h 






(km) 


(cm/seg 2 ) 




(cm 2 /seg 3 ) 


Granada 


M=4.9 


Santa Fe (T) 


40.4 


29.7 


0.67 


0.57 


1984 


H=12 km 


Santa Fe (L) 


40.4 


35.0 


0.67 


0.58 


Beznar (L) 


20.0 


19.80 


0.78 


0.13 


Dursunbey 

1979 


M=5.9 
H=7 km 


Dursunbey 


7.8 


175.2 


0.76 


14.02 






AYLO 


27.2 


98.1 


0.75 


4.17 






AMLO 


28.2 


274.7 


0.57 


36.52 






AMTR 


28.2 


333.5 


0.73 


47.42 






LALO 


50.2 


54.0 


0.61 


12.36 






LATR 


50.2 


59.8 


0.63 


14.28 






LHLO 


33.1 


192.3 


0.76 


17.16 






LHTR 


33.1 


412.0 


0.66 


91.46 






MALO 


20.0 


58.9 


0.65 


1.36 






MATR 


20.0 


186.4 


0.76 


10.03 






POLO 


42.1 


98.1 


0.34 


5.29 






POTR 


42.1 


107.9 


0.34 


5.29 


Mendoza 


M=5.7 


SMLO 


34.3 


78.5 


0.44 


5.78 


1985 


H=12 km 


SMTR 


34.3 


78.5 


0.49 


4.57 






TULO 


55.9 


107.9 


0.43 


9.16 






TUTR 


55.9 


107.9 


0.37 


12.47 






USLO 


77.6 


58.9 


0.40 


3.38 






USTR 


77.6 


107.9 


0.37 


7.13 






ARC 


8.6 


286.3 


0.44 


85.84 






DWY 


16.6 


334.0 


0.62 


94.36 






LWD 


24.0 


291.4 


0.66 


64.31 


Whittier 

Narrows 

1987 




LFS 


31.0 


131.5 


0.52 


13.36 


M=5.9 


PMC 


37.9 


122.1 


0.48 


16.02 


H=14.7 km 


CCY 


51.8 


113.1 


0.61 


11.03 




NWS 


62.0 


85.2 


0.44 


9.10 


Coalinga 














9/5 

Coalinga 

9/7 


M=5.3 


BNT 


12.4 


89.7 


0.74 


3.07 


M=5 


BNT 


12.4 


137.9 


0.58 


12.58 


M=6 


BNT 


12.4 


330.5 


0.57 


90.12 


Coalinga 


22/7 
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Un primer conjunto (llamado Ul) engloba acelerogramas de eventos de magnitud 
comprendida en el intervalo [4.9, 6]; en la tabla 1 se resumen las principales caracteristicas de 
estos eventos. Se detalla la magnitud (M ), profundidad ( H ) en km, distancia epicentral (A) en 
kilometros, valor pico de la aceleracion (A ma x) en cm/seg 2 , factor de Vanmarcke (8) e Intensidad 
de Arias (Ia) en cm 2 /seg 3 . Un segundo conjunto de datos (al que se denomino U2), ademas de los 
sismos moderados del grupo Ul, incluye registros correspondientes a sismos intensos de mayor 
magnitud, cuyas caracteristicas se muestran en tabla 2. 



5 FORMULAS DE REGRESION 



Evaluando las distintas alternativas y los distintos parametros que pueden ser objeto de 
una regresion, se juzgo conveniente optar por relacionar la energia de los sismos con la distancia 
epicentral y la magnitud. Se cuenta con expresiones de atenuacion de la aceleracion maxima de 
terremotos de corteza superficial para la zona en estudio. Se considero la hipotesis de que la 
energia responde a una ley de comportamiento similar a la de la aceleracion maxima, lo que se 
confirmo en el caso de la base de datos usada. 



Tabla 2 Registros de magnitud mayor a 6 



Sismo 


M-H 


Registro 


A 

(km) 


A 

max 

(cm/seg 2 ) 


8 


h 

(cm 2 /seg 3 ) 


Caucete 1977 


M=7.4 
H= 25 km 


SJLO 

SJTR 


70.0 

70.0 


186.4 

189.3 


0.51 

0.52 


107.35 

133.28 


Loma Prieta 
1989 


M=7 

H=17.6 

km 


Oackland (L) 
Piedmont (T) 
Saratoga (L) 
Saratoga (T) 
Santa Cruz (L) 
Santa Cruz (T) 


95.2 

93.3 
29.8 

29.8 

19.9 
19.9 


306.0 

110.0 
316.2 

494.5 

101.5 
433.1 


0.71 

0.68 

0.62 

0.65 

0.46 

0.45 


81.75 

5.12 

109.40 

146.23 

204.74 

263.73 


Northridge 

1994 


M=6.7 

H=18.4 

km 


Rinaldsi Stn 
(L) 

Rinaldi Stn (T) 

SylmarC.Stn. 

(L) 

SylmarC.Stn 

(T) 

Sylmar Stn E 
(L) 

Sylmar Stn E 
(T) 

Sylmar Hosp 
(L) 

Sylmar Hosp 
(T) 


10.0 

10.0 

12.0 

12.0 

13.0 

13.0 

16.0 
16.0 


804.4 

559.2 
618.0 

696.5 
735.8 
461.1 

863.3 

598.4 


0.81 

0.81 

0.73 

0.74 

0.70 

0.75 

0.73 

0.74 


796.89 

435.87 

425.22 

635.02 

409.93 

319.21 

177.14 

245.44 



5.1 Formula de regresion basada en expresiones de atenuacion existentes 



Se han propuesto varias relaciones de atenuacion para el oeste argentino, que describen la 
variacion de la aceleracion pico con la distancia. Bufaliza desarrollo una relacion de atenuacion 
usando datos registrados en la zona cuyana argentina, que se expresa como (INPRES, 1995) 



log(a max ) = -0.7837 + 0.353M - 1 .598 1 log R - 0.00099 R 
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donde M es la magnitud correspondiente a las ondas superficiales para M > 6 y a las ondas de 
cuerpo para M < 6; R es la distancia hipocentral en km y a ma x es la mayor de las aceleraciones 
pico de las dos componentes horizontales, expresada en fracciones de la aceleracion de la 
gravedad. En el caso de sitios cercanos al piano de falla, la ecuacion (9) deberia ser modificada 
para considerar la saturacion del movimiento del terreno, de la siguiente manera 

log(a max ) = -0.7837 + 0.353M - 1.59811og D - 0.00099 D 
donde D se expresa corno 



D = ^]r 2 + (4.32 10 0 073M ) 2 

Se realizo la regresion multilineal mediante la ecuacion 

log(/l ; . ) = a + bM + c log A + cl A 

donde a, b, c y d son los parametros a determinar. Se observa en la ecuacion (12) que son los 
momentos temporales del registro los que se determinan mediante la regresion ya que describen 
de algun rnodo la energia que contiene y el modo en que esta se distribuye, con lo que la 
informacion de la genesis, paso y sitio esta contemplada de modo sencillo y general, pero 
suficiente a los fines del calculo estructural. 



Figura 1 Relacion entre energia y distancia epicentral para sismos grupo U1 



5.1.1 Analisis de los resultados de la regresion 



En las tablas 3, 4 se resumen los resultados de la regresion para el conjunto de datos 
denominado Ul. En ellas Std son las desviaciones estandar y t es el valor de la prueba de 
Student. Se realizaron ademas las pruebas de significancia de Fischer (Herva et al., 1981) a fin 
de evaluar si al menos uno de los componentes de la formula de regresion utilizada (excepto la 
constante) explican el valor de los parametros que se desean estimar. 



Tabla 3 Resultados de la regresion 1 para grupo Ul 



Parametro 


a 


b 


c 


d 




Std 


Std 


Std 


Std 




t 


T 


t 


t 


An 


-9.59 


3.92 


-1.29 


0.013 


u 


3.99 


0.60 


0.94 


0.031 




-2.41 


6.53 


-1.37 


0.42 


A, 


-16.52 


5.05 


-0.53 


0.004 


1 


4.47 


0.67 


1.05 


0.035 




-3.70 


7.50 


-0.50 


0.10 


A 2 


-23.80 


6.24 


0.34 


-0.007 




6.03 


0.91 


4.42 


0.047 




-3.94 


6.86 


0.24 


-0.14 


t f 


-6.92 


1.38 


0.85 


-0.017 


J 


2.35 


0.35 


0.55 


0.018 




-2.94 


3.89 


1.53 


-0.93 


t p 


-4.28 


0.63 


0.99 


-0.012 


e 


2.37 


0.36 


0.56 


0.018 
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-1.80 


1.77 


1.76 


-0.66 




Tabla 4 Resultados de la regresion 1 para el grupo U2 




Parametro 


a 


b 


c 


d 




Std 


Std 


Std 


Std 




t 


T 


t 


t 


K 


0.23 


2.41 


-1.66 


0.010 




2.56 


0.26 


0.74 


0.022 




0.09 


9.17 


-2.24 


0.47 


JL 


-5.39 


3.10 


-0.42 


0.014 


1 


2.78 


0.29 


0.80 


0.024 




-1.94 


10.84 


-0.52 


-0.57 


X 2 


-10.63 


3.73 


0.99 


-0.040 




1.21 


0.37 


1.04 


0.031 




-2.84 


10.14 


0.88 


-1.28 


h 


-3.43 


0,61 


1.16 


-0.025 


1.24 


0,12 


0.35 


0.011 




-2.84 


4.90 


3.31 


-2.34 


t„ 


-4.41 


0.51 


1.36 


-0.027 


e 


1.24 


0.13 


0.36 


0.011 




-3.54 


3.02 


3.79 


-2.52 



Analizando los valores para el grupo U1 de datos, se tiene que el momento /l 0 tiene un 

valor F de 16.07 que es mucho mayor que F(0.0l; 3, 27)=4.60, donde 0.01 es la densidad de 
probabilidad, 3 es el numero de parametros a determinar (sin considerar la constante) y 27 son 
los grados de libertad (cantidad de datos de la muestra menos numero total de parametros a 
estimar). Se procedio del mismo modo con el resto de los parametros, obteniendose el valor mas 
bajo de F para la duracion efectiva (6.12). 

Acerca de la prueba de Student una evaluacion general de los parametros a,b,c y d 
considerando t (0.10; 27) =1.70 lleva a la conclusion que el termino con menor significancia en 
todos los casos es d , ya que para el mismo, el valor mas alto obtenido fue de 1 1 1 =0.93. En cuanto 
a la constante c , los valores correspondientes a y X 2 son malos; para los restantes 

parametros, no obstante sus bajos valores, se pueden aceptar. En general, la constante b explica 
bien el modelo propuesto en todos los casos. De este analisis puede inferirse que la distancia 
epicentral no explica muy bien el modelo para esta base de datos. 

Considerando ahora el grupo de datos U2 los valores de F para los momentos espectrales 
se ubican entre 35 y 45 mientras que para la duracion total y efectiva, F es aproximadamente 12; 
considerando el nivel F(0.0l; 3, 43) = 4.31, se concluye que en el modelo al menos alguna de sus 
variables es significativa para explicar los valores encontrados. 

En cuanto al valor de las constantes a,b,c y d, se observa que sus valores t 
aumentaron y teniendo en cuenta un mejor nivel de estimacion (t(0.02; 43)= 2.42) permite inferir 
que el termino que multiplica a d explica la duracion total y efectiva pero no asi los momentos 
temporales. En cuanto a la constante c , tal como en el caso anterior no explica los momentos de 
orden 1 y 2 pero si explica, desde el punto de vista estadistico, el valor de la energia, la duracion 
total y la efectiva. El valor de b es muy bueno en cuanto a su significancia en el modelo de 
regresion propuesto. Como conclusion, la magnitud y el logaritmo de la distancia epicentral 
explican muy bien el valor de la energia. Para la duracion efectiva y total todos los parametros 
son significativos, mientras que para los momentos y X 2 , solo la magnitud tiene relevancia. 
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5.2 Modelo de regresion propuesto 



En vista del analisis estadistico realizado a los resultados de la regresion hecha en base a 
las formulas de atenuacion existentes, se propone una expresion que determina la energia del 
proceso y la duracion efectiva en funcion de dos parametros mas la constante, dejando de lado la 
distancia epicentral (termino que multiplica la constante d ) por los bajos valores de significancia 
obtenidos en el caso anterior. La nueva formula de regresion tiene la expresiom 

log(/l 0 ) = a + b log (M ) + c log(A) 
log0 0 ) = a + b log (AT ) + c log(A) 

Con el modelo anterior, en la estimacion de los valores de los momentos temporales de 
orden 1 y 2 solo la magnitud tiene relevancia; por ello la formula de regresion para los momentos 
de orden superior se expresa en funcion de A 0 y t e de la siguiente manera : 

log(l ( ) = a + b log(l 0 ) + c log (t e ) 1 = 1,2 

donde A () y t son determinados mediante (13). 

Los resultados de la regresion efectuada para los dos grupos de datos se resumen en la 

tabla 4. 



Tabla 4 Resultados de la regresion con el modelo alternativo 



Parametros 


Grupo Ul 


Grupo U2 




a 


b 


c 


a 


b 


c 




Std 


Std 


Std 


Std 


Std 


Std 




t 


t 


t 


t 


t 


t 


/L 


-25.02 


21.25 


-0.91 


-14.44 


15.68 


-1.20 


U 


5.49 


3.19 


0.29 


3.03 


1.54 


0.25 




-4.56 


6.65 


-3.13 


-4.76 


10.18 


-4.73 


K 


-1.370 


1.125 


0.897 


-0.982 


1.063 


0.953 




0.174 


0.017 


0.037 


0.158 


0.013 


0.040 




-7.89 


65.90 


24.30 


-6.20 


79.46 


24.13 


A, 


-2.853 


1.259 


1.942 


-1.788 


1.118 


1.975 


z 


0.335 


0.033 


0.071 


0.298 


0.025 


0.074 




-8.51 


38.18 


27.24 


-6.00 


44.40 


26.56 




-5.81 


3.40 


0.64 


-6.42 


3.63 


0.69 




3.28 


1.91 


0.17 


1.61 


0.82 


0.13 




-1.77 


1.78 


3.65 


-3.98 


4.43 


5.11 



5.2.1 Analisis de los resultados del modelo alternativo de regresion 



Para el conjunto de datos Ul, los valores de F para el momento de orden 0 y la duracion 
efectiva son de 24.81 y 9.18 respectivamente, lo que considerando un nivel F(0.0l; 2, 28) = 5.45 se 
concluye que las variables del modelo tienen una alta significancia estadistica para explicar los 
parametros mencionados. En cuanto a A t y A 2 sus valores son 2,387 y 1,047 respectivamente, lo 
que estadisticamente significa un valor muy bueno de relacion entre los parametros involucrados. 
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Para el grupo de datos U2 los valores de F se incrementaron con respecto al caso anterior; 
se tiene para los momentos temporales las siguientes cantidades: 77, 3,489 y 1,363 
respectivamente, mientras que para la duracion efectiva es 19 y considerando un nivel F(0.0l; 2, 
42) = 5.18, se concluye que hay una relevante relacion entre los parametros estimados y las 
constantes obtenidas por la regresion. 

Analizando a,b y c con la prueba de Student para los datos del grupo Ul, los valores t 
de todas las constantes sobrepasan el nivel t(0.00l; 28)=3.66 excepto el caso de c (que solo esta 
a 0.005 por debajo de este valor) y las constantes a y b de la duracion efectiva (que satisfacen el 
nivel t(0.10; 28)=1.70). En cuanto al grupo U2, todas las variables sobrepasan el nivel 
t(0.00l;42)=3.55, lo que implica que el modelo es muy bueno desde el punto de vista estadistico. 
Asi, el logaritmo de la energia y de la duracion efectiva explican muy bien los momentos 
temporales de orden 1 y 2, lo que puede observarse del analisis de sismos reales. 

Como las regresiones se realizan sobre los logaritmos de la magnitud y distancia 
epicentral, pequenas variaciones en ellos provocan grandes diferencias en los resultados, por ello 
la regresion debe ser realizada con cuidado para obtener resultados confiables. 

Este modelo se denomina de “regresion 2” para diferenciarlo del anterior realizado en 
base a la formula (12) al que se identifica como modelo de “regresion 1”. Se hicieron simulaciones 
con ambas expresiones y en general se observa una mayor estabilidad en las formas de las 
funciones de amplitud determinadas en base a los momentos temporales obtenidos por el modelo 
de regresion 2 con respecto a las obtenidas por la regresion 1, en especial para registros de banda 
angosta que en general corresponde a campo cercano, mientras que las diferencias son 
insignificantes para el caso de eventos de banda ancha. 



6. PARAMETROS DE CONTROL DEL MODELO 



Los aspectos importantes a considerar para determinar los parametros que controlan el 
proceso de generacion de historias con determinadas caracteristicas se explicitan a continuacion. 
Se incluyen algunos ejemplos de aplicacion a fin de mostrar sus efectos e influencias en los 
registros obtenidos. 

a) Tipo de registro que se desea obtener: banda angosta o banda ancha. 

b) Distribucion de la energia en cada franja frecuencial. 

c) Diseno de la funcion de amplitud. 

d) Diseno de la funcion de densidad espectral de potencia. 

e) Control de la energia en cada franja. 

f) Duracion efectiva y total, y aceleracion maxima. 



6.1 Tipo de registro 



Es posible decidir si se desea obtener un acelerograma de banda angosta o de banda 
ancha. Tales caracteristicas son determinadas por el valor del factor de (Vanmarcke, 1980) del 
registro que se define como : 
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v KK 

donde los X i son los momentos espectrales del registro Ci(t) definidos por 

A.=pA A f (a>) i=0, 1, 2 

y A f (a>) es el espectro de amplitudes de Fourier. Para valores de 8 mayores a 0.66 se tienen 

registros de banda angosta y para valores menores a 0.66, registros de banda ancha. La base de 
datos de Mendoza arrojo como conclusion que para registros de banda angosta se obtenian 
mejores resultados en el proceso de simulation considerando una division de tres franjas en el 
campo frecuencial, mientras que para los de banda ancha la division en cinco franjas era la 
apropiada. Si bien esto puede limitar el universo a sismos similares solo a los del evento con que 
se cuenta mayor cantidad de datos, sigue siendo valido el concepto en el que se basa la division 
del espectro : considerar “paquetes de energia”, esto es dividir el espectro de energia del registro 
en tantas franjas como queden determinadas por los valores nulos o casi nulos del espectro de 
energia (considerese que esto es importante en el caso de utilizar el rnodelo que reproduce las 
caracteristicas de un sismo muestra). 

Es entonces posible “crear” un registro con un numero de paquetes de energia distinto a 
tres y cinco, si asi se desea. Se observo en el proceso de simular algunos sismos distintos a los 
cuyanos que para sismos muy impulsivos la division en dos bandas podria funcionar bien. 
Mientras que en sismos de banda muy ancha, con bajos valores del factor de Vanmarcke (menores 
de 0.45 por ejemplo), a veces es preciso recurrir a mas de 5 divisiones en el campo de la 
frecuencia. El sismo de Caucete (San Juan) es un ejemplo. 

En base a los datos con que se conto para este estudio, para distancia epicentrales 
mayores a 40 km se recomienda considerar historias de banda ancha, para distancias epicentrales 
menores de 30 km, se aconseja considerar banda estrecha, y para los situados entre 30 km y 40 
km se adopta uno u otro criterio de acuerdo al que mejor se adapte a los resultados que se desean 
obtener. Esto es solo una sugerencia, ya que se observo en la base de datos algunas variaciones a 
la tendencia general, pues las caracteristicas de los movimientos en el sitio de estudio dependen 
grandemente de las condiciones locales de suelo. 



6.2 Distribution de la energia en cada banda de frecuencia 



De todos los analisis numericos realizados se calculo el promedio del monto de la energia 
contenido en cada banda frecuencial, tanto para registros de banda angosta como de banda ancha. 
Esto no impide decidir “a priori” los porcentajes de la energia total con que se dotara a cada banda 
de frecuencia considerada: asi puede tenerse el caso de que se requiera probar el comportamiento 
dinamico no lineal de una estructura sujeta a una historia de aceleraciones en un rango de 
frecuencia determinado (cercano al periodo propio, por ejemplo), y al no contar con registros reales 
con suficiente energia en tal rango, es preciso simular historias en el tienrpo que respondan a esa 
caracteristica deseada. 

Para sismos de banda estrecha (para los que se requiere dividir el espectro en tres 
franjas), los valores medios de la fraction de energia total obtenidos para cada banda son: 0.77, 
0.20 y 0.05 con desviaciones estandar igual a 0.13, 0.13 y 0.04 respectivamente. Para sismos de 
banda ancha (con division en cinco franjas), se tiene que el promedio de las fracciones de la 
energia total en cada franja son 0.02, 0.24, 0.52, 0.20 y 0.03 respectivamente, siendo sus 
desviaciones estandar 0.02, 0.18, 0.28, 0.29 y 0.03. Se hace notar nuevamente que estos valores 
son solo indicativos. 
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6.3 Diseno de la funcion de amplitud 



A los fines de modular las amplitudes de los acelerogramas simulados, varios 
investigadores hicieron distintas propuestas de funciones de modulacion a lo largo del tiempo 
(Barbat et al, 1994). La mas simple es considerar un pulso rectangular, de duracion / ; . Esta 

funcion no describe apropiadamente las variaciones en el valor cuadratico medio de los 
acelerogramas, que es la caracteristica de la no estacionariedad en el tiempo. 

Se probaron varias funciones, a fin de reproducir del mejor modo posible los registros 
muestras. Los registros disponibles en la base de datos tendian a una funcion de amplitud tipo 
gamma, por lo que fue esta la adoptada para el modelo. 

Se define la funcion gamma de la siguiente manera 

m{t ) = at p e'*” 

en la que la variable t m =t — t 0 es la variable tiempo modificada por un parametro t 0 que tiene en 

cuenta el instante en que irrumpe el pulso, importante a tener en cuanta en el caso de sismos 
impulsivos, corno lo eran la mayoria de la base de datos mendocina. En lo que sigue, se obviara 
escribirlo de esta manera, pero siempre debe tenerse en cuenta este hecho. 

Figura 2 Funcion de amplitud para registros de banda ancha (Casos A, B, C y D) 



A fin de estimar los parametros a , J3 y y de las funciones de modulacion (pit) en cada 
franja, se impone la condicion de ergodicidad de los procesos no estacionarios (Bogdanoff y Kozin, 
1961) 

*f b 

A n = JV‘ (p 2 {t) dt = jV a 2 (t) dt 

o o 

con lo que se obtienen las siguientes expresiones algebraicas, en las que los parametros a , J3 y 
y estan en funcion de los momentos temporales obtenidos por regresion, 



a - 

P = 

7 = 



K 7 P+l 

TiP + l) 

% 

A 0 /t 2 ~ ^ 

A 0 A 2 -A] 



-1 



Los resultados de aplicar esta metodologia para un evento de banda ancha, de 32 km de 
distancia epicentral y una magnitud M = 5.4 grados se muestran como Caso A en la figura 2. 
Figura 3 Registros artificiales de banda ancha (casos A, B, C y D) 



Otra manera de estimar la forma de la funcion de amplitud es mediante relaciones entre 
los momentos de orden 1 y 2 (que justamente describen el modo en que la energia se distribuye) 
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con respecto al momento de orden 0 (que describe la energia contenida en rango frecuencial 
considerado) y cuyo valor es el determinado por la tecnica de regresion. 

Se definen entonces los parametros k l y k 2 corno 



K 

k 2 



K 

Aq 

a 2 



En la Tabla 5 se muestran los valores medios de los parametros k x , k 2 y t 0 para sismos 

cuyanos, obtenidos de los resultados del ajuste no lineal efectuado en cada franja con el modelo de 
simulacion que usa el sismo muestra como semilla de generacion. Como se aprecia el rango de 
valores en que tales parametros se ubican es amplio, como lo denotan los altos valores de sus 
desviaciones estandar (que es el valor que aparece entre parentesis). 

Los valores de k x y k 2 son fijados a voluntad, segun la apariencia de la funcion de 
modulacion que se desea lograr en cada franja frecuencial. Por ello, los valores medios se 
muestran a titulo de ejemplo. A fin de asegurar que la modulacion del proceso uniforme en cada 
franja sea mediante la funcion gamma, los valores de k x , k 2 y t Q deben satisfacer los siguientes 
requisitos, a fin de que los parametros de la funcion no resulten negativos, en especial el valor de 

P ■ 




Tabla 5 Valores medios de k x , k 2 y t 0 (sismos mendocinos) 



Lranja 


Banda estrecha 


Banda ancha 


K 


k 2 




k i 


k 2 




1 


2.51 

(0.59) 


9.04 

(3.24) 


1.61 

(0.78) 


6.84 

(3.02) 


71.82 

(59.61) 


2.49 

(2.24) 


2 


2.83 

(0.66) 


11.60 

(3.46) 


1.62 

(1.01) 


12.71 

(8.38) 


356.53 

(360.93) 


5.49 

(3.84) 


3 


2.88 

(0.59) 


11.53 

(3.78) 


1.87 

(1.01) 


10.04 

(5.25) 


195.09 

(183.03) 


4.40 

(2.70) 


4 








10.14 

(3.18) 


185.87 

(115.96) 


5.03 

(2.46) 


5 








11.84 

(4.26) 


287.58 

(224.81) 


5.04 

(2.50) 



La Tabla 6 muestra algunos ejemplos mas de pares de valores k x y k 2 que se agrupan 

en los denominados casos B, C y D y cuyas las funciones resultantes se aprecian en figura 2; las 
historias correspondientes y la evolucion temporal de la energia en cada ejemplo, se muestran en 
las figuras 3 y 4 respectivamente. 
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Tabla 6 Ejemplos de valores de k l y k 2 para sismos de banda ancha 



Franja 


Caso B 


CasoC 


Caso D 




K 


k 2 


K 


k 2 


K 


k 2 


1 


5.00 


43.17 


6.00 


59.51 


10.33 


140.24 


2 


4.00 


23.07 


15.00 


356.93 


7.95 


77.72 


3 


5.00 


36.52 


12.00 


195.09 


10.12 


135.00 


4 


5.00 


32.82 


11.00 


185.87 


8.79 


99.47 


5 


6.00 


52.70 


15.00 


287.58 


10.96 


168.38 



En todos los casos mostrados en la figura 2, se considero un valor nulo para el parametro 
t 0 . A fin de graficar su influencia en la forma de la funcion de amplitud buscada, se muestra en la 
figura 5 un detalle de las funciones de modulacion de la segunda y tercera franja correspondientes 
al caso D, en las que se vario el valor de t 0 . El intervalo de valores de este parametro adecuados 

para determinar una funcion gamma es reducidoi para el caso denominado Dl, t 0 e s de 1 segundo 

en ambas franjas; en el caso D2, los valores son de 2.5 y 2.80 segundos respectivamente y el caso 
D3 corresponde a valores de -1 segundo en cada una de ellas, manteniendo todas las demas 
condiciones constantes con respecto al caso D. La figura 6 muestra las historias resultantes. 

Por otro lado, en la figura 7 se observa el efecto de variar k l en la segunda y tercer franja 
frecuenciales. El caso denominado D4 corresponde a valores de k 1 cercanos a los limites inferiores 
del intervalo definido por (21) mientras que el caso D5 resulta de considerar valores muy 
cercanos al limite superior. La apariencia de los acelerogramas simulados se muestran en la 
figura 8. 



Figura 4 Intensidad de Arias para casos A, B, C y D. 



Figura 5 Efecto de variar t Q en la funcion envolvente de la segunda y tercera franja 
Figura 6 Registros artificiales de banda ancha (casos Dl, D2 y D3) 

Figura 7 Efecto de variar k , en la funcion envolvente de la segunda y tercera franja 

Figura 8 Registros artificiales de banda ancha (casos D4 y D5) 

En el caso que se consideren valores de t Q que no cumplan con la condicion (22) pueden 

obtenerse funciones envolventes del tipo que se muestra en la figura 9, para las dos franjas 
mencionadas. Si el valor del parametro es mayor que el limite superior, la funcion se hace muy 
estrecha en la punta y deja de ser del tipo gamma; la historia simulada con este caso se muestra 
en figura 10. 



Figura 9 Funcion envolvente obtenida cuando t Q no cumple la condicion (22) 

Analizando el caso de registros de banda angosta, en la Tabla 5 se muestran los valores 
medios y las desviaciones estandar de los parametros k x , k 2 y t 0 , mientras que en la tabla 7 se 

anotan los valores correspondientes a tres combinaciones de valores de las relaciones entre los 
momentos espectrales con las que se obtuvieron las funciones de modulacion de amplitud que se 
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detallan en la figura 11. Este ejemplo corresponde a un registro de magnitud 5.2 y una distancia 
epicentral de 15 km. Para la distribution de energia entre las franjas se adoptaron valores 
cercanos a los medios, esto es 70%, 20% y 10 % respectivamente. La figura 12 muestra los 
acelerogramas simulados segun las funciones de modulacion detalladas en la figura 11. Para el 
caso D los resultados son similares a los obtenidos con la regresion 1, por lo que se muestra solo 
uno de ellos. 

Las funciones obtenidas de aplicar en las formulas (19) los momentos temporales 
estimados directamente por regresion son las denominadas “Regresion 1” y “Regresion 2” en la 
figura 11. 



Tabla 7 Ejemplos de valores de k x y k 2 para sismos de banda angosta 



Franja 


Caso B 


CasoC 


Caso D 


K 


k 2 


K 


k 2 


K 


k 2 


1 


2.55 


9.06 


3.00 


14.50 


2.55 


7.00 


2 


2.90 


12.08 


3.00 


16.00 


2.80 


8.00 


3 


2.92 


11.59 


3.00 


17.00 


2.91 


8.00 



Figura 10 Registros artificiales de banda ancha (Caso D6) 



Los resultados obtenidos segun una u otra formula de regresion difieren notoriamente en 
este caso, lo que se detalla en las figuras 13 y 14, en las que se designa como regresion 1 las 
curvas envolventes obtenidas por (12) y como regresion 2 las calculadas segun (13) y (14). En el 
curso de las experimentaciones numericas llevadas a cabo durante el desarrollo de este trabajo, se 
observo que las mayores diferencias entre ambos modelos de regresion se dan al considerar como 
base de datos el conjunto denominado U1 y en especial para registros de banda angosta y de 
campo cercano; por otro lado, al realizar la regresion sobre el grupo U2 las funciones resultantes 
tienden a ganar estabilidad en su forma y a ser similares, tanto para registros de banda angosta 
como ancha. La figura 13 muestra un ejemplo de funciones envolventes para un sismo de banda 
angosta de magnitud 5, distancia epicentral de 10 km y profundidad de foco de 15 km mientras 
que la figura 14 hace lo propio para un sismo de magnitud 5.9, distancia epicentral de 30 km y 
profundidad de foco de 20 km. 



Figura 11 Funciones de modulacion para registros de banda angosta 
(Casos A, B, C y D) 



Figura 12 Registros artificiales de banda angosta 
Figura 13 Diferencias entre las funciones de amplitud resultantes de 
considerar U1 o U2 (M=5, A = 10 km, H=15 km) 



Figura 14 Diferencias entre las funciones de amplitud resultantes de 
considerar U1 o U2 (M=5.9, A = 30 km, H= 20 km) 



Se aclara que en cada ejemplo de este paragrafo se mantuvieron constantes todas las 
variables del problema salvo aquellas cuyo efecto se queria demostrar. En cada caso el proceso 
estacionario es el mismo, pues se mantuvieron constante los nucleos del generador de niimeros 
aleatorios. 



6.4 Diseno de la funcion de densidad espectral de potencia 
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La generacion de acelerogramas artificiales con este modelo de regresion puede realizarse 
utilizando como funcion espectral de potencia una funcion constante dada por la expresion (5) o la 
que incluye el espectro de fuente de Boore, expresada por (6). En general puede decirse que si se 
desean registros de banda angosta, como los sismos de tipo impulsivo por ejemplo, la funcion que 
incluye el espectro de fuente da mejores resultados, mientras que para los sismos de banda ancha, 
una distribucion constante para la densidad espectral es suficiente. 



6.5 Control de la energia 



La energia total del registro se determina por regresion lineal, en funcion de la magnitud 
y la distancia epicentral. Por medio de la division porcentual de esa energia en las franjas 
consideradas segun el tipo de registro que se desea, se determina la energia de cada banda y se 
disena la funcion de amplitud correspondiente en base a lo expuesto anteriormente. Una vez que 
el proceso de generacion es llevado a cabo banda a banda, la energia en cada una de ellas es 
controlada mediante la potencia (energia por unidad de tiempo) de la funcion de amplitud y la de 
la historia simulada correspondiente a esa banda frecuencial. Una vez realizada la correction si 
correspondiere, se suma la contribution de esta banda al registro total que se genera. Asi se logra 
que la energia del proceso generado es la que se obtuvo por regresion lineal. 



6.6 Duracion total y efectiva del registro generado 



Con respecto a la estimation de la duracion del registro que se busca generar, se pueden 
considerar distintos criterios. Asi por ejemplo, Crempien y Araneda (1993) hacen una regresion 
lineal sobre la duracion de la fase fuerte del movimiento, con la misma formula de regresion 
usada para los momentos. Hay varios aspectos que influyen en la duracion del registro en un 
determinado sitio, la que depende de factores como la frecuencia y la geomorfologia del lugar 
entre otros (Trifunac y Novikova, 1994) 

En este trabajo estos parametros tambien se determinan por regresion lineal, 
considerando las dos formulas propuestas. Los valores de las constantes obtenidos por medio de 
la regresion 1 se muestran en las tablas 3a y 3b para ambos grupos de datos. Por otro lado, en la 
tabla 4 se detallan los resultados de la regresion realizada segun las formulas del modelo 
alternative de regresion. 

Considerando los datos con que se cuenta, se juzga que esta simple metodologia dio 
resultados aceptables. De todos modos, es posible dar valores distintos a los calculados inferidos 
por necesidades de la aplicacion del registro generado (por ejemplo, para probar casos de 
licuefaccion de suelos se puede exigir sismos de mayor duracion total y efectiva), o por condiciones 
particulares del suelo y la topografia en el sitio de estudio. 

Una vez definida la duracion total del evento a simular, a los fines del algoritmo de 
simulation se debe definir el paso de discretizacion de la variable tiempo para determinar el 
numero de puntos del acelerograma artificial. La funcion temporal de amplitud no depende del 
intervalo de discretizacion At pero si hay diferencias en el campo frecuencial. Cuanto mas grande 
sea este paso de tiempo, el contenido espectral es menor y si la duracion total del acelerograma 
generado es pequena se perdera information. 

Como ejemplo, para simular acelerogramas con una duracion de 5 a 6 segundos, un paso 
de discretizacion temporal de 0.005 segundos dio buenos resultados; en el caso de registros de 10 
segundos de duracion total, un paso de 0.01 segundos dio resultados satisfactorios. Puede 
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sugerirse que si la duracion total es mayor a 20 segundos, un intervalo de tiempo de 0.02 
segundos es lo aconsejable. 



6. 7 Aceleracion maxima del registro generado 



Finalmente, con respecto a la aceleracion maxima deseada, esta puede ser definida en 
base a la obtenida por alguna formula de regresion, por condiciones geomorfologicas, exigencias 
de algun reglamento, etc. En general por este metodo se obtienen aceleraciones cercanas a las 
exigidas, dependiendo del buen diseno de la funcion de amplitud principal merited se admite una 
ligera correccion del pico, siempre y cuando esto no implique deformaciones inadmisibles de la 
forma del registro. 



7. EJEMPLO DE APLICACION 



Con los conceptos antes expuestos se ha confeccionado un programa de computadora 
(PREDICE) con el que se ha realizado una amplia experimentacion numerica aplicada 
principalmente a la zona en estudio, dando resultados satisfactorios tanto para sismos moderados 
como intensos. Se muestra un ejemplo que tiene como objetivo simular acelerogramas de un 
evento similar al ocurrido el 24 de Junio de 1984 en Granada (Espana). No es el principal 
proposito de este modelo generar acelerogramas concordantes con alguno tornado como muestra, 
pero se desea analizar los alcances de la aplicacion del modelo y ver de que manera es posible 
controlar sus parametros a fin de obtener historias con caracteristicas determinadas a priori. El 
espectro de Boore fue estimado con los siguientes valores ( Morales et al., 1996) 



Figura 15 Registros SFLO y SIM1 



Magnitude 5 

Distancia epicentraf 40.4 km 
Profundidad: 12 km 
Momento sismico : 3.0 10 16 N-m 
Caida de tension^ 0.100 10 8 N/m 2 
Frecuencia maxima^ 25 Hz 

Coeficiente de atenuacion dependiente de la frecuencia : 80 f° 65 



Figura 16 Espectros de Fourier de SFLO y SIM1 



La historia temporal de las aceleraciones registradas en Santa Fe corresponde a un 
evento de banda angosta por lo que se procedio a dividir el espectro en tres franjas frecuenciales, 
con un contenido energetico en cada franja de 85%, 10% y 5% respectivamente. En cuanto a las 
funciones de amplitud que rnodula el proceso estacionario en cada franja dieron buenos resultados 
las disenadas tanto en base a parametros determinados por regresion, como para relaciones entre 
momentos espectrales correspondientes al caso B de la tabla 7 y valores del parametro t Q de 
0.06, 0.05 y 0.10 segundos para cada franja respectivamente. 



Figura 17 Comparacion entre SFLO y SIM5 

En las figuras 15 a 19 se observan los resultados obtenidos para dos casos de seis 
simulaciones efectuadas. Se muestran las historias simuladas y sus respectivos espectros de 
amplitudes de Fourier, cotejandolas con el registro y espectro de Fourier originales. Las distintas 
simulaciones resultan de considerar distintos nucleos para el calculo de los procesos aleatorios 
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correspondientes. La figura 19 muestra los espectros de respuesta elasticos para otras tres 
simulaciones y puede observarse una buena correspondencia con el del registro original. 



Figura 18 Espectros de Fourier de SFLO y SIM5 



Figura 19 Espectros de respuesta elastica de SFLO y tres simulaciones 



8. CONCLUSIONES 



En cuento a las bases conceptuales sobre las que se define este modelo, se contempla la 
no estacionariedad temporal y frecuencial, mediante la aplicacion del concepto de espectro 
evolutivo donde la evolucion del contenido frecuencial en el tiernpo se materializa mediante una 
discretizacion del proceso en el campo frecuencial, considerando luego en el campo del tiernpo a 
cada franja componente corno un proceso estocastico uniformemente modulado, metodologia 
utilizada previamente por varios autores y aplicada a zonas que cuentan al rnenos con un registro. 

La determinacion de los parametros de la funcion de amplitud que rnodula cada proceso 
estacionario en las distintas franjas se los determina en funcion de los momentos temporales del 
registro a simular, obtenidos por medio de la tecnica de regresion multilineal en funcion de dos 
magnitudes macrosismicas de uso corriente en sismologia : magnitud y distancia epicentral. 

En cuanto a la base de datos utilizada en la regresion consta de registros cuyas 
magnitudes se ubican en el intervalo [4,7.4] los que fueron divididos en dos grupos : 

• El primero, al que se denomina Ul, contiene eventos de magnitudes comprendidas entre 4 y 6 
grados relacionadas a la frecuente actividad sismica moderada de la zona cuyana. Amen de 
los registros del sismo de Mendoza'85 se incluyen otros considerados del mismo tipo que los 
del evento mencionado, esto es sismos moderados, superficiales y resultantes de 
deslizamientos de fallas debido mayormente a un proceso tensional de compresion. 

• Para incluir los dos acelerogramas del sismo de Caucete, de magnitud 7.4, se anadieron 
algunos acelerogramas de sismos superficiales intensos del Oeste Americano, con lo que se 
conforma el conjunto U2 de datos. Ademas, teniendo en consideracion que la ciudad de 
Mendoza y sus alrededores esta asentada sobre varias fallas activas muy superficiales, se 
juzga factible que ocurran sismos de este tipo en la region bajo estudio, habida cuenta de que 
los informes de microzonificacion realizadas en esta area hablan de similitudes geotectonicas 
entre arnbas regiones. 

• Se comprobo que tanto la aceleracion maxima corno la energia de los registros de la base de 
datos utilizada siguen una distribution similar con respecto a la distancia epicentral y la 
magnitud del sismo. Por ello en primer lugar se hizo la regresion multilineal en base a la 
formula de atenuacion de aceleraciones maximas determinada por profesionales del Instituto 
Nacional de Prevencion Sismica de la Argentina, con la que se obtiene los valores de /t () , y 

X, . La misma expresion fue aplicada para determinar las duraciones efectiva t y duration 
total t j de los registros a simular, obteniendose buenos resultados. 

• Analizando las propiedades estadisticas de los resultados de la regresion se concluye que con 
esta formula la constante que multiplica a la distancia epicentral tiene en general poca 
significancia para explicar una relacion con los parametros que se desean estimar; se infiere 
que en el caso de los momentos de orden 1 y 2 y en especial para el conjunto de datos Ul solo 
la magnitud tiene relevancia en la definicion de estos parametros. 
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• A fin de mejorar estos aspectos se propone una segunda formula de regresion o mas 
propiamente dos formulas; primeramente se realiza una regresion para determinar el 
momento de orden 0 (energia del registro) y las duraciones efectiva y total, los parametros 
que estadisticamente tanto la magnitud como la distancia epicentral explican 
significativamente con las formulas de atenuacion. Los momentos de orden 1 y 2 se 
determinan tambien por regresion, mediante una expresion que es funcion de la energia y la 
duration efectiva determinadas por regresion. 

• Con estas expresiones el modelo (denominado de “regresion 2”) presenta mejores niveles en 
las pruebas de confidencia y significancia estadisticas que el modelo basado en las formulas 
de atenuacion existentes para la zona bajo estudio (identificado como modelo de “regresion 
1 ”). 

• Se realizaron simulaciones utilizando ambos conjuntos de datos y con ambas formulas de 
regresion. Se detectaron diferencias entre las funciones de amplitud logradas con uno u otro 
rnetodo al utilizar el grupo U1 de datos, sob re todo para registros de campo cercano. A1 
ampliar la base de datos y hacer uso del grupo U2, se observa mayor estabilidad en la formas 
de las funciones como asi tambien una reduction apreciable de las diferencias en los 
resultados segun ambos metodos, tanto para eventos de campo cercano como lejano. 

• Las regresiones se realizaron sobre los datos sin hacer especiales correcciones ni filtrados de 
ruidos o fallas de muestreo. Por ello, los momentos de orden 1 y 2 resultantes conllevan estos 
efectos que como se explico antes, dependen del intervalo considerado antes del primer 
impulso de la serial, esto es del parametro t Q y el ruido de la serial (parametro S ). 

• En este caso, el modelo no cuantifica automaticamente estos parametros por lo que las 
funciones de amplitud resultantes pueden ser defectuosas y no describir apropiadamente la 
variation temporal de la serial. Se determinaron formulas explicitas para corregir los 
momentos temporales de orden 1 y 2 a fin de incluir estos efectos de manera simple y clara en 
la determination de los misrnos. 

• Esta correction se impone en caso de generar sismos impulsivos de banda angosta y no es tan 
critica en el caso de eventos de banda ancha. 

Una vez determinados los momentos temporales totales, se calculan los misrnos para cada 
franja como fracciones de los respectivos momentos totales, considerandose los misrnos 
porcentajes del contenido energetico de la franja con respecto a la energia total del evento. 

Cabe tambien en el modelo la posibilidad de manipular los momentos de orden 1 y 2 a fin 
de permitir la construction de historias probables de ocurrir en un lugar determinado, con 
caracteristicas especificas de acuerdo a la aplicacion que se hara de las historias simuladas, 
manteniendo el momento de orden 0 (que representa la energia de la franja) inalterable. 

• Se proponen a tal fin los dos parametros y k 7 , que simplemente relacionan los momentos 
espectrales de orden 1 y 2 con el momento de orden 0 del registro, respectivamente. 

• En este caso, el momento temporal determinado por regresion es el de orden 0 y en funcion de 
este se calculan los de orden superior, en funcion de k l y k 2 . Se ha pensado esta alternativa 
para controlar de alguna manera el diseno de las funciones de amplitud en cada franja, 
variando el modo en que la energia se distribuye, pero no su intensidad. 

• A fin de asegurarse que la funcion de amplitud sea del tipo gamma los parametros k l , k 2 y 
t Q deben satisfacer ciertas condiciones que se explicitan en el trabajo. 
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• Se especifican los valores medios de los parametros k x , k 2 y t 0 obtenidos para sismos 

cuyanos, obtenidos de la experimentacion numerica llevada a cabo con el modelo de ajuste no 
lineal; sus desviaciones estandar denotan una gran variabilidad de los mismos. Los valores 
explicitados en el trabajo son solo indicativos y pueden ser tornados corno punto de partida 
para definir los que se ajusten efectivamente al problema a resolver en cada caso. 

• Se ha experimentado con varias combinaciones entre estos parametros, de lo que resultaron 
algunas reglas simples que se pueden sugerir a fin de disenar funciones de caracteristicas 
especificas determinada. Si se desean funciones gamma angostas se debe considerar valores 
bajos de k 2 y si se desean que las funciones gamma sean simetricas, se debe considerar un 
valor de k t cercano al limite superior del intervalo de valores que resulta de aplicar los 
requisites de funcion gamma. 

Con los conceptos antes expuestos, se ha confeccionado un programa de computadora al 
que se denomino PREDICE, en el que se requieren cuantificar los siguientes parametros que 
definen las caracteristicas de la serial generada: 

Numero de franjas frecuenciales para registros de la zona cuyana se sugiere considerar tres 
franjas para sismos de banda estrecha y cinco para sismos de banda ancha. En el caso de sismos 
cercanos (menos de 30 km de distancia epicentral) conviene optar por registros de banda angosta 
y para sismos lejanos es apropiado considerar registros de banda ancha. Como tambien las 
condiciones locales del suelo determinan el contenido espectral, se pueden tener en cuenta estas 
particularidades y por ende hacer otras consideraciones. 

Contenido energetico de las franjas frecuenciales la distribucion energetica en cada franja 
frecuencial para sismos moderados, cercanos y producidos por deslizamiento de fallas 
superficiales pueden ser: 68 %, 28% y 4% para sismos de banda estrecha (division en tres franjas) 
y del 4%, 23%, 51%, 19% y 3% para sismos de banda ancha (division en cinco franjas 
frecuenciales). Estos son los valores promedios de la base de datos usada pero pueden cambiarse 
a voluntad segun la necesidad de la simulacion. 

Diseno de la funcion de amplitud en cada franja, que puede realizarse mediante la aplicacion de 
formulas cerradas para la determinacion de los parametros de la funcion gamma tradicional en 
funcion de los tres momentos espectrales determinados por regresion lineal o bien en funcion de 
los dos parametros propuestos k t y k 2 que simplemente relacionan los momentos espectrales de 
orden 1 y 2 con el momento de orden 0 del registro, respectivamente. 

Eleccion de la funcion de densidad espectral de potencia, en cada franja uniformemente 
modulada ; puede ser considerada constante o en funcion del espectro de Boore. 

La aceleracion maxima, duracion efectiva y duracion total, se determinan mediante formulas de 
regresion con los parametros determinados para la base de datos usada o bien por consideraciones 
especificas del lugar de estudio y fines de la simulacion. En cuanto a la duracion de los sismos 
generados, los parametros determinados por regresion lineal condujeron a resultados 
satisfactorios en todos los casos analizados; de todos modos, se pueden determinar a priori y por 
otras consideraciones, valores distintos a los determinados por regresion, que corno se acoto 
repetidas veces, resumen caracteristicas de la base de datos utilizada en este trabajo. 

• En cada franja se controla que la energia del proceso uniformemente modulado generado 
tenga tenga minimas diferencias con respecto a la energia asignada a la banda frecuencial. 
Una vez obtenido el registro artificial se controla que la energia total del mismo sea similar a 
la obtenida por regresion. En general, se obtienen valores picos del orden de los determinados 
por la formula de atenuacion para la zona cuyana para suelo duro; el programa admite la 
posibilidad de una correccion que no altere las caracteristicas dinamicas basicas de la historia 
generada. 
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• Una vez generado el registro artificial se procede a su correction de linea base, tal como se 
realiza con los registros reales, mediante el clasico procedimiento de Berg y Housner. 

Se realizo una amplia experimentation numerica, aplicando el modelo principalmente a la 
zona bajo estudio. Dio resultados satisfactorios para predecir posibles movimientos en un 
determinado sitio, tanto moderados como intensos. Debido a la exigua base de datos con que se 
contaba al momento de realizar esta investigation, la mayor confiablidad corresponde a sismos 
moderados en el area de la capital mendocina y del Gran Mendoza, en la Republica Argentina. Se 
considera de todos modos que la metodologia propuesta es valida para cualquier sitio, con solo 
considerar los datos apropiados. 
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RESUMEN 



Con el apoyo financiero de la Direction Academica de Investigation (DAI) de la Pontificia 
Universidad Catolica del Peru (PUCP), se efectuo en el Laboratorio de Estructuras de la PUCP, 
un estudio experimental que tuvo por objeto impedir la falla por deslizamiento en los muros de 
albanileria armada construidos con bloques de concreto vibrado, cuando ellos se ven sujetos a 
terremotos. En primer lugar, se analizaron las posibles causas que pueden originar la falla por 
deslizamiento, para luego plantear soluciones economicas y sencillas de ejecutar en obra. 

Posteriormente, estas soluciones se aplicaron en la construction de cinco muros, los que 
luego fueron ensayados a carga lateral ciclica obteniendose resultados satisfactorios. 



1. INTRODUCTION 



En los edificios estructurados con muros de albanileria armada sometidos a terremotos, 
puede ocurrir, por una serie de razones, una falla por deslizamiento (cizalle o corte ’friction). Esta 
forma de falla es muy peligrosa, porque la section transversal del muro se reduce 
significativamente figura 1, debido a la action simultanea de las cargas sismicas coplanares y 
perpendiculares al piano del muro, y tambien, porque esta falla puede causar incluso la rotura del 
refuerzo vertical y el disloque del muro por encima del piano de falla, lo que genera una doble 
excentricidad de la carga vertical. De esta manera, el objetivo del proyecto es tratar de eliminar la 
falla por deslizamiento, o al menos retardar su aparicion hasta que los muros alcancen la 
distorsion angular inelastica maxima permitida para las construcciones de albanileria armada 
por nuestra norma sismica (0.005, Ref.l) cuando ocurren sismos severos, transformando este tipo 
de falla en otras menos peligrosas, como son: las fallas por flexion o por fuerza cortante. 



2. POSIBLES CAUSAS DE LA FALLA POR DESLIZAMIENTO 

2.1. Falta de adherencia bloque-mortero. 



Los bloques tienen una suction relativamente elevada, lo que conduce a que el mortero 
colocado sobre el endurezca rapidamente, perdiendose la adherencia con el bloque de la hilada 
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inmediata superior, originandose de este modoun piano potencial de deslizamiento (figura.2). 
Cabe indicar que los bloques se asientan en su estado natural (secos) y el unico tratamiento que 
reciben es la limpieza de particulas sueltas. 




Figura 1 Falla por deslizamiento. 




Figura 2 



2.2. Traslape del refuerzo vertical. 



Tradicionalmente, para facilitar el proceso constructive de la albanileria armada, se 
acostumbra dejar espigas ancladas en la cimentacion, para despues de terminar con la 
construccion de la albanileria, insertar las barras verticales traslapandolas con las espigas, 
creandose de esta manera un piano de debilidad en la zona donde terminan las espigas, por donde 
puede ocurrir la falla por deslizamiento figura 3. 

Para evitar el problema que ocasionan los traslapes, en otros paises se fabrica bloques H 
figura 4, que permiten el uso de refuerzo vertical continuo, con lo cual se descongestionan las 
celdas reduciendose la posibilidad de que se formen cangrejeras en el concrete liquido (“grout” o 
material con que se rellenan los muros). El emplear refuerzo vertical continuo y bloques 
convencionales (no en forma de H), retarda el proceso de asentado durante la primera jornada de 
trabajo, por el tiempo que demora insertar y bajar al bloque a traves del refuerzo vertical. 




Figura 3 
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Figura 4 



2.3. Falta de adherencia grout-cimentacion. 



Puede ser causado por una serie de razones a) tratamiento inadecuado de la superficie del 
cimiento (rayado poco profundo, presencia de lechada superior, falta de limpieza al iniciar la 
construction de la albanileria y tambien antes de vaciar el grout); b) segregacion del grout figura5 
producida por: S u alto contenido de agua, compactacion inadecuada, altura excesiva de vaciado; 
etc. 




Figura 5. Segregacion del grout. 



2.4. Contraction de secado del grout. 



Lo ideal es que quede integrado como una sola unidad la union bloque-grout-mortero- 
acero. De separase el grout del bloque y del acero, el refuerzo perderia efectividad ante las 
acciones de flexion, cizalle, etc. El polvo de la arena y la cantidad de agua que requiere el grout, 
incrementan este fenomeno. 



2.5. Distribution uniforme del refuerzo vertical. 



Tradicionalmente, se distribuye el refuerzo vertical a lo largo del muro. Puesto que el 
refuerzo en los extremos del muro es pequeno, la base del muro rota con facilidad. Este hecho 
provoca el aplastamiento de los talones y el corrimiento de las fisuras horizontales por flexion, las 
que al encontrarse por la action ciclica de la carga sismica, terminan transformandose en una 
falla por deslizamiento figura 6. 
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Figura 6 Transformation de la falla por flexion en falla por deslizamiento. 



2.6. Uso de cintas de mortero longitudinales y verticales. 



El grout no rellena el espacio entre las cintas, lo que debilita la union bloque~mortero, con 
la consecuente formation de fisuras. 




Figura 7 Cintas de mortero en las juntas. 



3. SOLUCIONES PARA. EVITAR LA FALLA POR DESLIZAMIENTO 



En este proyecto se trata de encontrar soluciones al problema del deslizamiento, que sean 
econo micas y faciles de ejecutar en obra. Sin embargo, en este acapite se describen todas las 
posible soluciones, ordenadas en funcion del costo y de la dificultad constructiva, a fin de tener un 
panorama amplio que permita eliminar las causas del problema, indicadas en el acapite 2, y 
seleccionar las soluciones mas adecuadas. 

1. Profundizar el rayado de la superficie superior de la cimentacion (minimo 6 mm), limpiarla de 
particulas sueltas y de lechada, para luego humedecerla antes de asentar los bloques y antes 
de vaciar el grout. El objetivo de esta operation es mejorar la adherencia y la resistencia a 
corte'friccion en las interfases mortero'timentacion y grout-cimentacion. El rnuro adoptado 
como patron de comparacion (Ml, traditional) tuvo un rayado conventional. 

2. Manteniendo la cuantia total de refuerzo vertical, concentrar la mayor area de acero en los 
extremos del muro, empleando en la zona central una cuantia minima igual a 0.001. El 
objetivo es disminuir el giro por flexion para controlar el corrimiento de las fisuras 
horizontales. El muro patron (Ml) tuvo una distribution uniforme del refuerzo vertical. 
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3. Invertir la cara de asentado de los bloques correspondientes a la primera hilada, creando 
aberturas en las 3 tapas transversales, de manera que el grout al pasar a traves de esas 
aberturas, rellene una mayor area en su zona de contacto con la cimentacion. 




4. Compactar al grout empleando una madera con section transversal de l”x 2” en reemplazo de 
la varilla lisa de 5/8” que se utiliza tradicionalmente. El objetivo es tratar que el grout se 
expanda lateralmente en mayor grado al apisonarlo con la madera. Esta solution es aplicable 
solo cuando se utiliza bloques de 19 cm. de espesor, con celdas grandes. 

5. Disminuir la cantidad de confitillo en el grout, de una proportion cemento : arena : confitillo 1- 
2 54: 1 14 (traditional) a 1- 3 : 1, y tamizar la arena a traves de la malla #200 para eliminar el 
polvo. El objetivo es atenuar el problema de segregation y contraction de secado. 

6. Similar a la solution 5, pero empleando cemento puzolanico y reduciendo el slump de 10” 
(traditional) a 9”. Se tiene conocimiento que el cemento puzolanico atenua la contraction de 
secado, mientras que al disminuir la cantidad de agua (al bajar el slump), se reducen los 
problemas de segregation y contraction de secado. 

7. Similar a la solution 5, pero agregando 1/10 de cal para retardar la contraction de secado. 

8. Tratamiento con agua para incrementar la adherencia bloqueunortero y atenuar la 
contraction de secado del grout. Con una brocha, humedecer las juntas 2 horas despues de 
haber asentado los bloques, una vez al dia, hasta el dia del vaciado del grout. Antes del 
vaciado, regar internamente las celdas, encofrar las ventanas de limpieza, vaciar y compactar 
al grout, luego curar al muro regandolo externamente una vez al dia durante una semana. 

9. Adicionar espigas en las celdas que carecen de refuerzo vertical, ancladas en la cimentacion. 
El objetivo es incrementar la resistencia a corte -friction en las dos primeras hiladas del muro 
y en la interfase muro’timentacion, ademas, que las fisuras por flexion se originen en las 
partes mas altas del muro, de tal manera que la cimentacion no estorbe su trayectoria. 

10. Dejar orificios en la cimentacion, en las zonas ocupadas por las celdas. El objetivo es crear 
llaves de corte con el grout, para transferir la fuerza cortante hacia la cimentacion. Esta 
solution es menos efectiva que la numero 9, debido a que la falla de deslizamiento se puede 
trasladar a las hiladas superiores. Tambien, es dificil hacer los orificios y limpiarlos. 

11. Emplear refuerzo vertical continuo, con bloques recortados en forma de H para facilitar el 
asentado. El objetivo es tener mas espacio en las celdas para poder compactar al grout, 
evitando la posibilidad que se formen cangrejeras, asi como evitar los problemas que causan 
la discontinuidad del refuerzo vertical. Esta solution mejora al 100% las juntas verticales al 
quedar llenas de grout, pero debe anadirse otras soluciones para mejorar la union horizontal 
bloquemiortero, grout-cimentacion y atenuar la contraction y segregation del grout. 

12. Usar grout fino (cemento-arena-cal V 3 : %) en la mitad inferior del muro y grout grueso 
(cemento-arena-confitillo-cal V 3 : V- 0.1) en la mitad superior. Con esta solution se evitaria el 
problema de segregation del grout en la base del muro, pero se incrementaria el problema de 
contraction de secado en la mitad inferior. Esta solution es de dificil aplicacion en obra, por 
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tener que variarse la mezcla en un mismo piso, sin embargo, no se descarta ya que el grout se 
vacia siempre en dos etapas, con un lapso de % hora entre ellas. 

13. Uso de aditivo en polvo Wacker (culminal mas vinnapas) en el mortero, en reemplazo de la 
cal. Este aditivo mejora notoriamente la adherencia bloqueunortero (37%, segun la referenda 
4), pero es costoso (segun la referenda 4, el costo de la albanileria se eleva en 18%). 

14. Uso de aditivo Intraplast (de Sika) en el grout para mejorar su plasticidad y disminuir su 
contraccion de fragua. En el mercado existen aditivos que permiten reducir la cantidad de 
agua en el grout, o superplastificantes (Rheobuild 1000 de MBT) para producir grout 
autocompactado (no necesitan vibradora ni chuzo). Pero, estas soluciones son costosas. 

15. Uso de vibradora de aguja de !4” para compactar al grout. Para esto se inserta apagada la 
vibradora a traves de las celdas de los bloques, se vacia el grout, se enciende la vibradora y se 
la retira lentamente, debiendose hacer esta operacion celda por celda. Esta solucion se 
descarta por ser poco practica. 

16. Low Lift Grouting. Consiste en construir la primera mitad del muro, llenar con grout al dia 
siguiente, al tercer dia construir la segunda mitad y al cuarto dia volver a llenar con grout. 
Esta solucion atenuaria la posibilidad de segregacion del grout, ya que el vaciado se hace 
desde una altura de 1.2 m, pero no se utiliza por la posibilidad de que se generen juntas frias 
entre las capas de grout y tambien porque el proceso constructive se retarda un dia. Cabe 
indicar que para limpiar los desperdicios de mortero que caen a traves de las celdas durante 
la construccion de la segunda mitad del muro, habria que abrir ventanas en los bloques 
ubicados a media altura del muro, esto no es posible porque este tipo de muro es caravista. 



4. PROBETAS PRELIMINARES DE GROUT 



El objetivo que se busca en este acapite es definir el grout a utilizar en los muros 
mejorados, tratando de reducir los problemas de contraccion de secado y de segregacion. Para esto 
se fabricaron 10 mezclas distintas (Gi en la figura 8), empleando como molde un pilar de bloques 
U (medios bloques), con una altura de 2.4 m (figura 10). Estos moldes permitieron la 
transferencia natural de agua desde el grout hacia los bloques. Por otro lado, se utilizo encofrados 
transparentes con el objeto de observar lo que ocurria con el grout durante su vaciado y 
compactacion. Al retirar el encofrado, pudo observarse si existian los problemas mencionados. 




Las mezclas (Gi) utilizadas fueron. 
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• Gl. Grout tradicional, cemento normal-arena normal-confitillo ; 1 ; 2 Vi- 1 Vi, slump 10”. 

• G2. Cemento normal-arena tamizada-confitillo : 1- 3 : 1, slump 10”. 

• G3. Similar a G2, pero con 0.1 de cal aerea normalizada. 

• G4. Similar a G2, pero con aditivo Intraplast. 

• G5. Similar a G2, pero con cemento puzolanico IP. 

• G6. Similar a G5, pero con 9” de slump. 

• G7. Similar a G6, pero con 0.1 de cal. 

• G8. Similar a G6, pero sin confitillo (cemento-arena 1- 3) en la mitad inferior y con confitillo 
(cemento-arena-confitillo V 3 : l) en la mitad superior. Slump 9”, sin cal. 

• G9. Similar a G8, pero con 14 de cal en la mitad inferior y 0.1 de cal en la mitad superior. 

• G10. Similar a G5, pero con tratamiento con agua (regado interno de celdas antes del 

vaciado). 

Externamente se pudo observar que el grout tradicional (Gl) tuvo los mayores problemas 
de segregacion y contraction de secado, mientras que estos problemas se redujeron en las mezclas 
G5, G6, G7 y G10 (con cemento puzolanico y arena tamizada), pero, las mezclas G6 y G7 (con 9” 
de slump), se notaron muy secas, lo que podria provocar cangrejeras internas en los muros, por lo 
que es conveniente utilizar 10” de slump. Para revisar internamente al grout, el pilar fue izado y 
volcado con un puente grua, para luego cortar su base con una amoladora a la altura de la 
primera hilada. De la base del pilar se extrajeron muestras de grout, con el objetivo de observar 
su estructura interna, cabe destacar que la mayor concentration de piedras (segregacion en la 
base) se encontro en la probeta Gl (figura 10). Como conclusion de estos ensayos preliminares, se 
utilizo en los muros mejorados: G2 para M2, G5 para M3, y G10 para M4 y M5. 




Figura 9 Preparation del grout y construction del pilar. 
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Figura 10 Extraccion de muestra de grout de la base del pilar y muestras G5 y Gl. 



5. CONSTRUCCION DE LOS MUROS 



En total se construyeron cinco muros a escala natural (ver las figuras 11 a 15), 

empleandose la tecnica constructiva tradicional para el muro Ml y manteniendose constante los 

siguientes parametros (la mayor parte de ellos se emplean tradicionalmente en obras reales). 

• Geometria del muro las dimensiones fueron 2.4 x 2.4 x 0.14 m. 

• Aparejo se utilizo un aparejo de soga (14 cm. de espesor). 

• Calidad del bloque se utilizo bloque tipo Bill, serie 15.14, procedente de la fabrica FIRTH, con 
ranuras en las tapas transversales que permitieron el paso del refuerzo horizontal. 

• Tratamiento de los bloques las superficies de asentado fueron limpiadas con una brocha de las 
particulas sueltas antes de asentarlos en su estado natural (secos). 

• Calidad del mortero cemento-cabarena R Vi- 4. Se uso arena gruesa tamizada entre las mallas 
ASTM #4 y #200. El cemento fue Portland tipo 1, no se uso cemento puzolanico porque 
proporcionaria un color marron a las juntas. 

• Calidad de la cal se utilizo cal aerea hidratada y normalizada de Agregados Calcareos. 

• Grosor de las juntas las juntas horizontales y verticales (brunadas) tuvieron 1 cm. de grosor. 

• Calidad del acero de refuerzo se utilizo acero corrugado de grado 60 de SIDERPERU. 

• Batido del grout se utilizo una mezcladora. El slump medido en el cono de Abrams fue 10”. 

• Construccion de la albanileria se construyo en 2 jornadas con una mano de obra calificada. 

• Tecnica de limpieza de las celdas los bloques de la primera hilada tuvieron ventanas de 
limpieza de 3”x4” en cada celda. Terminada cada jornada de trabajo, se limpio internamente 
las celdas con una varilla, sin tocar al muro. Antes del vaciado del grout, la base se rasqueteo 
con un tortol, las celdas se limpiaron con aire comprimido y se humedecio la cimentacion. 

• Tecnica de vaciado del grout fue vaciado en dos capas, cada una fue compactada con una 

varilla lisa de 5/8”, transcurrido 5 minutos, cada capa fue recompactada. El intervalo de 

vaciado entre capas fue Vi hora. El borde superior de cada capa quedo aproximadamente a 10 
cm. por debajo de la superficie superior correspondiente a la ultima hilada de cada jornada. 
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• Cuantia de refuerzo. El refuerzo horizontal tuvo una cuantia de 0.0025 (l f 3/8” @ 20 cm.) y la 
cuantia de refuerzo vertical total fue 0.0024. 

• Vigas. Las vigas de cimentacion y las vigas soleras tuvieron la misma seccion transversal, la 
misma calidad de concreto y el mismo refuerzo. 

Matriz de Variantes en la tabla 1, se indica las variantes que tuvieron los 5 muros, 
indicandose que el muro Ml, construido con la tecnica tradicional, corresponde al patron de 
comparacion, asimismo, las mejoras varian progresivamente, desde M2 hasta M5. 



Tabla 1. Matriz de Variantes. 



VARIANTE 


Ml 


M2 


M3 


M4 


M5 


Rayado de la 
cimentacion 


Normal 


Profundo 


Profundo 


Profundo 


Profundo 


Distribucion del 
refuerzo 

vertical 


Distribuido 


Concentrado 


Concentrado 


Concentrado 


Concentrado 


Tamizado de la arena 
del grout 

por mallas #4 y #200 


no 


si 


si 


si 


si 


Tipo de junta vertical 


Cintas 


Llena 


Llena 


Llena 


Cintas 


Traslape del refuerzo 
vertical 


100% en la 

misma 

seccion 


alterno 
h = 70 cm. 
y 95 cm. 


alterno 
h = 70 cm. 
y 95 cm. 


alterno 
h = 70 cm. 
y 95 cm. 


continuo, 

sin 

traslapes 


Grout 

cemento'arena - 

confitillo 


l: 2 y 2 : i y 2 


l: 3: 1 


l: 3: 1 


l: 3: 1 


l: 3: 1 


Tipo de cemento en el 
grout 


Portland 


Portland 


puzolanico 


puzolanico 


puzolanico 


Tratamiento con 
agua en las 
juntas, celdas 
internas y muro 


no 


no 


no 


si 


si 


Espigas adicionales 


no 


no 


no 


si 


si 


Bloque de hilada 1 


normal 


normal 


invertido 


invertido 


H 





Figura 11 . Caracteristicas comunes de los 5 muros. 
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Figura 12. Refuerzo vertical en Ml y M2. 





Figural3 Refuerzo vertical, bloque de la primera hilada y espigas adicionales en M3, M4 y M5. 
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Figural4. Construccion de muros Ml (izquierda), M3 (centro) y M5 (derecha). 




Figura 15. Tratamiento con agua en M4 y M5. Ver el parrafo 3.8. 



6. ENSAYOS DE CONTROL 



6.1. Grout 



Empleando como moldes a los bloques forrados con papel filtro figura 16, para permitir la 
transferencia de agua del grout hacia los bloques, se fabricaron 9 probetas de 10x10x19 cm., tres 
provinieron del muro Ml (tradicional), las otras tres de M2, mientras que las tres restantes 
fueron hechas con el grout de M3, M4 y M5. Estas probetas no fueron curadas, y luego de 28 dias 
fueron ensayadas a compresion, obteniendose los siguientes resultados promedios. 
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• Grout V 2 54 : 1 14 (cemento Portland tipo P muro Ml): f'c = 240 kg/cm 2 

• Grout l: 3 -1 (cemento Portland tipo I - muro M2): f'c = 193 kg/cm 2 

• Grout l: 3 : 1 (cemento puzolanico tipo IP- muros M3, M4 y M5) : f'c = 197 Kg./cm. 2 

Para los tres casos, la resistencia supero al valor minimo especificado por la Ref.2 (140 Kg./cm 2 ). 

6.2. Prismas de albanileria simple (no reforzada) 



Se construyeron un total de 6 pilas y 6 muretes, los cuales fueron ensayados a compresion 
axial y a compresion diagonal, respectivamente. Las pilas (14x39x60 cm.) estuvieron compuestas 
por 3 hiladas y los muretes por 4 hiladas (14x80x80 cm.). Tres pilas y tres muretes fueron 
rellenados con el grout tradicional de Ml (l: 2 54: 1 54), mientras que el resto de prismas fueron 
rellenados con grout similar al empleado en M4 y M5 (l ; 3 : 1, con cemento puzolanico), ademas de 
recibir el tratamiento humedo especificado en 3.8 para M4 y M5. En la tabla 2 se indica los 
resultados, mientras que en la figura 17 se aprecia que los muretes rellenos con grout tradicional, 
tuvieron una falla escalonada a traves de las juntas, en tanto que los rellenados con grout 1:3:1, 
tuvieron algunas grietas que atravesaron los bloques, lo que implica una mejora en la adherencia 
mortero-bloque, producto del tratamiento que tuvieron esos prismas, sin embargo, el grout de esos 
prismas se noto que aun estaba humedo, pese a que tenia 28 dias de edad, lo que debe haber 
repercutido en los resultados, principalmente en la resistencia a compresion de las pilas. 



Tabla 2. Resistencia de los Prismas de Albanileria Simple 



Prisma 


Resistencia 


Grout l: 2 54 '■ 154 Cemento I 


Grout 1 ; 3 ; 1 Cemento IP 


Pilas 


Promedio fm 


104 kg/cm 2 


84 kg/cm 2 


Dispersion 


3.0 % 


6.0 % 


Caracteristica f m 


101 kg/cm 2 


79 kg/cm 2 


Muretes 


Promedio vm 


13.1 kg/cm 2 


12.5 kg/cm 2 


Dispersion 


8.6 % 


3.5 % 


Caracteristica v'm 


12.0 kg/cm 2 


12.1 kg/cm 2 




Figura 17. Falla escalonada en muretes rellenos con gout i: 2 54 : 1 14 (cemento tipo I, izquierda) y falla 
mixta en muretes rellenos con grout i: 3 : 1 (cemento IP, derecha). 
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7. TECNICA DE ENSAYO DE CARGA LATERAL CICLICA 



Los cinco muros fueron ensayados a carga lateral ciclica figura 18, controlando el 
desplazamiento horizontal de la viga solera en 10 fases de distinta amplitud tabla 3 y figura 19. 
La carga lateral (V) se aplico mediante un actuador dinamico, mientras que los desplazamientos 
fueron medidos con LVDT (D en la figura 18). Las fisuras de cada fase se pintaron de distinto 
color. 



Tabla 3. Fases del Ensayo de Carga Lateral Ciclica. 



Fase 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


D1 

(mm) 


0.25 


1.00 


2.50 


5.00 


7.50 


10.00 


12.50 


15.00 


17.50 


20.00 


No. 

Ciclos 


1 


2 


3 


3 


3 


3 


4 


4 


4 


4 


Color 


-- 


rosa 


lila 


negro 


Rojo 


azul 


marron 


verde 


naranja 


celeste 




Figura 18. Vista global de los dispositivos utilizados en el ensayo de carga lateral ciclica. 

8. COMPORTAMIENTO DE LOS MUROS 



A continuacion se describe el comportamiento que tuvieron los muros durante el ensayo de 
carga lateral ciclica (figura 20), poniendo enfasis en el problema de deslizamiento. Cabe indicar 
que la norma sismica (referenda l) permite una distorsion inelastica maxima de 0.005 para los 
entrepisos de los edificios de albanileria armada, que en este caso se produce en la fase 7 del 
ensayo (distorsion = Dl/h = 12.5 mm/ 2500 mm = 0.005, figura 19‘derecha), por lo que si la falla 
por deslizamiento se produjera en una fase mayor, careceria de importancia; sin embargo, debido 
a que en el experimento no se contemplo una serie de factores reales, como la accion dinamica con 
una gran cantidad de ciclos que imprime un sismo real, se considero prudente tomar como limite 
maximo un desplazamiento D1 = 15 mm (fase 8), asociado a una distorsion de 0.006. 
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PI Ml 



01 (mat 




D1 (mm) 

Figura 19 Lazos histereticos V-Dl de los muros Ml y M4 (arriba) y envolvente V'Dl del 

primer ciclo para los 5 muros (derecha). 



En el muro Ml se presentaron 2 pianos de deslizamiento, el primero estuvo ubicado en la 
base y se produjo en la fase 5, mientras que el segundo se ubico en la base de la quinta hilada, 
donde terminaban las espigas, esta falla se presento en la fase 7. En M2, el deslizamiento se 
produjo a traves de la base durante la fase 5. En M3, la falla por deslizamiento en la base se inicio 
en la fase 7, en tanto que M4 no tuvo problemas de deslizamiento. En M5, la falla por 
deslizamiento se produjo durante la fase 9 en la base de la tercera hilada, donde terminaban las 
espigas. 



8.1. Observaciones Post-Ensayo 



Terminado el experimento, se procedio a inspeccionar internamente el estado de cada 
muro (figuras 21 y 22). Esta operacion se hizo con el objetivo de detectar cangrejeras, segregacion 
del grout, el grado de adherencia mortero-bloque-grout-acero, el estado del refuerzo, etc. 
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Figura 20 Fase 10 D1 = 20 mm 
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Figura 22 Pandeo del refuerzo vertical por deslizamiento en la base de M3 y eliminacion del 
pandeo al evitarse el deslizamiento en M4. 



9. CONCLUSIONES 

El objetivo principal de este proyecto, que fue controlar o eliminar la falla por 
deslizamiento transformandola en otra menos peligrosa como es la falla por flexion o por corte, 
pudo lograrse. Para ello es imprescindible el uso de espigas de 3/8”, colocandolas en aquellas 
celdas donde no exista el refuerzo vertical calculado por flexion, y que estas espigas tengan 
distintas longitudes por encima de la cimentacion (se propone 30 y 50 cm.), intercalandolas cada 
20 cm., que es el espaciamiento entre las celdas del bloque, y anclandolas 30 cm. en la 
cimentacion. 

Cabe destacar que ninguna mejora que se haga en este tipo de muro sera efectiva o se 
traducira en incrementos de resistencia, si es que no se evita la falla por deslizamiento; por 
ejemplo, de nada servira incrementar el refuerzo horizontal, o la adherencia bloque-grout-acero, o 
la eliminacion de cangrejeras internas (como en M5), si no se elimina la falla por deslizamiento. 



9.1. Proceso Constructive 

1. El rayado profundo de la superficie de la cimentacion (M2 a M5), mejoro la union mortero- 
cimentacion. Sin embargo, en la base del muro las fisuras horizontales se trasladaron a la 
union bloque-mortero, por la baja adherencia que presentaron esos materiales en todas las 
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hiladas. De requerirse mayor resistencia sismica, el problema de adherencia bloque-mortero 
puede solucionarse (referencia 4) agregando al mortero el aditivo en polvo Wacker. 

2. El rayado profundo de la superficie de la cimentacion y el empleo en la primera hilada de 
bloques recortados en sus tapas transversales para que circule el grout ocupando una mayor 
area, no incremento la resistencia a corte-friccion en la base del muro, debido a que la base 
del grout fue cubierta por una capa lisa de lechada de cemento proveniente del mismo grout. 




3. La disminucion del confitillo en el grout, manteniendo el slump en 10”, de una proporcion 
tradicional cemento-arena-confitillo V 2 W- 114 (usada en Ml), a una dosificacion V 3- 1, 
disminuyo el problema de segregacion del confitillo en la base del grout. 




4. El tratamiento con agua hecho a los muros M4 y M5 (curado de juntas con brocha, asi como el 
regado interno antes de vaciar el grout y regado externo despues del vaciado) mejoro 
ligeramente la adherencia bloque-grout, al atenuar la contraccion de secado del grout. 

5. El empleo de cemento puzolanico IP en el grout del muro M3, redujo ligeramente el problema 
de contraccion de secado del grout, observandose una mejor adherencia bloque-grout que la 
existente en los muros Ml y M2, donde se utilizo cemento Portland tipo I. Esta mejora no se 
reflejo en incrementos de resistencia porque M3 fallo por deslizamiento. 

6. En este proyecto no se presentaron cangrejeras en el grout en la zona de traslapes de las 
barras verticales, debido a que esas barras tuvieron poco diametro (%”), por lo que no 
causaron mayor congestion en las celdas de los bloques. 

7. Todas las barras horizontales quedaron desplazadas hacia uno de los lados del muro, debido a 
que se instalaron sueltas y se movieron durante el vaciado y compactacion del grout. Esto 
podria causar con el tiempo problemas de corrosion y falta de adherencia con el grout. Para 
solucionar este problema se sugiere utilizar barras verticales continuas, con lo cual, las barras 
horizontales pueden amarrarse con la verticales, esta solucion exigiria el empleo de bloques H 
para facilitar el asentado, o que las fabricas produzcan bloques con menor separation entre 
las ranuras existentes en sus tapas transversales (lugar donde se aloja el refuerzo horizontal). 

8. El espacio existente entre las cintas verticales de mortero, empleadas en el muro Ml 
(tradicional), no fue rellenado por el grout, lo que dio lugar a una zona de debilidad. En el 
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resto de muros, se utilizo junta vertical llena, pero esta mejora no se reflejo en la resistencia 
de los muros (excepto en M4), debido a que fallaron por deslizamiento. 

9. El empleo de arena seca colocada en la base del muro durante su construction, para evitar 
que los desperdicios de mortero queden pegados sobre la cimentacion formando juntas frias en 
la interfase grout-cimentacion, deberia descartarse, debido a que la arena seca absorbe el 
agua del mortero debilitando en la base del muro las uniones bloquemiortero y mortero - 
cimentacion. 

En reemplazo de la arena, es conveniente utilizar retazos de plasticos. 



9.2. Espigas Verticales 



En este proyecto se utilizaron dos tipos de espigas verticales ancladas en la cimentacion y 
son l) Las necesarias para traslapar las barras verticales requeridas por flexion (tipo A); y, 2) las 
adicionales para incrementar la resistencia a corte-friccion en la base del muro (tipo B, de 3/8” de 
diametro y 30 cm. de longitud), obteniendose las siguientes conclusiones. 

1. El empleo de espigas tradicionales tipo A de la misma longitud (barras traslapadas en la 
misma section transversal), causo en el muro Ml dos pianos de deslizamiento : en la base y 
donde terminan las espigas, lo que resulta muy peligroso, especialmente cuando el muro este 
sujeto a cargas sismicas perpendiculares a su piano, en vista que se comportaria como una 
barra triarticulada (inestable). Incluso, los bloques ubicados en los extremos de Ml llegaron a 
cizallarse horizontalmente en su parte intermedia, en la zona donde terminaban las espigas. 

2. El empleo de espigas tipo A con distintas longitudes (barras traslapadas en distintas 
secciones transversales) ubicadas alternadamente, evito en los muros M2 y M3 el segundo 
piano de deslizamiento; en estos muros la falla por deslizamiento se concentro en la base. 

Cabe destacar que en M3 la falla por deslizamiento se initio a partir de la fase 7 del 
ensayo (D1 = 12.5 mm), coincidente con la distorsion maxima permitida por la norma E - 030 para 
sismos severos, mientras que en Ml el deslizamiento se initio en una etapa temprana del ensayo 
(fase 5, D1 = 7.5 mm), es decir, en muros reales como el Ml, la falla por deslizamiento podria 
presentarse ante la action de sismos moderados. 

3. En el muro M4, donde se utilizo espigas tipo A con distintas longitudes y se adiciono espigas 
tipo B en aquellas celdas que carecian de refuerzo vertical, el comportamiento sismico mejoro 
notablemente. En este muro no se produjo la falla por deslizamiento, sino mas bien una falla 
mixta por flexion y corte, elevandose la resistencia maxima en 21% con relation a la 
alcanzada por Ml. Cabe destacar que en M4, no se produjo el pandeo de las barras verticales 
extremas, a diferencia del resto de muros que tuvieron falla por deslizamiento. 

4. En el muro M5, donde se utilizo barras verticales continuas y espigas tipo B en aquellas 
celdas que carecian de refuerzo vertical, se evito la falla por deslizamiento en la base, pero 
esta se traslado a la base de la tercera hilada, donde terminaban las espigas, por lo que es 
necesario que estas espigas tengan distintas longitudes. Cabe destacar que a diferencia de 
M4, el muro M5 carecio de las espigas tipo A, por lo que en la zona de traslape de las barras 
verticales, M5 tuvo la mitad de refuerzo vertical que es el elemento que controla la falla por 
deslizamiento. 

Adicionalmente, el deslizamiento se initio despues que los talones se trituraron durante la 
fase 9, y la trituration se debio a que el grout aun estaba humedo producto del tratamiento con 
agua que tuvo M5, pese a que habian transcurrido 28 dias hasta el instante del ensayo. 
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RESUMEN 



Se presenta un modelo matematico para identificar las matrices dinamicas de rigidez y 
flexibilidad de estructuras reticulares a partir de datos de mediciones de vibracion forzada. En el 
planteamiento se considera que el comportamiento del sistema estructural es lineal, elastico y no 
amortiguado. 

La definition de las matrices dinamicas de rigidez y flexibilidad se consigue al identificar 
mediante mediciones de vibracion: la matriz modal, las frecuencias naturales de vibracion y la 
frecuencia de excitation. Es necesario proponer la matriz de masas del sistema por no haber sido 
posible deducirla experimentalmente en esta etapa de la investigation. 

Se muestra la aplicacion de la metodologla desarrollada a un modelo estructural sencillo y 
se comentan los resultados obtenidos. 

Palabras Clave : identification, mediciones, vibracion, rigidez, flexibilidad 



ABSTRACT 



A mathematical model to identify the dynamic matrix of stiffness and flexibility for 
reticular structures from measurement data of forced vibration is presented. The model considers 
that the system is linear, elastic and undamped. 

The dynamic stiffness and the flexibility matrix are defined once the following data is 
obtained by means of vibration measurements: the modal matrix, the natural vibration 
frequencies, and the frequency of the driving force. It is necessary to propose the structure’s mass 
matrix since it is not possible to deduce it experimentally in this stage of the research. Some 
results for a simple structural model are included and discussed. 

Keywords: identification, measurements, vibration, stiffness and flexibility 



1. INTRODU C CION 



En general los sistemas estructurales estan sometidos a cargas que varian en el tiempo 
tal como las que se originan por viento, oleaje, sismos y maquinaria. El efecto de estas cargas 
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sobre la estructura origina una respuesta dinamica que es de suma importancia para el ingeniero 
disenador debido a que la respuesta dinamica puede exceder a la respuesta estatica en un factor 
significante. 

El gran avance de la computation como ayuda de calculo ha permitido el desarrollo de 
numerosos modelos matematicos para analizar el comportamiento dinamico de las estructuras 
reticulares. Sin embargo, la complejidad de estos sistemas hace necesario confirmar y mejorar 
estos modelos mediante tecnicas experimentales. La utilization de metodos para identificacion de 
sistemas permite que los resultados experimentales puedan usarse para definir con mejor 
precision sus caracteristicas dinamicas y calibrar sus modelos matematicos. 

Debido a la importancia de la infraestructura industrial para la economia de Mexico se ha 
hecho necesario el desarrollar algunas metodologias para la determination de las caracteristicas 
dinamicas de estructuras sometidas a carga de maquinaria. Entre estas metodologias destacan 
aquellas que se basan en el analisis de la serial de vibration en el dominio de la frecuencia. A 
partir de estas metodologias es posible evaluar el comportamiento dinamico de las estructuras 
industriales para contar con una estimation realista de los parameters dinamicos que 
caracterizan al sistema estructural y con base en ello determinar las fuerzas dinamicas 
generalizadas que actuan sobre la estructura en condiciones de operation. 

En este trabajo se desarrolla una metodologia para utilizar los resultados de mediciones 
de vibration efectuadas directamente en campo y a partir de estos datos hacer predicciones del 
comportamiento de estructuras sometidas a vibration forzada. 



2. IDENTIFICACION DEL SISTEMA ESTRUCTURAL 



2.1 Analisis de senales 

Mediante un analisis en el dominio de la frecuencia es posible determinar el sentido y la 
interdependencia estadistica de los movimientos de una estructura (Goldman, 1999). Para tal fin 
es necesario estudiar estos movimientos en varios puntos en las frecuencias de resonancia del 
sistema (Espinoza Barreras, 1999). El procedimiento experimental consiste en analizar dos 
senales simultaneamente, en donde uno de los acelerometros se coloca como referenda mientras 
que el otro se situa en puntos estrategicos sobre la estructura con el proposito de definir mejor el 
modo que se pretende medir. 

Como primer paso se evalua la Transformada Finita de Fourier de la serial (empleando el 
algoritmo conocido como Transformada Rapida de Fourier) mediante la siguiente expresion 
(Brighan, 1997): 






(i) 



En esta expresion x (f) es la serial registrada, t define la ventana de medicion, j es (-l) 1 ' 2 y f 
es la frecuencia de analisis. La ecuacion (l) define la forma en que una serial se distribuye con la 
frecuencia. Si al resultado de esta expresion se le calcula su modulo y se eleva al cuadrado se 
obtiene el espectro de potencia de x (t). 

Por otra parte, durante el proceso de medicion se pueden registrar varias senales 
simultaneamente, por lo que se puede calcular su funcion de correlation cruzada, esto es : 

R xyW = Y Jo x ( t )y( t + x ) dt 



( 2 ) 
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En la Relacion 2, y (t+x) representa otra de las senales registradas siendo x un tiempo de 
retraso. Esta funcion permite estimar en que orden una senal se correlaciona con otra, en que 
medida las senales fueron originadas por la misma fuente y con que retraso, y detectar la 
presencia de ruido o senales extranas. Cuando la Ecuacion 2 se aplica sobre la misma senal se 
obtiene la funcion de autocorrelacion^ 

RxxW^xttWt^t (3) 

Utilizando el concepto de Transformada Finita de Fourier, ecuacion (l), en las Ecuaciones 
2 y 3, se obtienem 



F x, ( f ) = J 0 ' R xy (t )e _J2rf> dt 


(4) 


F„(f) = jX(tU”'dt 


(5) 



Las Expresiones 4 y 5 son, respectivamente, las transformadas de las funciones de 
correlacion y autocorrelacion. Una aplicacion de la funcion F X y(f) consiste en el calculo del angulo 
de fase entre las componentes imaginaria y real de una senal. Este angulo se define de la 
siguiente manera : 



M f ) =tair ' 



' Q.,(f) s 



( 6 ) 



En donde Q xy (f) y C xy (f) representan respectivamente las partes imaginaria y real de F xy (f). 
La interpretacion fisica de este concepto se base en el hecho de que para cada frecuencia natural, 
la funcion de transferencia de fase correspondiente (relacion de fase entre ambas senales) es 
cercana a 0 o 180 grados. 

Las transformadas de las funciones de correlacion y autocorrelacion resultan de utilidad 
en el calculo de las funciones de coherencia y transference 



2 lf) MO 2 

T " J F„(f)F„(f) 



(7) 




( 8 ) 



En estas expresiones y xy (f) es la funcion de coherencia, H xy (f) es la funcion de transferencia 
y Fyy(f) denota la transformada de la funcion de autocorrelacion de la senal y(t). La principal 
utilidad de la funcion de coherencia radica en que proporciona una medida del grado de linealidad 
entre dos senales. Los valores de esta funcion van de 0 a 1. Cuando se alcanza el valor de 1 se dice 
que existe una perfecta relacion lineal entre ambas senales. 
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2.2 Determination de parametros estructurales 



La principal utilidad de los conceptos descritos, dentro de un contexto de analisis 
dinamico estructural, radica en la estimation de las frecuencias y modos naturales de vibracion 
que corresponden a los sistemas analizados. Para ello se procede de la siguiente manera : 

a) De los espectros de potencia, ( | x(f) | ) 2 , se determinan las frecuencias en las cuales se generan 
los picos mas grandes. Un pico esta asociado a un valor maximo de (|x(f)|) 2 . Estos picos 
representan estimaciones preliminares de las frecuencias naturales de vibracion. 

b) Se calculan las Transformadas Finitas de Fourier de las funciones de correlation entre dos 
senales registradas para una misma medicion. A estas transformadas se les obtiene su modulo y 
se eleva al cuadrado. Los resultados se dibujan en funcion de la frecuencia. En estas figuras se 
detectan los valores de frecuencia en que ocurren los picos o valores maximos de | Fxy(f) | 2 . 
Aquellos que coincidan con los detectados en el inciso anterior se podran considerar como 
asociados a las frecuencias naturales de vibracion. 

c) Se obtiene la funcion de Transferencia de Fase entre las senales utilizadas en el inciso (b) y se 
dibujan en funcion de la frecuencia. Para aquellos valores detectados como frecuencias naturales 
de vibracion se obtiene la magnitud de la Funcion de Transferencia de Fase. De acuerdo con ello 
se define el modo dominante de vibracion (lateral, vertical, torsional) asociado a la frecuencia pico. 

d) Se dibujan las Funciones de Coherencia entre las senales analizadas con el objeto de estimar la 
magnitud de los posibles ruidos registrados durante el proceso de medicion de las senales. 

e) Conocidas las frecuencias de interes, se procede a calcular la forma relativa de los modos de 
vibracion con ayuda de la expresion (Bendat and Piersol, 1961) 

( Pi(y,)= F y, yj ( f .) (9) 



Donde F v es la transformada de la funcion de autocorrelation para la frecuencia fi 

ynj 

asociada al punto de medicion j. Esta ultima expresion equivale a dibujar el espectro de potencia, 
para el modo de interes, en diferentes puntos de medicion de la estructura. (yj) y despues unir los 
picos que correspondan a las frecuencias naturales detectadas (fi). 



2.3 Matriz dinamica de flexibilidad 



En este apartado se desarrolla un modelo para determinar la matriz dinamica de 
flexibilidad de una estructura reticular sometida a vibracion forzada. Se acepta que el sistema 
estructural es lineal, elastico y no amortiguado. En consecuencia, por definition de modo de 
vibracion (Clough and Penzien, 1975) 

j[K]-p, 2 M]{U,M0} (10) 

En donde [K] es la matriz estatica de rigidez, pi es la frecuencia natural de vibracion (pi = 
2nfi), [M] es la matriz de masas, {Ui} es el vector modal asociado a pi. A partir de la Expresion 10 
puede demostrarse que (Ruiz Sibaja, 1991) 

{U, } T [K]{u , } = Pj 1 {u, } T [M]{Uj } 



( 11 ) 
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Por otra parte para una estructura sometida a vibracion forzada la ecuacion de 
movimiento es (Paz, 1997) 



Jk]-co 2 [m]}{x 0 }={F 0 } (12) 

En la Expresion 12 {xo} es el vector de amplitudes de desplazamiento, a> es la frecuencia de 
excitacion y {Fo} es el vector de amplitudes de la carga de excitacion. La amplitud de la respuesta 
dinamica se obtiene como (Clough and Penzien, 1975) : 

( x o} = £ a i( U i} (13) 

i=i 

En la cual ai es la amplitud modal y {Ui} el vector modal. A1 sustituir la Ecuacion 13 en la 
Ecuacion 12 y considerar la Ecuacion 11 se obtiene la expresion para la amplitud modal (Lopez 
Rios, 1990) 



„ _ i<J,) T {F„} 1 

' {U,} T [M]{U,}(p, J -a> 2 ) 



Al incorporar la Ecuacion 14 en la Ecuacion 13 resulta 



u 

1=1 



x 0 ( = 



En — n /TT lTr „v TT K u i fe 

tr( P i 2 -(0 2 j{U i } T [M]{U i } 



Se establecen las siguientes definiciones 



T> 2 2 

B i =Pi 



n, ={u,} T Wu, 



(14) 



(15) 



(16) 

(17) 



Se obtiene, a partir de la Ecuacion 15, la expresion para la matriz dinamica de flexibilidad 



fcj=£ 



1 

BA 



(u.Ku.r 



(18) 



En esta expresion es necesario conocer la frecuencia de excitacion, la frecuencia natural de 
vibracion y los vectores modales para determinar la matriz dinamica de flexibilidad de una 
estructura reticular, [Fdin]. 



2.3 Matriz dinamica de rigidez 

A partir de la Ecuacion 13 se puede demostrar que (Ruiz Sibaja, 1991) 

{u, } T [M]{x 0 } = a, {U, } t [mKu, } (19) 

En esta expresion puede incorporarse la expresion para la amplitud modal, Ecuacion 14, 
de lo cual resulta 
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! U ,) T[M] b„)=|^ (20) 

Que tambien puede escribirse 

B i {U, } T [M]{x 0 } = {U, } t {F 0 } (21) 

A1 premultiplicar por [M] y luego por {Ui} la Ecuacion 21 puede llevarse a la forma (Ruiz 
Sibaja, 1991) 



[m]{u , Xn , 1 ] B , {u , } t [m ]{x 0 } = {F„ } (22) 

E incorporando la Ecuacion 12 en la Ecuacion 22 y considerando la participacion de todos 
los modos de vibracion resulta 

Z[M]{U i }(fi i l'B 1 {U i } T [M]=[K-(0 ! M] (23) 

i=l 

Finalmente 

[K a „, ] = £ Mu i Xn, r B , {u, } T [m] (24) 

i=l 

Mediante la Ecuacion 24 se puede obtener la matriz dinamica de rigidez, [KdiJ, de una 
estructura reticular. A1 igual que en la Ecuacion 18 es necesario conocer la frecuencia de 
excitacion, la frecuencia natural de vibracion, la matriz de masas y los vectores modales para 
determinar [KdiJ. 



3. CALIBRATION DEL MODELO DESARROLLADO 



3.1 Modos de vibracion medidos 



Figura 1 Modelo estructural experimental. 
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A fin de realizar una calibration preliminar del modelo desarrollado previamente, se 
llevaron a cabo pruebas de vibracion forzada sobre una viga de acero doblemente articulada 
(figura l). Esta viga, ademas de su propio peso, carga cinco pesos concentrados tambien de acero. 
La vibracion forzada de la viga se consiguio mediante un motor electrico de 117 V, el cual hace 
girar un vastago con un peso de 20 grs. concentrado en el extremo del mismo. En la prueba inicial 
la frecuencia de excitacion se establecio en 12.20 Hz. 

Se utilizo un equipo de medicion a base de dos acelerometros piezoelectricos y un 
analizador de espectros con dos canales de entrada. Los acelerometros utilizados tienen una gama 
muy extensa de frecuencias (de 1 a 1000 Hz) y muy buena linealidad en todas ellas (Endevco, 
1971). Estos acelerometros se conectaron a cables blindados para enviar la serial de vibracion al 
analizador de espectros. El analizador que se uso cubre el intervalo de frecuencias de 0.02 Hz a 25 
Khz. (Hewlett Packard, 1979). A la serial registrada en cada uno de los acelerometros se le hizo un 
muestreo y se proceso mediante un filtro digital antes de almacenarla en la memoria del aparato. 
Estos datos almacenados en la memoria se analizaron mediante la Transformada Rapida de 
Fourier. 

El procedimiento de medicion consistio en analizar dos senales simultaneamente. Uno de 
los acelerometros se coloco fijo corno referenda (canal A), mientras que el otro se situo en puntos 
estrategicos sobre la viga (canal B), con el fin de definir la respuesta dinamica a rnedir. 



Figura 2 Posiciones de medicion de los acelerometros. 



De esta manera se seleccionaron cuatro posiciones diferentes de los acelerometros (figura 
2), para realizar series de mediciones sobre la viga de acero. En las figuras 3 y 4 se muestran, 
para la primera posicion de medicion, los espectros de transferencia de fase y de coherencia en el 
dominio de la frecuencia, en el intervalo de 0 a 100 Hz, despues de efectuar 32 promedios RMS. 
En todos los casos se tomo como referenda el acelerometro del canal A. 

Se observa que en los 6 Hz el valor de la funcion de transferencia de fase (Figura 3) es 
cercano a 0° (-4°), lo cual indica que las senales de los canales A y B estan en fase en esa 
frecuencia. Ademas el valor de la funcion de coherencia (Figura 4) cercano a 1.0 (0.94) senala que 
se trata de la misma respuesta dinamica (en este caso el primer rnodo de vibracion). 
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Figura 3 Funcion de transferencia de fase para la primera posicion de medicion. 

En 24.10 Hz el valor de la funcion de transferencia de fase (figura 3) es de 180°, indicando 
asi que la serial registrada por el acelerometro del canal B esta desfasada en 180° con respecto a 
la serial del canal A. El valor de 0.93 de la funcion de coherencia (Figura 4) confirma que se trata 
del segundo modo de vibracion. 



Figura 4 Funcion de coherencia para la primera posicion de medicion. 
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En la tabla 1 se resumen los valores obtenidos, a partir de los espectros en el dominio de 
la frecuencia, para los cinco primeros modos de vibracion que pudieron medirse en la viga. De 
acuerdo con los datos de la tabla 1 se pueden formar los vectores caracteristicos que corresponden 
a los modos de vibracion de la estructura, esto se consigue al considerar las amplitudes 
registradas en cada posicion y el comportamiento de las funciones de transferencia de fase y de 
coherencia. 



3.2 Matriz dinamica de flexibilidad medida 



En la figura 1 se resumen las caracteristicas geometricas y de material del modelo 
estructural experimental utilizado durante las pruebas de vibracion forzada. De acuerdo con los 
datos de esta figura se calculo una matriz teorica de masas concentradas (ecuacion 29). Mediante 
la aplicacion de la ecuacion 18 se obtuvo la matriz dinamica de flexibilidad medida. En el calculo 
se considero la participacion de los modos medidos de vibracion 1,2, 3, 4 y 5 y de la matriz de 
masas teorica, el resultado se muestra en la ecuacion 25 

Tabla 1 Resumen de resultados de mediciones de vibracion 
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1 

-0.264 




1,2, 3,4, 5 



SIMETRICA 



2 

-4.387 

-9.475 



La norma Frobenius de esta matriz es 



[fJ 



1,2, 3.4,5 || 2 



3 4 5 



-9.306 


-11.058 


-7.039" 


-15.383 


-17.424 


-11.059 


-16.731 


-15.698 


-9.658 




-9.357 


-4.503 






-0.444 


IV =5 


.530x10 2 





xlO 3 



(25) 



(26) 



3.2 Matriz dinamica de rigidez medida 



Mediante la aplicacion de la expresion 24 se obtiene, considerando la participation de los 
modos medidos de vibration 1, 2, 3, 4 y 5 y de la matriz de masas teorica 



1 


2 


3 


4 


5 




746.728 


-739.631 


303.722 


-100.454 


40.133 ' 


1 




1095.434 


-834.956 


333.080 


-83.104 


2 






1155.042 


-877.672 


353.329 


3 


SIMETRICA 






1159.323 


-783.481 


4 










722.416 


5 



Cuya norma Frobenius es 



[kJ, 



,2, 3, 4, 5 



£k ;j 2 =3302.933 



(27) 



(28) 



3.3 Matrices de flexibilidad y rigidez teoricas 



Se llevo a cabo el analisis modal del modelo estructural para encontrar las frecuencias 
naturales teoricas de vibration. De esta manera se obtuvieron los modos naturales de vibration 
cuyas frecuencias respectivas son pi = 6.172 Hz, p 2 = 24.698 Hz, p3 = 55.355 Hz, p4 = 95.664 Hz y 
ps = 136.678 Hz. A1 considerar la matriz de masas teorica (Ecuacion 29), los modos de vibration y 
sus frecuencias, se obtuvieron las matrices dinamicas de flexibilidad y rigidez teoricas, [F’din] y 
[K’din], 




1 2 

4.120 0 

0 4.257 

0 0 

0 0 

0 0 



3 4 

0 0 

0 0 

4.192 0 

0 4.179 

0 0 



0 

0 

0 

0 

4.193 



xlO 3 



(29) 
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A estas matrices se les calculo la norma Frobenius para comparar el resultado contra lo 
que se presenta en las Ecuaciones 26 y 28. En las Ecuaciones 29 y 30, se muestran los valores 
obtenidos para | | [F'duJ I 1 2, y I I [K'diJ I 1 2, respectivamente. 



feWIWIX =5.607 xlCT 2 (30) 



[KL], aj ,« || 2 = Ik, =3447.673 (31) 



A1 comparar los valores dados en las ecuaciones 26 y 28 (| | [Fdin] I I 2, y II [Kdin] I 1 2), 
contra lo que se presenta en las Ecuaciones 30 y 32 ( | | [F'diJ I 1 2, y I I [K'din] I 1 2) se observa una 
buena aproximacion entre lo que se obtiene mediante mediciones de vibracion forzada y lo 
pronosticado por la teoria. 



3.3 Comportamiento del modelo para diferentes frecuencias de excitation 



Para realizar una calibration mas amplia del modelo de analisis propuesto se estudio su 
respuesta ante diferentes frecuencias de excitation. Se aplicaron seis diferentes frecuencias de 
excitation: «i = 76.655 rad/seg, (02 = 219.597 rad/seg, «3 = 439.195 rad/seg, 004 = 658.792 rad/seg, 
0)5 = 878.389 rad/seg y«6= 1097.987 rad/seg. Para cada prueba de vibracion forzada se obtuvieron 
las matrices [Fdin] y [Kdin], es decir, se determinaron seis matrices [Fdin] y seis matrices [Kdin]. Se 
calculo la norma Frobenius para cada una de estas matrices y los valores obtenidos se escalaron 
respecto de los mayores valores calculados para | | [Fdin] I 1 2 y I I [Kdin] I 1 2 , respectivamente, de 
forma tal que se fijo una escala horizontal entre 0 y 1. 

En la figura 5 se muestra el comportamiento en funcion de la participation modal, de la 
norma | | [Fdin] I 1 2 escalada, para diferentes frecuencias de excitation. 

A partir de la figura 5 pueden considerarse los siguientes intervalos de comportamiento de 

I I [Fdin] | | 2 



0 < CO; 


<^(Pl +P 5 ) 


(32) 


(Pi +P 5 ) 


<®i <|(Pi +P 5 ) 


(33) 


(Pi +P 5 ) 


<®i <|(Pi +Ps) 


(34) 



Para el primer intervalo se observa que | | [Fdin] | 1 2 tiende a mantenerse constante a 
partir de la participation del 60% de la matriz modal (los tres primeros modos de vibracion). Esto 
significa que para este intervalo de excitation, la definition de [Fdin] no presenta variaciones 
apreciables a partir de la participation del tercer modo de vibracion. En el segundo intervalo, el 
comportamiento de | | [Fdin] I 1 2 tiende a mantenerse constante a partir de la participation del 
80% de la matriz modal, es decir, que se requiere de los cuatro primeros modos de vibracion para 
que la definition de [Fdin] no presente variaciones notorias. Conforme la frecuencia de excitation 
incursiona hacia valores superiores se hace necesaria la participation del 100% de la matriz 
modal, esto se advierte en el tercer intervalo, en el cual se necesitan de los cinco modos de 
vibracion para que [Fdin] alcance una buena definition. 
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Figura 5 Comportamiento de | | [Fdin] I I 2 para diferentes frecuencias de excitacion. 



Figura 6 Comportamiento de | | [Kdin] I I 2 para diferentes frecuencias de excitacion. 
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En la figura 6 se muestra el comportamiento en funcion de la participacion modal, de la 
norma | | [KdhJ I 1 2 escalada, para diferentes frecuencias de excitacion. Pueden considerarse los 
siguientes intervalos de comportamiento de | | [Kdin] I I 2 : 



3 

4 



0<co i <-(pi+p 5 ) 
(Pl+P 5 )<®i <|(Pi+Ps) 



(35) 

(36) 



Para el primer intervalo, se advierte que es necesario contar con la participacion del 100% 
de la matriz modal para alcanzar una buena definicion de [Kdin]. En el segundo intervalo, el 
comportamiento de | | [K di n] I 1 2 se mantiene casi constante a partir de la participacion del 80% 
de la matriz modal, por lo tanto, solo se necesitan cuatro modos de vibracion para que la 
definicion de [Kdin] no presente variaciones apreciables. 



4. CONCLUSIONES 



Se presento una metodologia para deducir las matrices dinamicas de flexibilidad y rigidez 
a partir de datos de mediciones de vibracion forzada. En esta etapa de la investigacion, es 
necesario proponer la matriz de masas del sistema por no haber sido posible deducirla 
experimentalmente. Aunque el estudio de un modelo de laboratorio con cinco masas concentradas 
no justifica una generalizacion, los resultados obtenidos con la metodologia propuesta muestran 
que es posible determinar, con una aproximacion adecuada, las matrices dinamicas de flexibilidad 
y rigidez, para estructuras sometidas a vibracion forzada. 

En las pruebas realizadas lograron identificarse cinco modos de vibracion en flexion del 
modelo estructural experimental. Las calibraciones del modelo matematico muestran que se 
obtiene una mejor definicion de [Fdin] cuando la frecuencia de excitacion se encuentra en la 
vecindad de las primeras frecuencias naturales de vibracion. Por el contrario, se observo una 
mejor aproximacion hacia [Kdin] cuando la frecuencia de excitacion incursiona en el intervalo de 
las frecuencias naturales de vibracion superiores. Los resultados obtenidos muestran, en general, 
que puede obtenerse una buena aproximacion de [Fdin] cuando se conoce el 60% de la matriz 
modal. En cambio, la definicion de [Kdin] requiere del conocimiento del 80% de la matriz modal. Se 
advirtio una gran similitud entre lo que se obtiene mediante mediciones de vibracion forzada y lo 
pronosticado por la teoria, probablemente porque se trata de una estructura sencilla. 
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RESUMEN 



Se hacen comentarios sobre los llamados criterios de diseno sismico por desempeno, 
indicando que partes de estos criterios, a juicio del autor, constituyen aportaciones novedosas, y 
que partes son ideas antiguas con nuevos nombres. Se hace ver que ciertas porciones de los 
criterios de diseno por desempeno han estado incluidas desde hace muchos anos en las normas 
convencionales. Posteriormente, se senalan algunos inconvenientes practicos que hay en la 
aplicacion de los criterios de diseno por desempeno, muchos de los cuales no son exclusivos de este 
nuevo enfoque, sino que son, desde hace tiempo, limitaciones del diseno sismico a secas. 
Finalmente, se indican algunos de los puntos en que las normas mexicanas estan actualmente 
anticuadas y se esbozan algunos caminos para mejorarlas, no necesariamente en el marco del 
diseno por desempeno sino en el marco del diseno racional. 



1. INTRODUCTION 



Las normas de diseno por sismo en Mexico tienen mas de 50 anos de antigiiedad. En 
varios momentos de su historia, las normas mexicanas han aportado ideas y practicas novedosas, 
algunas de las cuales han sido adoptadas despues en normas de otros paises (Fukuta, 1991): en 
1942, los factores de importancia; en 1957, la distribucion de fuerzas sismicas con la altura del 
edificio, el metodo dinamico de analisis, la limitation de desplazamiento lateral en las estructuras 
y la obligatoriedad de uso de coeficientes mayores en la zona de lagoi en 1976, factores novedosos 
de reduction por ductilidad. 

En la actualidad, sin embargo, muchas de las normas mexicanas son relativamente 
antiguas, por lo que no reflejan todos los avances recientes en la materia. 

En el presente trabajo se hace una breve reflexion sobre el estado actual de la 
normatividad en nuestro pais y sobre que modificaciones serian pertinentes a la luz de los 
llamados criterios de diseno por desempeno. 

Pero antes de discutir sobre las normas mexicanas, a continuation se hacen algunos 
comentarios sobre los propios criterios de desempeno, con el fin de hacer ver algunos 
inconvenientes practicos para su aplicacion y con la intention de senalar que algunos de estos 
criterios han estado presentes en las normas mexicanas desde hace mucho tiempo. 



2. QUE NO ES NUEVO EN EL DISENO POR DESEMPENO 

2.1. Tener objetivos de diseno 



Diversos autores han senalado (e.g, Hadjian, 2002) que el diseno siempre ha tenido 
objetivos de desempeno. Las estructuras siempre se han disenado para que “hagan algo”, es decir, 
para que tengan un cierto desempeno que alguien ha visualizado de antemano. Por esta razon, el 
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nombre “diseno por desempeno” es desafortunado y para muchos fines podria sustituirse 
simplemente por “diseno”. 

Los objetivos de diseno no siempre se han logrado; muchas estructuras han fallado a lo 
largo de la historia. Pero que no se hayan logrado siempre no se debe a que no hayan existido. Se 
debe a ignorancia, dolo, mala suerte, pero no a ausencia de objetivos de diseno. Esto nos conduce a 
una cuestion que el autor considera central siempre han existido objetivos de desempeno; lo que 
ha fallado es el como. La diferencia principal esta entre el diseno bueno y el diseno malo. 



2.2 El concepto de “nivel de desempeno” 



Una lectura somera (o profunda) de los textos seminales del diseno por desempeno nos 
hace ver que “nivel de desempeno” es un sinonimo perfecto de “estado limite”. No hare mas 
comentarios, salvo insistir en lo innecesario e inconveniente que es llamar a la misma cosa de dos 
maneras distintas. Propongo que mantengamos el acostumbrado “estado limite” y evitemos 
confusiones, aun a costa de perder cierta sensacion de modernidad y avance. 

2.3 Diseno multinivel o diseno para varios estados limite 



A continuacion se reproducen algunos articulos del Reglamento para las Construcciones 
del DF vigente en la actualidad. Estos dos articulos, en esencia, han estado presentes en los 
reglamentos del DF desde 1966. Las cursivas son del autor: 

> Art. 172 (Titulo Sexto del Reglamento de Construcciones del DF). “Este titulo contiene los 
requisitos que deben cumplirse... para lograr un nivel de seguridad aceptable i^que ocurra 
poco free uen temen te .9 contra fallas estructurales Q colapso ?) asi como un comportamiento 
estructural aceptable Q opei n ti vida d?) e n condiciones normales de operacion Q/recuente?) 

> Art. 182 “Toda estructura y cada una de sus partes deberan disenarse para cumplir con 
los requisitos basicos siguientes: 

> I Tener seguridad adecuada Q fallas poco frecuentes?) contra la aparicion de todo estado 
limite de falla posible (^colapso?) ante las combinaciones de acciones mas desfavorables 
que puedan presentarse durante su vida esperada, y 

> II No rebasar ningun estado limite de servicio Q o pei 'a ti vida d?) ante combinaciones de 
acciones que correspondan a condiciones normales de operacion Qcargas free uen tes?)’ 

Como se observa, la necesidad de que un diseno satisfaga simultaneamente varios estados 
limite asociados a diferentes frecuencias de ocurrencia ha sido desde hace muchos anos una 
aspiracion explicita de los criterios de diseno. Hadjian (2001) hace ver algo parecido, aunque mas 
dramatico : la aspiracion explicita del documento Vision 2000 (Hamburger et al, 1995), 
considerado la piedra de toque del diseno por desempeno, es practicamente igual a la expresada 
por la SEAOC en 1967. 

Nuevamente, el problema ha sido el como, no la ausencia del deseo de disenar para 
satisfacer simultaneamente varios estados limite. 



2.4 Disenos mas conservadores para estructuras importantes 

Desde hace mucho tiempo las estructuras mas importantes han sido disenadas para 
satisfacer estados limites mas estrictos. Los llamados factores de importancia (los que se aplican a 
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estructuras del Grupo A, del Grupo B o del Grupo C) estan presentes en los reglamentos 
mexicanos por lo menos desde 1966. 

No esta claro que el uso de factores de importancia conduzca siempre a disenos adecuados; el uso 
de un factor de importancia independiente de la sismicidad de la zona parece tener poco sentido. 
Es, nuevamente, un problema de ignorancia sobre el fenomeno sismico y dificultad para hacer 
muchos calculos, pero no una falla esencial de los criterios de diseno. 



3. CONTRIBUCIONES REALES DEL LLAMADO DISENO POR 
DESEMPENO 

3.1 El deseo de usar mejores indicadores de desempeno 



Ha sido muy importante para el progreso del diseno por sismo la insistencia, 
relativamente reciente, sobre usar mejores indicadores del desempeno estructural. En particular, 
la difusion del uso de los desplazamientos como medidas de la respuesta ha sido muy benefica. 
Esto esta relacionado con un deseo, correcto desde mi punto de vista, de saber en que condiciones 
va a quedar la estructura despues de sismos que ocurren con frecuencias dadas. 



3.2 El recordatorio de que satisfacer el estado limite de colapso no lleva a satisfacer 
todos los demas. 



Aunque el diseno multinivel —satisfacer varios estados limite simultaneamente - ha 
existido desde hace mucho tiempo, tambien ha existido la percepcion de que si una estructura 
puede lo mas, tambien puede lo menos. Es facil demostrar que esto, en muchos casos, no es cierto : 
una estructura puede tener un comportamiento aceptable ante el sismo de colapso y un 
desempeno pobre ante sismos frecuentes de servicio; las estructuras flexibles de la ciudad de 
Mexico son un ejemplo notable de esta situacion. Por ello, la insistencia en disenar para, 
realmente, satisfacer varios estados limite es muy pertinente. 



3.3 El senalamiento de las limitaciones del uso de factores de importancia 



Se ha hecho ver que solo con incrementar las cargas de diseno para estructuras mas 
importantes no necesariamente se consigue un comportamiento satisfactorio. En efecto, para 
estructuras importantes no basta con incrementar las cargas de diseno : es necesario especificar 
detalladamente los estados limite deseados asociados al desempeno deseado. El diseno por 
desempeno ofrece un excelente marco para resolver esta limitacion de los factores de importancia. 



3.4 Tasas de excedencia versus tasas de falla 



En el estado actual del diseno por desempeno, se especifica la frecuencia de ocurrencia de 
las demandas asociadas a los estados limites considerados. Como se vera a continuation, esto no 
garantiza que se llegue a estructuras de confiabilidad uniforme, es decir, estructuras cuyas fallas 
acontezcan con la misma frecuencia promedio. Esto se debe a que, en general, las fallas pueden 
ocurrir aun con temblores de demanda nominal menor a la especificada en el codigo y, 
eventualmente, no ocurrir con temblores con demanda nominal superior. 
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Para ilustrar lo anterior proponemos el siguiente ejemplo. Supongamos que se tiene un 
tipo de estructuras que se disenaran para el estado limite a=c, donde a es la demanda y c la 
capacidad. Sea v A (a) la tasa de excedencia de la demanda, es decir, el numero medio de veces en 
que, por unidad de tiernpo, la demanda excedera un valor particular, a. Supongamos que Va(^) 
tiene la siguiente forma funcional (Esteva, 1970) : 

v A (a) = Ka ( 1 ) 

donde los parametros K y r (K>0, i< 0) dependen de la sismicidad regional. Se desea que este 
estado limite no sea excedido mas de, digamos, v 0 veces por ano, lo cual se supone puede lograrse 
escogiendo como capacidad de diseno una cuyo valor numerico es igual al de la demanda asociada 
a una tasa de excedencia igual a v 0 (ver figura l). Llamaremos co a esa capacidad de diseno. 

Si la capacidad real de la estructura fuera determinista e igual a co, entonces, en efecto, la 
estructura fallaria cuando a>co, lo cual aconteceria exactamente v„ veces por ano. En otras 
palabras, la tasa de falla de la estructura seria precisamente la tasa de excedencia de la demanda 
asociada al estado limite. 




Sin embargo, por diversas razones, la capacidad estructural no es determinista, sino que 
existen incertidumbres en su calculo. Por sencillez, modelaremos esta incertidumbre de la 
siguiente manera : supondremos que la capacidad real, c, es una variable aleatoria con 
distribucion lognormal, con mediana igual a la capacidad nominal de diseno, co, multiplicada por 
un factor de sobrerresistencia, F, y una cierta desviacion estandar del logaritmo natural de c igual 
a <J. La aplicacion del teorema de la probabilidad total nos lleva a que la tasa de falla cuando se 
adopta para diseno una capacidad nominal Co, Vc(co), esta dada por: 



v c( c o ) = -{ p r( c < a\c 0 ) 
0 



dv A ( a) 



da 



da 



( 2 ) 



donde Pr(c<a Ico) es la probabilidad de que la capacidad (real) sea menor que la demanda dado 
que se adopto para diseno una capacidad nominal co. Con la hipotesis simplificatorias 
mencionadas, esta probabilidad se determina de la siguiente manera : 
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Pr( c < a\c 0 ) = J 



1 ~(ln—) 2 

e 2a Fc ° dc 

lucre 



(3) 



Sustituyendo las ecuaciones 3 y 1 en la 2, y resolviendo las integrales correspondientes, 
obtenemos que: 



V , 



1 2 2 J 

dc 0 ) = K(Fc 0 }e r " = v(Fc n )e‘ 



= y 0 F' 



(4) 



Se observa, como ya hablan senalado otros autores (e.g., Esteva y Ruiz, 1989) resolviendo 
casos mas realistas, que las tasas de falla pueden ser significativamente diferentes de las tasas de 
excedencia prefijadas, excepto en el caso determinista en que F=1 y <j=0, con lo cual Vc(co)= Vo- 

La ecuacion 4 ilustra varios hechos importantes. 

1) Las tasas de falla, las que verdaderamente importan, son diferentes de las tasas de 
excedencia de la demanda, que son las que actualmente se estipulan en las propuestas 
reglamentarias. 

2) La tasa de falla de una estructura en particular depende, aun en nuestro ejemplo 
simplificado, de su sismicidad regional. Dos estructuras iguales localizadas en diferentes 
partes del pais (r diferente, en nuestro ejemplo simplificado) tendran confiabilidades 
diferentes aun si ambas se disenan para una demanda asociada al mismo periodo de 
retorno. 

3) La tasa de falla depende de la incertidumbre que exista en el calculo de la capacidad 
estructural o, de manera mas general, de su vulnerabilidad. Dos estructuras construidas 
en el mismo sitio tendran tasas de falla diferentes si sus caracteristicas de vulnerabilidad 
(ejemplificadas con F y a en el ejemplo simplificado) son diferentes, aun cuando hayan 
sido disenadas ambas para la demanda asociada al mismo periodo de retorno. 

Esta situation no debera ser ignorada en la normatividad futura. La solution a secas es 
simple^ que los reglamentos estipulen parejas estado limite-tasa de falla y no parejas estado 
limite-tasa de excedencia de la demanda, como se hace en la actualidad. La solution practica, sin 
embargo, es sumamente compleja. 

En efecto, para determinar aceptablemente la tasa de falla de la estructura no bastaria 
con calcular su probabilidad de falla ante el sismo con tasa de excedencia Vo (la tasa de excedencia 
de la demanda de diseno) sino que, de acuerdo con la ecuacion 2, seria necesario determinar las 
probabilidades de falla condicionadas a la ocurrencia de varios niveles de demanda y, 
posteriormente, sumar ponderadamente dichas probabilidades, usando como factores de peso las 
frecuencias de ocurrencia de cada uno de los valores de demanda elegidos (En el fondo, la 
ecuacion 2 no es mas que una suma pondera da de probabilidades de falla). Esto es impracticable 
con las herramientas y los costos de hoy. 

Sin embargo, el resultado expresado en la ecuacion 4, aun tratandose de un caso 
simplificado, sugiere posibles caminos para acercarse al objetivo de diseno para confiabilidades 
predeterminadas y constantes : 

l) Utilizar factores de sobrerresistencia o factores de seguridad diferentes (el factor F en 
nuestro ejemplo) dependiendo del tipo de estructura y la sismicidad regional. 



2) Utilizar periodos de retorno diferentes para la demanda de diseno (diferentes Vo) 
dependiendo del tipo y localization de la estructura. 
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Esto, desafortunadamente, eliminaria la aparente sencillez y transparencia de los planteamientos 
conceptuales del llamado diseno por desempeno. El tema amerita reflexion. 



3.5 ^Tasas de falla uniformes? 



Hemos visto que, aun si se estipulan periodos de retorno uniformes a lo largo de un pais 
para las demandas de diseno, no necesariamente se consiguen estructuras con confiabilidad 
uniforme. Pero cabe la siguiente pregunta : nadie duda que es util disenar estructuras con 
confiabilidad conocida, pero <^es deseable tener estructuras de diferentes tipos con confiabilidad 
uniforme a lo largo de regiones extensas?. Para responder a esta pregunta, que a mi juicio sigue 
abierta, se han esgrimido argumentos de muy diferente indole : equidad, justicia, optimacion de 
costos, y otros. 

En opinion del autor, que una confiabilidad sea adecuada para una estructura en una 
parte del pais no implica que tambien lo sea para otra estructura en otro lado. 



Los analisis de optimacion de costos indicaron desde muy temprano en el desarrollo de la 
ingenieria sismica (Esteva, 1970) que la estrategia de confiabilidad uniforme conduce casi 
inevitablemente a disenos sub-optimos. La razon es clara : el costo de dar resistencia lateral a una 
estructura crece mas rapidamente que, digamos, el coeficiente adoptado para diseno. Esto 
significa que en las zonas de baja sismicidad, con niveles de diseno relativamente bajos, es mas 
barato “sobre- disenar” una estructura que en una zona de alta sismicidad, por lo cual uno se 
puede dar el lujo de ser mucho mas conservador. 



Ilustraremos esto siguiendo el planteamiento de Esteva (1970), hecho hace mas de 30 
anos. Supongamos que el costo de construccion, Cc, varia con la resistencia nominal ante cargas 
laterales, co, de la siguiente manera : 



Cr( C 0 ) = 



K, 



k 1 + k 2C ; 



si c 0 < c K 
si c 0 > c K 



(5) 



donde K 2 es una constante de normalizacion y Ki es el costo fijo de construccion que se tendria 
aun sin diseno ante cargas laterales. El exponente a es mayor que 1, indicando que el costo de 
suministrar resistencia lateral crece mas rapidamente que el valor de esta resistencia lateral. 

Por otra parte, puede demostrase que el costo esperado — actualizado a valor presente- de 
las perdidas por sismo, Cp, suponiendo que estas ocurren de acuerdo con un proceso de Poisson y 
que se utiliza para diseno una resistencia co, esta dada por: 



C P ( c 0 ) = D 



K(Cp) 

Y 



( 6 ) 



donde D es el costo de la falla. El costo total de la decision de adoptar co como valor de diseno, Ct, 
es entonces la suma del costo de construccion mas el costo esperado de las perdidas por sismo. La 
mejor decision es la que minimiza el costo total, es decir, aquella para la cual 



dC T ( c 0 ) _ q 
dc 0 



(7) 



El coeficiente optimo de diseno, cop, resulta ser 
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/ KrD ' 



c op ~ 



, 2 2 

r+a _ iL ++ 
0 2 r+a 



KjOF' 

que resulta estar asociado a una tasa de falla, v OP , dada por 



( 8 ) 



V OP 



r K a/r ayK 1 V+« 

rDF a 



1 r cj a 

2 r+a 



(9) 



Como puede observarse, la tasa optima de falla depende de todos los parametros que 
hemos hecho intervene sismicidad, vulnerabilidad estructural, costo de construccion y costo de la 
falla. La tasa optima de diseno (el periodo de retorno de la falla) depende de la localizacion, el tipo 
estructural, el costo de construccion y el costo de la falla. No hay, pues, tasas optimas de falla 
generales. 

Analicemos un caso simple^ dos estructuras iguales (que tienen, por tanto, iguales F, a, K h 
a y D) seran construidas en dos sitios diferentes, tales que en el sitio A la sismicidad es 
exactamente el doble que en el sitio B. (igual rpero diferente A). Supongamos que en el sitio B -el 
de menor sismicidad- se ha calculado que 0.001/ano es la tasa optima de diseno (periodo de 
retorno optimo para la falla = 1000 anos). En ese caso, la tasa optima de diseno en el sitio A, 
aplicando la ecuacion 9, resultaria sen 



V OPA V OPB 



K, 



r+a f l\ 

= 0.001 - 

\2 J 



( 10 ) 



Si 05=1.5 y r =-2.5 (valores tipicos de estos parametros), tendriamos que la tasa optima de 
diseno en el sitio A -el de mayor sismicidad- seria de 0.0028/ano, asociada a un periodo de retorno 
optimo para la falla de solo 353 anos. En otras palabras, para el sitio de mayor sismicidad no es 
optimo ser tan conservadores como para el sitio de menor sismicidad. 

Desde este punto de vista, no parece razonable, por ejemplo, que si se juzga adecuado 
disenar una estructura para 150 anos de periodo de retorno en Acapulco sea razonable tambien 
hacerlo para el mismo periodo de retorno en Aguascalientes o Ciudad Juarez. El analisis de costos 
— y en parte el sentido cornun- indican que en Aguascalientes o Ciudad Juarez deberia aspirarse a 
mayores confiabilidades puesto que resulta barato. 

Quiza la estrategia mas adecuada sea imponer una confiabilidad minima, socialmente 
aceptable, para las zonas de alta sismicidad y aspirar a confiabilidades progresivamente mayores 
conformes la amenaza sismica vaya disminuyendo’ . 



4. ESTADO ACTUAL DE LA NORMATIVIDAD EN MEXICO Y MEJORAS 
FUTURAS 



Como se sabe, la normatividad de la construccion en Mexico es de competencia municipal. 
Por ello, en principio, cada uno de los cerca de 2500 municipios de Mexico podria tener sus propias 
normas de diseno y construccion 1 . A esto hay que anadirle las normas que algunos sectores 
industriales desarrollan para sus propias obras (CFE, por ejemplo). 



1 Esto no ha ocurrido asi en la practica. Ha sido mas comiin que se desarrollen normas de alcance estatal que 
posteriormente son adoptadas por todos los municipios del estado. 
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En el pasado se han hecho esfuerzos por desarrollar normas tipo que, adaptadas a la 
sismicidad y condiciones del suelo locales, constituyan la base de las normas estatales o 
municipales. Todos estos esfuerzos han fracasado. 

Existian hasta 2000 por lo menos 40 reglamentos diferentes en Mexico (Sanchez, 2000). 
En terminos generales, todas las normas mexicanas de diseno por sismo actualmente vigentes 
estan basadas en dos documentos principales : la Normas Tecnicas Complementarias para Diseno 
por Sismo del DF, y el Capitulo de Diseno por Sismo del Manual de Obras Civiles de la Comision 
Federal de Electricidad. Las normas del DF datan de 1985 (con actualizaciones en 1987 y 1992) 
mientras que las normas de CFE fueron emitidas en 1993. 

Para fines practicos, las normas de CFE fungen como una especie de reglamento nacional, 
puesto que, sin ser obligatorias, son usadas en muchas obras construidas en sitios para los que no 
se dispone de reglamento propio y, en algunos casos, aun en sitios donde si hay reglamento local. 
Puesto que no se preve un cambio drastico de esta situacion a mediano plazo, hay que resaltar la 
gran importancia de que las normas de CFE se mantengan al dia. 

En los siguientes parrafos haremos mencion de algunos inconvenientes que presentan la 
mayor parte de las normas mexicanas actualmente vigentes e iremos senalando algunas de las 
mejoras que deberian hacerse. 



5 NIVELES DE PELIGRO SISMICO ASOCIADOS A LOS ESPECTROS DE 
DISENO 



Salvo excepciones, los reglamentos mexicanos no especifican los niveles de peligro 
implicitos en los espectros de diseno. En otras palabras, rara vez se sabe a que periodo de retorno 
estan asociadas las ordenadas espectrales de diseno. El que los periodos de retorno no se 
especifiquen tiene el inconveniente de que el responsable del diseno o el dueno de la obra no 
saben, ni siquiera de manera aproximada, que tan frecuentemente se excederan los estados 
limite. Por esta razon, el disenador y el dueno tienen solo una vaga idea del desempeno futuro de 
la estructura. Pero la ausencia explicita de referencias a los periodos de retorno en los 
reglamentos quiza esconde un problema mayor: no tenemos criterios robustos y de aceptacion 
general para fijar esos periodos de retorno. 

Como se ha senalado en la introduccion, el problema de asociar periodos de retorno a los 
estados limite no es trivial. Da la impresion de que no es razonable usar el mismo periodo de 
retorno para todo el pais ni para todos los tipos estructurales. Los acuerdos sobre este asunto que 
se han alcanzado en otros paises nos ayudan poco : que un periodo de retorno de 500 anos sea 
considerado adecuado para cierto estado limite en una region de los Estados Unidos y los tipos 
estructurales que ahi se construyen, no implica que tambien lo sea para todo Mexico, con sus 
propios tipos estructurales y un nivel de desarrollo econo mico muy diferente. 

A pesar de su complejidad, es necesario avanzar rapidamente en la solucion de este 
problema. Si no, veremos en el futuro la proliferacion de normas o estudios especificos de sitio con 
espectros asociados a muy variados niveles de riesgo tolerable. Esta situacion traera confusion y, 
muy probablemente, uso sub-optimo del dinero. 



6. INCLUSION DE LOS EFECTOS DE SITIO 



Como se sabe, los efectos de amplification de las ondas sismicas producidos por la geologia 
local, tambien conocidos como efectos de sitio, son probablemente el principal factor de riesgo 
sismico para las construcciones. A pesar de ello, las normas mexicanas (con la exception de las del 
Distrito Federal) dan al asunto relativamente poca importancia. 
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A medida que la instrumentacion slsmica ha ido progresando, con lo que ha aumentado 
considerablemente el numero de acelerogramas registrados, ha ido quedando mas claro que 
clasificar los suelos en “firmes, intermedios y blandos” es decididamente inadecuado. 

En nuestro pals, el caso mejor estudiado es la Ciudad de Mexico. Quiza se trata de un caso 
extremo, pero las variaciones que hay en las caracterlsticas del movimiento del suelo en unos 
cuantos cientos de metros es notable. Clasificar toda la zona del lago como “suelo blando” no 
ayuda a lograr mejores disenos para resistir temblores. 

Ademas de contar con muchos mas acelerogramas registrados, en la actualidad se dispone 
de mejores herramientas teoricas y metodos de canrpo relativamente baratos que permitirlan 
mejorar nuestra capacidad para estimar movimientos fuertes en un mayor numero de sitios y con 
esto poder incluir en los reglamentos los efectos de sitio de mejor manera. 

En este sentido, podrla avanzarse mucho llevando a cabo estudios de microzonificacion 
slsmica de las ciudades mas expuestas en Mexico. En su mayorla, los estudios realizados hasta la 
fecha en nuestro pals son inadecuados, porque se limitan a clasificar los suelos segun los antiguos 
criterios (firme, intermedio y blando) y practicamente no dan indicaciones cuantitativas utiles 
para la estimacion de movimiento fuerte. Algunos palses, como Colombia, han vuelto obligatoria 
la elaboracion de estudios de microzonificacion slsmica para ciudades con ciertas caracterlsticas; 
quiza deberlamos seguir el ejemplo. 

A falta de buenos mapas de microzonificacion, las normas futuras podrlan reglamentar el 
uso de espectros especlficos de sitio y exigir su ejecucion para estructuras de cierta importancia. 



7 FORMA Y TAMANO DE LOS ESPECTROS DE DISENO 



En general, los espectros de diseno en las normas mexicanas no reflejan el verdadero 
tamano ni la forma real de los espectros de los temblores para los que se quiere disenar. Aun 
cuando las fuerzas y desplazamientos de diseno sean correctamente calculados, utilizando 
factores correctivos adecuados, el tener espectros de diseno que no se parezcan suficientemente a 
los de respuesta genera problemas de falta de claridad y perdida de perspectiva durante el 
proceso mismo de diseno. Algunos ejemplos : un espectro de diseno para suelo blando que tenga 
una meseta demasiado ancha, no transmite correctamente la idea de que hay que evitar que el 
periodo fundamental de la estructura este cerca del periodo predominante del suelo; o un espectro 
de aceleraciones que no tenga una caida en periodo largo compatible con las observaciones y la 
dinamica estructural, llevara a malas estimaciones de los desplazamientos. 

Un colofon de esta consideracion es que las amplificaciones (con respecto al terreno firme) 
observadas por efectos de sitio en una ciudad deberian estar reflejadas en los espectros de diseno. 
En otras palabras, si se observa que, para cierta ordenada espectral las aceleraciones en sitios 
blandos son mucho mayores que en terreno firme, las ordenadas de diseno deberan ser tambien 
mucho mayores. 

Una de las razones que han hecho que en la practica mexicana los tamanos de los 
espectros de diseno no reflejen los tamanos de los espectros reales es que en las normas 
generalmente se han incluido reducciones a las fuerzas de diseno de las que el disenador no esta 
suficientemente advertido. Estas reducciones impbcitas, que han tratado de tomar en cuenta 
efectos como la sobrerresistencia estructural, probablemente no llevan a disenos inseguros, pero si 
producen cierta falta de claridad que no es deseable. Por esta razon, en las normas futuras debera 
procurarse que todas las reducciones a los espectros elasticos de diseno se hagan explicitas, aun 
cuando los factores de reduccion usados no puedan justificarse completamente con bases teoricas 
o empiricas. 
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8. DEFINICION MAS CLARA DE LOS ESTADOS LIMITE 



Todos los reglamentos mexicanos vigentes estipulan que las estructuras deben disenarse 
para satisfacer dos estados limite, generalmente llamados de colapso y de servicio. El estado 
limite de colapso, con las salvedades que se han senalado anteriormente (periodo de retorno 
desconocido, formas espectrales inadecuadas, reducciones por sobrerresistencia implicitas), esta 
bien definido : la estructura debe tener propiedades de rigidez y resistencia tales que la demanda 
global de ductilidad sea inferior a la capacidad ductil de la estructura cuando ocurre un temblor 
caracterizado por el espectro de diseno. El estado limite de servicio, sin embargo, no esta 
claramente especificado. La revision requerida por los reglamentos mexicanos consiste en calcular 
los desplazamientos de la estructura ante las cargas correspondientes al estado limite de colapso 
y despues compararlos con unos desplazamientos permisibles, independientes del sistema 
estructura, que no son los que provocarian la falla de la estructura ni aquellos para los cuales se 
iniciaria el agrietamiento. Esta revision esta disenada para que los calculos sean simples y para 
que solo tenga que hacerse un analisis estructural; los desplazamientos permisibles estan 
calibrados para que el diseno lleve a rigideces adecuadas ante sismos de baja intensidad, cuyo 
tamano o periodo de retorno tampoco se especifica. 

El empleo de estos procedimientos de diseno lleva, generalmente, a estructuras 
razonablemente seguras. Sin embargo, es muy poco transparente y hace muy dificil que el 
disenador tenga una idea clara del desempeno de la estructura durante los sismos de diseno. 

Es deseable, entonces, que los estados limite se especifiquen con mayor claridad aun 
cuando esto suponga incrementos moderados en el trabajo de analisis estructural. En la nueva 
propuesta de Normas Tecnicas Complementarias para Diseno por Sismo del DF se ha avanzado 
en este aspecto. En dichas normas, se hace explicita la existencia de dos estados limite de 
desplazamiento, con distorsiones de entrepiso permisibles que reflejan mejor el desempeno 
estructural que se quiere obtener en los dos estados limites considerados. Es deseable que estas 
modificaciones se extiendan al resto de los reglamentos mexicanos. 



9. FACTORES DE REDUCCION POR RESISTENCIA 



Los criterios contemporaneos de diseno admiten que la estructura tenga incursiones en el 
rango no lineal ante el sismo que caracteriza el estado limite de colapso. Esto permite limitar las 
demandas de fuerza en los elementos estructurales (para, por tanto, utilizar resistencias de 
diseno menores), a costa de que se presenten demandas de ductilidad limitadas y cierto nivel de 
danos provocados por la fluencia de algunas secciones de la estructura. 

Para modelar el comportamiento no lineal, practicamente todos los reglamentos del 
mundo estan basados en analisis de un sistema de un grado de libertad con comportamiento 
elastoplastico. Es con este modelo que se determina la resistencia necesaria para limitar las 
demandas de ductilidad a un valor especificado, que llamaremos Q. Es usual expresar la 
resistencia necesaria para lograr una demanda de ductilidad dada, C(T, Q), como una fraccion de 
la resistencia necesaria para tener una demanda de ductilidad unitaria (comportamiento 
elastico), C(T,1 '), notese que ambas resistencias dependen del periodo estructural, T. Llamaremos 
Q ’ a este cociente: 



Q\t,q) 



C (T ,1) 
C(T,Q) 



( 11 ) 



La forma de la funcion Q ha sido ampliamente estudiada en los ultimos anos (Krawinkler 
et al., 1992; Miranda, 1993; Miranda y Bertero, 1994). En particular, Ordaz y Perez Rocha (1998) 
observaron que, bajo circunstancias muy generales, Q’ depende del cociente entre el 
desplazamiento espectral y el desplazamiento maximo del suelo. Esto implica que, como 
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observaron Meli y Avila (1988) al analizar registros del sismo de Michoacan obtenidos en la zona 
del lago del Valle de Mexico, los factores Q dependen fuertemente de las caracteristicas del suelo 
en que se desplanta la estructura. En vista de esto, es claro que las reglas que actualmente se 
utilizan en los reglamentos mexicanos, basadas en las derivadas por Newmark y Hall, son 
inadecuadas para suelos blandos por lo que es necesario sustituirlas por otras mas realistas. 



10. FACTORES DE IMPORTANCIA PARA ESTRUCTURAS DEL GRUPO A 



Los reglamentos mexicanos vigentes clasifican a las estructuras en grupos dependiendo de 
su importancia. La mayor parte de los requisites de diseno son los mismos para todos los grupos, 
excepto que se multiplican por un factor de 1.5 las fuerzas de diseno cuando se trata de 
estructuras importantes, generalmente llamadas del grupo A. 

Esto tiene al menos dos inconvenientes. l) El incremento en confiabilidad estructural que 
se logra al disenar para fuerzas mayores depende de la sismicidad regional. En otras palabras, 
dependiendo del regimen sismico del sitio, un incremento de 50% en las fuerzas de diseno puede 
significar un incremento de confiabilidad mucho mayor que para la misma estructura localizada 
en otra parte. En general, seria preferible estipular una confiabilidad objetivo para estructuras 
del grupo A y despues ver que fuerzas de diseno deberian usarse para conseguir esa confiabilidad. 
Esto, como se ha hecho ver anteriormente, conduciria a factores de importancia variables por tipo 
estructural, tipo de suelo y localizacion geografica de la estructura. 2) No esta claro que en todos 
los casos se consiga un desempeno satisfactorio de una estructura del grupo A mediante el solo 
incremento de las fuerzas de diseno. Seria deseable, ademas de especificar una confiabilidad 
superior, especificar claramente objetivos especiales de desempeno para estructuras importantes, 
tal como lo senala el documento Vision 2000. 



11. CONCLUSIONES 



La normatividad en materia de diseno por sismo en Mexico presenta dos problemas 
principales. Por una parte, la ausencia de uniformidad en los criterios adoptados en las distintas 
normas puede conducir paulatinamente a estructuras con confiabilidades muy diferentes 
dependiendo de las normas que se hayan empleado en su diseno. Para combatir este problema es 
necesario que los diferentes responsables de la elaboracion de normas (municipios, gobiernos 
estatales, empresas paraestatales, grupos universitarios) acuerden criterios generales de diseno 
comunes, alrededor de los cuales podrian elaborarse las normas locales. Por otra parte, ciertos 
criterios basicos de diseno necesitan cambios, drasticos en algunos casos, para poner las normas 
al dia e incluir los avances que recientemente se han tenido en diferentes aspectos del diseno por 
sismo. Se han mencionado en este trabajo los cambios que, a juicio del autor, son prioritarios. En 
algunos casos se han senalado posibles soluciones; en otros, sin embargo, solo se ha planteado el 
problema. 
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RESUMEN 



Despues de un terremoto fuerte, los danos en el area afectada son tan generalizados que 
no es posible que todas las evaluaciones de dano de los edificios sean hechas por ingenieros 
expertos. Por esto profesionales voluntaries con poca o ninguna experiencia en evaluacion del 
dano sismico deben evaluar gran parte de los edificios afectados. Los evaluadores inexpertos 
presentan la tendencia a agravar o subestimar el nivel real del dano. En estos casos la 
informacion es subjetiva e incompleta y se utilizan calificaciones linguisticas. Por lo anterior, se 
desarrollo un sistema experto para la evaluacion post-sismica del dano en edificios, utilizando 
tecnicas como las redes neuronales artificiales y la teoria de los conjuntos difusos. Esta 
herramienta facilita la participacion de inexpertos en el proceso de evaluacion. Este modelo de 
Inteligencia Artificial tiene en cuenta los posibles niveles de danos en elementos estructurales y 
no estructurales, los efectos del sismo en el suelo, asi como las condiciones preexistentes que 
pueden hacer mas vulnerable el edificio, como la calidad de los materiales de construccion, 
irregularidades en planta y altura y una mala configuracion estructural. El sistema decide sobre 
la habitabilidad y reparabilidad del edificio aplicando una base de reglas difusas a la informacion 
obtenida del edificio. 

1 INTRODUCTION 



En el caso de un sismo fuerte, debido a la magnitud del desastre, es necesario que el 
proceso de evaluacion del dano en los edificios sea hecho por profesionales con experiencia en esta 
tarea y que esten familiarizadas con la ingenieria estructural. Desgraciadamente el personal que 
cumple con estas caracteristicas es muy poco, por lo cual en casos de dano generalizado en un 
centro urbano, es necesario involucrar ingenieros y arquitectos inexpertos. Como consecuencia, es 
comun la sobreestimacion o sub estimation del dano, llevando a la ocupacion peligrosa o la 
demolition innecesaria de las edificaciones. Motivado en esta problematica, este trabajo propone 
un sistema experto para apoyar el proceso de evaluacion de danos, usando herramientas como las 
redes neuronales artificiales y la teoria de conjuntos difusos. 



2 EVALUACION DEL DANO DESPUES DE UN TERREMOTO 



Como resultado de los sismos ocurridos en diferentes paises ubicados en zonas sismicas, 
ha sido necesario que estos desarrollen manuales para la evaluacion del dano en edificios, con el 
proposito de decidir si los edificios pueden ser usados o no lo mas pronto posible. Despues de un 
sismo moderado o severo, es importante identificar las edificaciones que han sufrido serios danos 
y que pueden representar peligro para sus ocupantes o la comunidad en general, asi como 
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aquellas que pueden ser utilizadas como refugios temporales para personas que han sido 
evacuadas. Algunos paises han desarrollado manuales y procedimientos para evaluar el dano en 
edificios, algunos son: Japon, Macedonia (antes Yugoslavia), Italia, Grecia, Estados Unidos, 
Mexico y Colombia, entre otros. Las evaluaciones de dano no son solo utiles para la toma de 
decisiones sobre las estructuras, sino que tambien sirven para hacer mejoras a los codigos sismo- 
resistentes, identificando los tipos de falla en los sistemas estructurales y las razones por las que 
se presentaron. Ademas, con la informacion a cerca del dano en los edificios es posible conformar 
funciones de vulnerabilidad empirica. 



2.1 Problemas con la evaluation del dano 



Cuando el dano en el area afectada es generalizado, los expertos locales en ingenieria 
estructural son insuficientes para hacer la totalidad de las evaluaciones, asi que profesionales 
voluntaries pero inexpertos deben hacer parte del trabajo. A pesar de la existencia de manuales 
de evaluation con completas descripciones de los diferentes niveles de dano, la informacion que se 
maneja es altamente subjetiva y depende de la perception del evaluador en cada caso, por lo cual 
los evaluadores inexpertos tienden a calificar erroneamente los danos. Pueden cometerse errores 
como permitir la ocupacion de edificios que debian ser evacuados, o dar orden de demolition a 
estructuras que eran reparables. Estos errores llevan a serias complicaciones, especialmente en el 
caso de edificios esenciales, como los hospitales. Los niveles de dano son clasificados en los 
manuales con calificaciones linguisticas como leve, moderado, severo o fuerte, estos conceptos 
pueden tener diferentes significados dependiendo de la persona que los maneje, y no existen 
limites claramente definidos entre ellos. 

3 MODELO PARA LA EVALUACION POST-SISMICA DEL DANO 



El modelo propuesto utiliza la teoria de conjuntos difusos motivado en el caracter 
subjetivo e incompleto de la informacion. La evaluation post-sismica del dano utiliza 
calificaciones linguisticas cualitativas que pueden ser manejadas apropiadamente con la teoria de 
los conjuntos difusos. Por otro lado, se utiliza una red neuronal artificial para calibrar el sistema 
utilizando el criterio de especialistas. Esto posibilita el uso de estas herramientas para la 
evaluation del dano por inexpertos. 

Para el desarrollo del modelo se tuvieron en cuenta varios manuales de evaluation 
existentes, tales como los desarrollados en Mexico (Rodriguez, 1995), (SMIS, 1998), Japon 
(CENAPRED, 1996), Estados Unidos (ATC, 1989), Italia (Goretti, 2001), Macedonia (antes 
Yugoslavia) (IZIIS, 1985), la metodologia usada despues del sismo de 1999 en Colombia (Campos, 
1999) y la metodologia desarrollada para la ciudad de Bogota tambien en Colombia (AIS, 2002). 
Ademas miembros de la Asociacion Colombiana de Ingenieria Sismica apoyaron tecnicamente 
este trabajo. El modelo ha sido implementado en Visual Basic 6.0. El programa es de uso 
amigable al usuario, y ofrece ayudas al evaluador, tales como fotografias y descripciones de los 
tipos de danos. 

Este modelo ha sido implementado para ser una herramienta de uso oficial para la 
Alcaldia de la de ciudad de Manizales, Colombia y se realizo en el marco del proyecto “Evaluation 
de edificationes despues de un sismo”, en el cual se desarrollo un “Manual de campo para la 
inspection de edificios despues de un sismo” acompanado de un formulario de evaluation. 



3.1 Estructura de la red neuronal artificial 



La red neuronal tiene tres capas. Las neuronas en la capa de entrada son agrupadas en 4 
tipos, elementos estructurales (EE), no estructurales (NE), condiciones del suelo (CS), y 
condiciones preexistentes (CP). Cada una contribuye con informacion a una neurona intermedia, 
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la cual corresponde al grupo de variables al que pertenece. El numero de neuronas de entrada o 
de variables en el modelo, no es constante. Este numero depende del sistema estructural que se 
evalua y de la importancia de los diferentes grupos de variables seleccionados para la evaluation. 
En la tabla 1 se ve como las variables de entrada cambian de acuerdo con el sistema estructural. 
Se asigna una calificacion dependiendo del dano observado, usando cinco posibles niveles de dano 
que son tratados como conjuntos difusos. Para elementos estructurales y no estructurales los 
niveles corresponden a las siguientes calificaciones linguisticas : ninguno o muy leve (N), leve (L), 
moderado (M), fuerte (F) y severo (S) o destruction. En la figura 1 puede observarse algunos 
niveles de dano que pueden presentarse en diferentes elementos estructurales. La figura 2 ilustra 
las funciones de pertenencia para estas calificaciones linguisticas. Los conjuntos difusos son 
basados en indices de dano seleccionados, de los cuales se hablara mas adelante. 



Tabla 1 . Elementos estructurales de acuerdo con el sistema estructural 



Sistema Estructural 


Elementos Estructurales 


Porticos en concreto reforzado (con) muros de 
cortante 


Columnas/ muros, vigas, nudos y 
entrepisos 


Porticos de acero o madera 


Columnas, vigas, conexiones y entrepisos 


Mamposteria No reforzada/ Reforzada/ 
Confinada 


Muros de soporte y entrepisos 


Muros de bahareque o tapia 


Muros de soporte y entrepisos 



El dano en los elementos no estructurales no pone en peligro la estabilidad del edificio, 
pero puede representar un peligro para la seguridad de los ocupantes. Los elementos no 
estructurales estan divididos en dos grupos : elementos comunes y elementos opcionales. La tabla 
2 presenta esta clasificacion. Algunos ejemplos de los danos que pueden observarse en elementos 
no estructurales pueden observarse en la figura 3. 



Tabla 2. Elementos no Estructurales 



Elementos Comunes 


Muros divisorios 


Muros de fachada 


Escaleras 


Elementos 

Opcionales 


Cielo rasos y 
luminarias 


Instalaciones 


Cubierta 


Tanques elevados 
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c) d) 

Figura 1. Danos en elementos estructurales, a) Dano severo en un nudo de concreto reforzado, b) Dano 
moderado en una viga de concreto reforzado, c) Dano fuerte en un muro de mamposteria d) Dano fuerte en un 

muro de tapia. 

Figura 2. Funciones de pertenencia para las calificaciones linguisticas. 




a) 



b) 
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c) 



d) 



Figura 3. Danos en elementos no estructurales, a) Dano severo en un muro divisorio, b) Dano severo en 
una escalera, c) Dano moderado en un muro de fachada d) Dano severo en un tanque. 



Las variables de condiciones del suelo y preexistentes son valorados utilizando 
calificaciones linguisticas, de acuerdo con el estado observado en el momento de la evaluation. 
Las calificaciones linguisticas son: muy buena (MB), buena (B), regular (R), mala (M), y muy 
mala (MM). En las condiciones del suelo se incluyen efectos que pueden inducir una falla en la 
estabilidad de la edification que se estudia, tales como los asentamientos, subsidencias y 
licuacion del suelo y la ocurrencia de deslizamientos o fallas en el terreno; ejemplos de estas se 
pueden observar en la Figura 4. Las condiciones preexistentes se refieren a la calidad de los 
materiales de construction, irregularidades en planta y altura del edificio, y la configuration 
estructural (ejemplos en la figura 5). 

En la capa oculta o intermedia, se obtiene un indice para cada grupo de variables por 
medio de la agregacion y desfusificacion. Teniendo en cuenta los cuatro indices, es posible definir 
en la capa de salida el dano en el edificio utilizando una base de reglas difusas con las 
valoraciones del dano estructural y no estructural. La habitabilidad del edificio es obtenida 
teniendo en cuenta las condiciones del suelo, y finalmente utilizando las condiciones 
preexistentes se define el nivel de reparacion necesario para el edificio. De este modo, este 
sistema puede dar recomendaciones sobre la habitabilidad y reparabilidad del edificio. Asi con 
este sistema puede llegarse a resultados como “habitable con reparaciones menores”, “restringido: 
utilizable despues de reparacion”, o “peligroso : posible demolition”. 




a) 



b) 



Figura 4. Condiciones del suelo, a) Asentamientos y licuacion del suelo, b) Deslizamientos de tierra y 

grietas en superficie. 
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c) d) 

Figura 5. Condiciones preexistentes, a) Mala calidad de construction, b) Irregularidades en altura, piso 
blando; c) Irregularidad en planta, d) Mala configuration estructural, no hay continuidad vertical en los 

elementos. 



3.1.1 Capa de entrada 



Los conjuntos difusos para cada variable i (columnas, muros o vigas), en la capa de 
entrada, son obtenidos de las calificaciones linguisticas del dano dadas por los evaluadores Di en 
cada nivel j y su extension wj. La extension del dano (porcentaje de elementos de un tipo en cada 
nivel de dano) varia de 0 a 100 y es normalizado, Ecuacion (l). La calificacion acumulada del 
dano en cada variable Di es obtenida con la union de los conjuntos difusos escalados, teniendo en 
cuenta las funciones de pertenencia del dano pDj(Dj) y su extension o peso (asignada por el 
evaluador). La union entre conjuntos difusos esta definida por la maxima dependencia, esto es 
mas claro en las ecuaciones (2) y (3). Por medio de la desfusificacion, usando el metodo del 
centroide de area (COA), se obtiene un indice de calificacion C, para cada variable de cada grupo 
de neuronas. En la Ecuacion (4) describe el metodo de desfusificacion utilizado. En la Figura 6 se 
ilustra el ingreso de datos al sistema. 



D 



w i = 



I* 



] —,Yw.=i 



D i = (D N u D L D M u D f kj D s ) 
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Pd, (D) = max (vy V i * p Dk (d n . w S i *p 0j [D Si )) 
C, = [max (vv W i * p„ v [d nj w S l * {d S i ))\ centmide 




Figura 6. Ingreso de la extension del dano de cada nivel para un elemento estructural. 



Cada variable tiene predefinidas las funciones de pertenencia basicas de los conjuntos 
difusos correspondientes a los cinco niveles de dano posibles. Las calificaciones linguisticas 
cambian en cada caso como se menciono antes. 



3.1.2 Capa intermedia 



En esta capa se tienen cuatro neuronas, que corresponden a cada uno de los grupos de 
variables: elementos estructurales, elementos no estructurales, condiciones del suelo y 
condiciones preexistentes. En la figura 7 se presenta un esquema general del proceso de 
evaluacion. 

En este modelo de red neuronal las neuronas de entrada no tienen influencia sobre todas 
las neuronas intermedias, solo lo hacen sobre la neurona que corresponde al grupo de variables al 
que pertenecen. Los parametros de entrada de estas 4 neuronas son las calificaciones C, 
obtenidas para cada variable de cada grupo de neuronas y su peso Wi, o el grado de importancia 
dados por el evaluador de acuerdo con su criterio. Estos pesos son normalizados y calibrados por 
medio de una funcion de aprendizaje de la que se hablara mas adelante. Para el proceso de 
entrenamiento, los valores iniciales han sido definidos teniendo en cuenta el criterio de expertos. 
Con estas calificaciones y pesos de cada variable i se obtiene un indice global para cada grupo k, 
por medio de la union y desfusificacion de los conjuntos difusos escalados. Esto se hace con las 
funciones de pertenencia //«-, (Cki) y los pesos ILw como se indica en las ecuaciones (5) y (6). 



Pc££ (C) 

I C’C’ 



max {W EE1 * p c EE , (C EE1 )>•••> ^ EEi * p Cff . ip EEi )) 

= [max ( W EE j - P Ceei {C eei )>-> W eei * M- c EEi i^EEi ))\ centroide 
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Las variables de condiciones del suelo y preexistentes pueden ser tomadas o no en cuenta 
dentro de la evaluacion, dependiendo de su importancia en los resultados, esto hace que en 
algunos casos la habitabilidad y reparabilidad puedan ser obtenidas unicamente a partir del dano 
estructural y no estructural. 

Figura 7. Estructura de la red neuronal propuesta 



3.1.3 Capa de salida 



En esta capa se obtienen indices globales para los elementos estructurales, no 
estructurales, condiciones del suelo y preexistentes, corresponden a una calificacion linguistica 
final en cada caso. El nivel de dano es obtenido de acuerdo con la proximidad del valor del indice 
obtenido a la funcion de dano de referenda. Los niveles para los elementos estructurales y no 
estructurales son definidos de igual forma que antes, como^ ninguno, leve, moderado, fuerte y 
severo. Tambien para las condiciones preexistentes y del suelo se definen asi : muy buenas, 
buenas, regulares, malas y muy malas. En esta capa tiene lugar el proceso de entrenamiento, los 
indices que identifican cada nivel cualitativo (centro de cluster) cambian de acuerdo con los 
indices calculados en cada evaluacion que se hace con la red, teniendo en cuenta una tasa de 
aprendizaje. Esta ultima determina la velocidad de entrenamiento de la red. El procedimiento de 
aprendizaje es descrito en el siguiente apartado. Despues de hallar las calificaciones finales 
utilizando bases de reglas difusas se determinan el dano global en el edificio, su habitabilidad y 
reparabilidad. En secciones posteriores se explicara este procedimiento. 

3.2 Aprendizaje de la red neuronal 



La red neuronal es calibrada en la capa de salida, cuando se definen las funciones de dano 
en relation la matriz de indices de dano. Para dar inicio a este proceso se le dio un punto de 
partida a la red, se definio un valor initial para cada nivel de dano. Para la determination de 
estos indices iniciales han sido considerados los propuestos por el ATC-13 (ATC, 1985), Park, Ang 
y Wen (Park et al., 1984), las curvas de fragilidad usadas por Hazus-99 (FEMA, 1999) y los 
indices utilizados por Sanchez-Silva (Sanchez-Silva et al. 2001). Los valores de los indices 
iniciales corresponden a los centroides de area de cada funcion de pertenencia a cada nivel de 
dano. La tabla 3 presenta una comparacion entre los indices aqui propuestos, los de Park, Ang y 
Wen, y los de Sanchez-Silva y Garcia. 



Tabla 3. Tabla comparativa de los indices de dano 



Nivel de dano 


Park, Ang y Wen 


Sanchez Garcia 


Propuestos 


Muy leve 


<0.1 

0.07 


0.10 


0.07 


Leve 


0.10-0.25 

0.175 


0.20 


0.17 


Moderado 


0.25 - 0.40 
0.325 


0.35 


0.33 


Severo 


0.40-0.80 

0.6 


0.60 


0.55 


Destruction 


>0.80 

0.8 


0.90 


0.76 



Los indices iniciales para la matriz de dano estan basados en los indices de Park, Ang y 
Wen, dado que estos han sido calibrados con diferentes estudios y danos de terremotos reales y 
son los que cuentan con mayor aceptacion. Sus autores consideran la existencia de colapso desde 
0.8, mientras que Stone y Taylor (Stone et al. 1993) lo proponen a partir de 0.77. Considerando 





Estimation de la habitabilidad y reparabilidad de edificios despues de un sismo 



103 



esto se tomo 0.76 como valor para el indice del nivel de destruction o colapso. En la election de los 
indices de los niveles severo y moderado se ha intentado ser conservativos, dado la intensa 
discusion sobre estos para reducir su valor. 



La calibration se hace por nivel de dano y solo para los indices correspondientes a los grupos 
de variables considerados en cada evaluation que calibra el sistema. En la capa de salida el 
aprendizaje utiliza el algoritmo de Kohonen, como se ve en la ecuacion (7). 

hi (t+ ] )= hj (0 + - hj J 



donde I k j es el valor del mdice de un grupo de valuables k que presenta un nivel j, que se recalcula 
considerando una tasa de aprendizaje a y la diferencia entre el mdice resultante de la presente 
evaluacion y el mdice previo en cada nivel de dano j. En la ecuacion (8) se define la tasa de 
aprendizaje de la ecuacion (7). 



a (t) = 0.1* Exp (-0.1* t) 



donde t es el nurnero de veces que el indice o peso que se calibra ha sido usado. Para el 
entrenamiento de la red se utilizaron las evaluaciones hechas despues del terremoto del Quindio 
en Colombia en 1999. En el entrenamiento no se tuvieron en cuenta sistemas estructurales como 
porticos en madera y acero, esto debido a la ausencia de estos sistemas estructurales en la zona 
afectada por el sismo. Se utilizaron pocos edificios con porticos en concreto reforzado con muros 
de cortante, esto debido al reducido nurnero de evaluaciones de estructuras de este tipo. 



3.3 Bases de reglas fuzzy 



La habitabilidad y reparabilidad del edificio que se evalua son estimadas basandose en el 
nivel de dano de los elementos estructurales y no estructurales, el estado del suelo y las 
condiciones preexistentes. La Figura 8 ilustra el uso de reglas difusas para estas estimaciones. El 
nivel de dano global del edificio es deducido a partir de los resultados de dano estructural y no 
estructural; y tiene cinco posibles calificaciones : ninguno, leve, moderado, fuerte y severo. El 
estado global del edificio y la habitabilidad son estimados teniendo en cuenta el dano global del 
edificio y las condiciones del suelo en dos diferentes bases de reglas fuzzy. La habitabilidad puede 
ser clasificada en cuatro niveles : habitable, uso restringido, prohibida la entrada y peligroso. La 
reparabilidad depende de otra base de reglas que contempla las condiciones preexistentes de la 
edification, aqui se contribuye a definir la factibilidad tecnica y economica de una reparacion. La 
reparabilidad tiene cuatro posibles niveles : ninguna o reparaciones menores, reparacion, 
refuerzo, y posible demolition, estos niveles son identificados usualmente con colores en los 
manuales de evaluacion, asi : verde, amarillo, naranja y rojo. Para cada edificio que se evalua con 
este sistema es posible ayudar en dos decisiones criticas que son basicas en la fase de respuesta a 
la emergencia despues de ocurrido un sismo fuerte. En la figura 8 se ilustra la presentation de los 
resultados dada por el programa. 



4 CONCLUSIONES 



Despues de una completa revision de los diferentes manuales para evaluacion post-sismica de 
dano, se propone un sistema experto innovador que apoye el proceso de evaluacion. Para el 
desarrollo de esta herramienta se han tenido en cuenta las ventajas y desventajas de cada uno 
de los metodos. 
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El sistema experto fue desarrollado utilizando tecnicas de inteligencia artificial como las redes 
neuronales artificiales y la teoria de conjuntos difusos. Este tipo de herramientas es apropiado 
en esta practica, debido al caracter subjetivo de la evaluacion de los danos y que la 
information suele ser incompleta. 

Evaluaciones hechas por ingenieros expertos despues del terremoto del Quindio, en Colombia, 
en 1999, han sido utilizadas para el entrenamiento del sistema. 



Figura 8. Procedimiento para determinar la Habitabilidad y Reparabilidad 



Se confirmo la necesidad de tener un metodo preestablecido para la evaluacion del dano 
en edificios en caso de un sismo fuerte en ciudades ubicadas en zonas de sismicidad 
considerable. 



El uso de herramientas de inteligencia artificial en la ingenieria civil ha tenido muy poca 
difusion, es recomendable promover su uso para dar soluciones adecuadas y rapidas a 
problemas en este campo del conocimiento. 




a) 



b) 
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c) 



d) 



Figura 9. Presentation de los resultados de una evaluation de danos realizada con el sistema “EDE - 

Manizales” 



5. AGRADECIMIENTOS 



Estimation de la habitabilidad y reparabilidad de edificios despues de un sismo 



105 



Este trabajo ha sido parcialmente patrocinado por el Ministerio de Ciencia y Tecnologla 
de Espana, N° de referenda del proyecto REN2002-03365/RIES y por la Comision Europea, 
Proyecto Risk-UE, contrato EVK4-CT-2000-00014, y el programa de ayudas academicas del 
ProVention Consortium del Banco Mundial. 



REFERENCIAS 



1. - AIS, Asociacion Colombiana de Ingenierla Slsmica (2002). Gula Tecnica para Inspection de 

Edificaciones Despues de un Sismo”. Manual de Campo. Bogota, D.C., Colombia. 

2. - Applied Technology Council (1989). Procedures for postearthquake safety evaluation of buildings, 

ATC-20. Redwood City, CA. 

3. - Applied Technology Council, (1985). Earthquake damage evaluation data for California, ATC-13. 

Redwood City, CA. 

4. - Campos, A., (1999). Memoria Tecnica del Censo de Inmuebles Afectados por el Sismo del 25 de 

Enero de 1999 en El Eje Cafetero. Ministerio de Desarrollo Economico. In Spanish.. Bogota D.C., 
Colombia. 

5. - Centro Nacional de Prevention de Desastres (CENAPRED) (1996). In Spanish. Mexico. Norma 

para la Evaluation del Nivel de Dano por Sismo en Estructuras y Guia Tecnica de Rehabilitation 
(Estructuras de Concreto Reforzado). Cuadernos de Investigation, Numero 37. 

6. - FEMA, Federal Emergency Management Agency, (1999). “Earthquake Loss Estimation 

Methodology HAZUS”, Technical Manual, Vol. I, II and III, First edition 1997, National Institute 
of Buildings Sciences of Federal Emergency Management Agency, Washington. 

7. - Goretti A., (2001). Post- earthquake building usability: a Assessment. Technical Report 

SSN/RT/01/03. 

8. - Instituto de Ingenieria Sismica y Sismologia -IZIIS. (1984). In Spanish. Metodologia y 

Procedimiento para la Evaluation de Danos Producidos por Terremotos. Universidad “Kiril y 
Metodij” - Skopje, Yugoslavia. 

9. - Park Y.J., Ang A. and Wen Y. (1984). Seismic Damage Analysis and Damage-Limiting Design of 

R.C. Buildings. Structural Research Series, Report No 516, University of Illinois at Urban- 
Champaign, Urbana 

10. -Rodriguez, M. y Castrillon, E. (1995). In Spanish. Instituto Nacional de Ingenieria UNAM. 

Manual de Evaluation Postsismica de la Seguridad Estructural de Edificaciones. Basado en 
investigaciones Realizadas para el Departamento del Distrito Federal. Series del Instituto de 
Ingenieria 569. 

11. - Sanchez-Silva, M. and Garcia, L. Earthquake Damage Assessment Based on Fuzzy Logic and 

Neural Networks. EERI Earthquake Spectra, Vol. 17, N. 1, February, pp. 89-112. Oakland, 
California, 

12. -Sociedad Mexicana de Ingenieria Sismica (1998). A.C. Secretaria de Obras y Servicios Gobierno 

del Distrito Federal. In Spanish. Manual de Evaluation Postsismica de la Seguridad Estructural 
de Edificaciones. 

13. - Stone, W.C., Taylor A.W. (1993). Seismic Performance of Circular Bridge Columns Designed in 

Accordance with AASHTO/CALTRANS standards, NIST Building Science Series 170, National 
Institute of Standards and Technology, Gaithersburg MD 




Revista Internacional de 



Vol. 10-1, 109-118(2005) 



Ingenierfa 
de estructuras 



Evaluacion rapida de la Deriva Maxima de piso en 
Edificios de Hormigon Armado 

Roberto Aguiar 

ESPE - Centro de Investigaciones Cientfficas - 
Av. El Progreso s/n Valle de los Chillos - 
Quito - Ecuador - raguiar@espe.edu.ec 

Carlos Bobadilla 

ESPE - Centro de Investigaciones Cientfficas - 
Av. El Progreso s/n Valle de los Chillos - 
Quito - Ecuador - lioncarl007 @ hotmail.com 



RESUMEN 

Se comparan las distorsiones maximas de piso, obtenidos del analisis de 72 
edificios de hormigon armado compuestos por vigas y columnas ante la accion de 9 
registros sfsmicos cuya aceleracion maxima del suelo es mayor a 0.1 g., al utilizar un 
analisis no lineal paso a paso, con las que se encuentran al aplicar una formula empfrica 
que considera cinco parametros. 

Uno de los parametros de la formula empfrica ha sido obtenido del analisis de 63 
acelerogramas registrados en Colombia, Peru, Chile y Argentina y otro de los 
parametros toma en cuenta el deterioro de rigidez, deterioro de resistencia y efecto de 
cierre de grietas en el modelo constitutive. 



ABSTRACT 

The maximum story lateral drifts obtained from analyses of 72 reinforced 
concrete frame buildings, on the action of 9 seismic registers which maximum ground 
acceleration is higher than O.lg using a nonlinear step by step analysis, are compared 
which those determinated using an empirical formulae which considered 5 parameters. 

One of the empirical formulae’s parameter has been obtained from analysis of 
63 accelerographs registered in Colombia, Peru, Chile and Argentina, other parameters 
considerate stiffness, strength degradation and the pinching effect. 
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1 INTRODUCTION 

En la ultima decada del siglo pasado se dio gran importancia al desarrollo de 
metodos simplificados que reportan el punto de desempeno de una estructura ante una 
determinada accion sfsmica, punto que esta determinado por un desplazamiento lateral 
maximo o una deriva de piso maxima. Para el efecto se desarrollaron varios algoritmos 
para el Metodo del Espectro de Capacidad que fueron recomendados por el ATC-40 
(1996), los mismos que fueron criticados debido a que en algunos casos, estos 
algoritmos no eran convergentes. 

Varias son las contribuciones realizadas para determinar el punto de 
desempeno empleando el Metodo del Espectro de Capacidad, entre las que se puede 
citar a las formuladas por Bertero (1995), Reinhorn (1997), Chopra y Goel (2002), entre 
otras. Lo cierto es que en todas estas metodologfas se debe suministrar, a un programa 
de ordenador, la geometrfa de la estructura, el armado de los elementos y la calidad de 
los materiales. Informacion que demanda un tiempo considerable, pero es un poco 
menor a la que se suministra cuando se utilizan programas de analisis no lineal, paso a 
paso. 



Con el proposito de encontrar el punto de desempeno en una forma mas rapida 
se han propuesto ecuaciones empfricas, las mismas que estan orientadas, 
fundamentalmente a la evaluacion de la vulnerabilidad sfsmica de un conjunto de 
estructuras, las mismas que demandarfa demasiado tiempo si se pretende analizar con el 
Metodo del Espectro de Capacidad o peor aun con los algoritmos de analisis no lineal, 
paso a paso. 

Una de ecuaciones empfricas es la propuesta por Miranda (1997) y sirve para 
encontrar la deriva maxima de piso ^ , empleando la siguiente ecuacion. 



AAAA 

' TT 



( 1 ) 



Donde ^ es el valor de paso del sistema de un grado de libertad al sistema de 

multiples grados de libertad; es un factor de amplificacion que permite determinar la 
distorsion maxima de entrepiso a partir de la distorsion global de la estructura. La 
distorsion global se define como la relacion entre el desplazamiento lateral maximo en 

el tope del edificio dividido para la altura total del edificio; es un factor que permite 
calcular los desplazamientos laterales maximos con comportamiento inelastico a partir 

de los maximos desplazamientos laterales con comportamiento elastico ; es un factor 
que sirve para determinar el cociente entre la distorsion maxima de entrepiso y la 
distorsion global pero calculado en una estructura con comportamiento inelastico con 
relacion a la misma relacion pero calculada con comportamiento elastico; H e s la altura 
5 

total del edificio y d es el desplazamiento espectral elastico asociado al perfodo 
T 

efectivo c de la estructura. 
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Una vez que se determina la deriva maxima de piso, ante una accion sfsmica 
dada se puede inferir el desempeno de la misma, utilizando los valores recomendados 
por Ghobarah (2004) o por el comite VISION 2000 (1995). 

2 MODIFICACION A LA ECUACION EMPIRICA 

Los modelos constitutivos de los materiales que han sido utilizados por la 
mayor parte de investigadores (Miranda 2000, Whittaker et al 1998, Nassar y 
Krawinkler 1991, Shimazaki y Sozen 1984, entre otros) para determinar alguno de los 
parametros indicados en la ecuacion (1) no han contemplado el deterioro de rigidez en 
la descarga, la perdida de resistencia y el efecto de cierre de grietas. En efecto la mayor 
parte de ellos consideraron un modelo elasto perfectamente plastico o un modelo 
bilineal que contempla rigidez post fluencia. 



Por este motivo es que se incorporo el parametro , Aguiar (2005) que toma 
en cuenta la perdida de rigidez, de resistencia y efecto pinching, quedando la ecuacion 
(1) de la siguiente manera: 



_ A Pi P± Pa 

' TT 



( 2 ) 



Donde se obtiene en funcion de la demanda de ductilidad ^ y los valores 

se indican en la tabla 1. Se destaca que FEMA 273 recomienda un solo valor de y es 
de 1.2 



Tabla 1 Valores de en funcion de la demanda de ductilidad. 



Ductilidad 


1 


2 


3 


4 


5 


6 


A 


1.00 


1.14 


1.17 


1.19 


1.22 


1.23 



La ecuacion (2) es la que se utilizara para encontrar la distorsion maxima de 
piso, bajo el tltulo de metodologfa rapida y cuyos valores seran comparados con los que 
reporta el programa ID ARC al realizar un analisis no lineal, paso a paso. 

3 PARAMETRO & 

Se han propuesto ecuaciones para el calculo de ^ 3 por parte de Miranda 
(2000), Nassar y Krawinkler (1991), entre otros, las mismas que fueron obtenidas en 
base a sismos registrados en los Estados Unidos de Norte America, razon por la cual 
Aguiar y Guerrero (2005) en base al analisis de 63 acelerogramas registrados en 
Colombia, Peru, Chile y Argentina, recomiendan las siguientes ecuaciones para calcular 

, que son las que mejor se ajustan a los resultados obtenidos y que se indican en las 
figuras 1 y 2. 
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P, 



H/'-o+r 



c(T,a) 



T e 2 - 07 0,381 

l + T e zm+ T e 



c(T,a ) 



T e 1247 ^ 0.248 

i+r/- 247 + T e 



para a = 0.0 



para a = 0.05 



(3) 

(4) 

(5) 



Siendo a la relacion entre la rigidez post fluencia con respecto a la rigidez 
elastica de la curva de capacidad sfsmica resistente que relaciona el cortante basal con el 
desplazamiento lateral maximo. Siempre se debera considerar a mayor que 1. 




PERiODO (s) 

— DUCnUDAD 2 DUCTILIDAD 3 DUCTILIDAD 4| 



Figura 1 Valores medios de para a ~ 

4 ESTRUCTURAS Y SISMOS DE ANALISIS 

En Aguiar et al (2005) se describe con detalle la geometrfa de los 72 edificios 
de hormigon armado de 1 a 6 pisos que fueron analizados. La distribucion en planta de 
todos los edificios es la indicada en la figura 3. La cuantfa de la armadura longitudinal 
de columnas vario entre el 1 % y el 2 %; y, en las vigas la cuantfa vario del 0.5 % al 1.5 

%. El refuerzo transversal para vigas y columnas esta constituido por 1 ^ de 8 mm., de 
diametro, espaciado a 10 cm., en los extremos del elemento y a 20 cm., en el centra de 
luz. 
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La dimension de las columnas variaron desde 20/20 cm., en las construcciones 
de 1 piso hasta 45/45 en los edificios de 6 pisos. Las vigas van desde 20/20 cm., en las 
de un piso hasta 60/30 en las de 6 pisos; la primera dimension representa la base y la 
segunda la altura de la seccion transversal. La altura de cada entrepiso es constante y 
vale 3.0 m. Las dimensiones y armadura de los edificios considerados corresponde a la 
forma como se construye en el Ecuador. 




0.00000 0.50000 1.00000 1.50000 2.00000 2.50000 3.00000 

PERIODO (s) 

— DUCTILIDAD 2 DUCTILIDAD 3 DUCTILIDAD 4| 

Figura 2 Valores medios de para a = 0-05 




Figura 3 Planta tipo de edificios analizados. 
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Tabla 2 Datos relev antes de los sismos considerados en el estudio. 



Reg. 


Estacion 


Pais 


Fecha 


Magnitud 


Comp. 


Distancia 

Epic. 

(km.) 


Aceler. 

Maxima 

(Gal.) 


1 


Armenia (Carme) 


Colombia 


25/01/99 


6.20 


E-W 


13 


518,4 


2 


Armenia (Carme) 


Colombia 


25/01/99 


6.20 


N-S 


13 


-580.1 


3 


Finlandia (C Flan) 


Colombia 


25/01/99 


6.20 


E-W 


33 


554.9 


4 


Finlandia (C Flan) 


Colombia 


25/01/99 


6.20 


N-S 


33 


-477.8 


5 


Pereira (C Per) 


Colombia 


25/01/99 


6.20 


E-W 


48 


-207.8 


6 


Pereira (C Per) 


Colombia 


25/01/99 


6.20 


N-S 


48 


141.5 


7 


Ricaurte (C Rica) 


Colombia 


1 7/07/98 


3.07 


E-W 


637 


132.4 


8 


Ricaurte (C Rica) 


Colombia 


1 7/07/98 


3.07 


N-S 


637 


-124.8 


9 


Pereira (C Mazpa) 


Colombia 


1 9/02/97 


6.10 


E-W 


154 


127.7 



La carga vertical vario desde 400 kg/m2 en el de un piso hasta 700 kg/m2 en el 
de 6 pisos. Finalmente el hormigon utilizado fue de 210 kg/cm2 y el acero una fluencia 
de 4200 kg/cm2. 



En la tabla 2 se indican los datos mas relevantes de los sismos considerados en 
el estudio y estos son: el codigo del registro, la estacion en la cual fue registrada, el 
Pais, la fecha, la magnitud, la componente del registro, la distancia epicentral y la 
aceleracion maxima del registro. Todos los eventos fueron registrados en Colombia y la 
mayor parte de ellos corresponden al sismo que afecto al eje cafetero el 25 de enero de 
1999. 

5 ANALISIS NO LINEAL Y EVALUACION DEL DRIFT 

Se analizaron 648 casos que corresponden a las 72 estructuras sometidas a los 
9 registros indicados en la tabla 2. El analisis se realizo de dos maneras, la primera 
utilizando el programa ID ARC y la segunda empleando la ecuacion (2). 



Para el primer caso se encontro la respuesta no lineal, paso a paso, de un 
portico central del edificio, por medio del programa ID ARC, de este estudio, interesaba 
ver unicamente la distorsion maxima de piso y la demanda de ductilidad la misma que 
en la mayorfa de los casos fue superior a 4. 



Para la evaluacion con la ecuacion (2); el parametro fue obtenido con un 
programa de analisis sfsmico elastico, este parametro no es mas que el factor de 
participacion del primer modo multiplicado por el valor modal correspondiente al 

ultimo piso. El parametro que relaciona la deriva maxima de piso con la deriva 
global del edificio fue obtenido de los resultados del programa ID ARC. El parametro 

se hallo con las ecuaciones (3) y (5) teniendo en cuenta que este valor sea mayor a 

la unidad. El parametro @ 4 se encontro con la siguiente ecuacion. 



/?4 — 1 + + 

30 



N 

200 



( 6 ) 
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Donde N es e ] numero de pisos. El parametro se determino en base a los 
valores indicados en la tabla 1. Para la determinacion de los valores ^ 

considero ^ ~ ^ . 



se 



T 

Por otra parte, el perfodo efectivo e se obtuvo del espectro de capacidad de 

5 

cada una de las estructuras, Gomez y Zambrano (2005). Finalmente, para el valor de d 
de los espectros elasticos de desplazamiento que se indican en la figura 4. 

En Bovadilla (2005) se indica con detalle los valores obtenidos del programa 
ID ARC y los valores de cada uno de los parametros que intervienen en la ecuacion (2) 
para cada uno de los 648 casos. 




Sismo 1 



Sismo 2 



Sismo 3 



Espectro de Respuesta al Desplazamiento 


1 l5 ° E ’ 01 
5) 1.00E+01 






































0.2 0.4 0 


Periodo (T) 



Sismo 4 



Espectro de Respuesta al Desplazamiento 




Sismo 5 



Espectro de Respuesta al Desplazamiento 




Sismo 6 




Sismo 7 Sismo 8 Sismo 9 



Figura 4 Espectros de desplazamiento elasticos utilizados en el estudio. 
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6 RESULTADOS 

Se denomina y "> a la deriva maxima de piso obtenida con el programa ID ARC 
y Y a la deriva maxima de piso encontrada con la ecuacion (2). Para poder visualizar 

los resultados en una sola figura se define la relacion Y "■> ^ y . S i esta relacion tiende a 1 
significa que la aproximacion de la deriva maxima de piso que se encuentra en forma 
rapida con la ecuacion (2) reporta muy buenos resultados. En la figura 5 se presenta esta 
relacion para todos los 648 casos analizados, se aprecia que la mayor parte de valores 
tienen un valor inferior a 2 pero existen casos aislados en que esta relacion es alta. 



a 



s 

'2 

SJ 

3 

Su 

SC 



10,000 

9.000 

8.000 

7.000 

6.000 

5.000 

4.000 

3.000 

2.000 
1,000 
0,000 




Piso 



-PISO 1 
PISO 2 
PISO 3 
PISO 4 
-PISO S 
-PISO 6 
-MEDIA 



Figura 5 Relacion entre la deriva maxima de piso encontrada con ID ARC y la evaluation rapida. 



La lmea continua de la figura 5, que une a todos los casos, corresponde al valor 
medio encontrado con los casos de cada piso. Esta lmea tiene un valor promedio que 
esta alrededor de 2 que se considera satisfactorio si se toma en cuenta el tiempo de 
calculo que demanda encontrar la deriva maxima de piso con la evaluation rapida es 
pequeno en comparacion con lo que demanda la preparation de datos para el programa 
ID ARC. En la figura 6 se muestra la desviacion estandar encontrada para cada piso, en 
forma general se aprecia que estan alrededor de 2. 




Piso 



Figura 6 Desviacion estandar encontrada para cada piso. 
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7 CON CLU SIONES 

Se ha comparado la deriva maxima de piso que se obtiene al utilizar el 
programa IDARC al encontrar la respuesta no lineal, paso a paso, con la que se 
encuentra al aplicar la metodologfa rapida, la misma que ha sido modificada con la 

incorporacion del parametro y con la determinacion del parametro que se hallo 
en base a sismos registrados en Colombia, Peru, Chile y Argentina. 

De la comparacion de los 648 casos se desprende que la metodologfa rapida, 
con los parametros definidos en el estudio, reporta resultados que se consideran 
satisfactorios, tomando en cuenta que con una ecuacion tan sencilla se obtiene el punto 
de desempeno que al aplicar otros metodos obligatoriamente se debera recurrir al uso de 
un programa de ordenador con la consiguiente demanda de tiempo. 
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